• Chaotic dynamics in a macrospin spin-torque nano-oscillator with delayed feedback(1709.04310)

March 11, 2019 nlin.CD, cond-mat.mes-hall
A theoretical study of delayed feedback in spin-torque nano-oscillators is presented. A macrospin geometry is considered, where self-sustained oscillations are made possible by spin transfer torques associated with spin currents flowing perpendicular to the film plane. By tuning the delay and amplification of the self-injected signal, we identify dynamical regimes in this system such as chaos, switching between precession modes with complex transients, and oscillator death. Such delayed feedback schemes open up a new field of exploration for such oscillators, where the complex transient states might find important applications in information processing.
• Chaotic Properties of Single Element Nonlinear Chimney Model: Effect of Directionality(1708.01092)

Feb. 21, 2019 nlin.CD
We generalize the chimney model by introducing nonlinear restoring and gravitational forces for the purpose of modeling swaying of trees at high wind speeds. We have derived general equations governing the system using Lagrangian formulation. We have studied the simplest case of a single element in more detail. The governing equation we arrive at for this case has not been studied so far. We study the chaotic properties of this simple building block and also the effect of directionality in the wind on the chaotic properties. We also consider the special case of two elements.
• Trajectory-free approximation of phase space structures using the trajectory divergence rate(1705.07949)

Feb. 9, 2019 nlin.CD
This paper introduces the trajectory divergence rate, a scalar field which locally gives the instantaneous attraction or repulsion of adjacent trajectories. This scalar field may be used to find highly attracting or repelling invariant manifolds, such as slow manifolds, to rapidly approximating hyperbolic Lagrangian coherent structures, or to provide the local stability of invariant manifolds. This work presents the derivation of the trajectory divergence rate and the related trajectory divergence ratio for 2-dimensional systems, investigates their properties, shows their application to several example systems, and presents their extension to higher dimensions.
• Edgeworth expansions for slow-fast systems with finite time scale separation(1708.06984)

Feb. 4, 2019 nlin.CD, math.PR, math.DS
We derive Edgeworth expansions that describe corrections to the Gaussian limiting behaviour of slow-fast systems. The Edgeworth expansion is achieved using a semi-group formalism for the transfer operator, where a Duhamel-Dyson series is used to asymptotically determine the corrections at any desired order of the time scale parameter $\varepsilon$. The corrections involve integrals over higher-order auto-correlation functions. We develop a diagrammatic representation of the series to control the combinatorial wealth of the asymptotic expansion in $\varepsilon$ and provide explicit expressions for the first two orders. At a formal level, the expressions derived are valid in the case when the fast dynamics is stochastic as well as when the fast dynamics is entirely deterministic. We corroborate our analytical results with numerical simulations and show that our method provides an improvement on the classical homogenization limit which is restricted to the limit of infinite time scale separation.
• A Tube Dynamics Perspective Governing Stability Transitions: An Example Based on Snap-through Buckling(1705.05903)

Feb. 1, 2019 nlin.CD
The equilibrium configuration of an engineering structure, able to withstand a certain loading condition, is usually associated with a local minimum of the underlying potential energy. However, in the nonlinear context, there may be other equilibria present, and this brings with it the possibility of a transition to an alternative (remote) minimum. That is, given a sufficient disturbance, the structure might buckle, perhaps suddenly, to another shape. This paper considers the dynamic mechanisms under which such transitions (typically via saddle points) occur. A two-mode Hamiltonian is developed for a shallow arch/buckled beam. The resulting form of the potential energy---two stable wells connected by rank-1 saddle points---shows an analogy with resonance transitions in celestial mechanics or molecular reconfigurations in chemistry, whereas here the transition corresponds to switching between two stable structural configurations. Then, from Hamilton's equations, the analytical equilibria are determined and linearization of the equations of motion about the saddle is obtained. After computing the eigenvalues and eigenvectors of the coefficient matrix associated with the linearization, a symplectic transformation is given which puts the Hamiltonian into normal form and simplifies the equations, allowing us to use the conceptual framework known as tube dynamics. The flow in the equilibrium region of phase space as well as the invariant manifold tubes in position space are discussed. Also, we account for the addition of damping in the tube dynamics framework, which leads to a richer set of behaviors in transition dynamics than previously explored.
• Isotropic polarization of compressible flows(1509.02988)

Jan. 28, 2019 physics.flu-dyn, nlin.CD
The helical absolute equilibrium of a compressible adiabatic flow presents not only the polarization between the two purely helical modes of opposite chiralities but also that between the vortical and acoustic modes, deviating from the equipartition predicted by {\sc Kraichnan, R. H.} [1955 The Journal of the Acoustical Society of America {\bf 27}, 438--441.]. Due to the existence of the acoustic mode, even if all Fourier modes of one chiral sector in the sharpened Helmholtz decomposition [{\sc Moses, H. E.} 1971 SIAM ~(Soc. Ind. Appl. Math.) J. Appl. Math. {\bf 21}, 114--130] are thoroughly truncated, leaving the system with positive definite helicity and energy, negative temperature and the corresponding large-scale concentration of vortical modes are not allowed, unlike the incompressible case.
• Backward Simulation of Stochastic Process using a Time Reverse Monte Carlo method(1708.08045)

The "backward simulation" of a stochastic process is defined as the stochastic dynamics that trace a time-reversed path from the target region to the initial configuration. If the probabilities calculated by the original simulation are easily restored from those obtained by backward dynamics, we can use it as a computational tool. It is shown that the naive approach to backward simulation does not work as expected. As a remedy, the Time Reverse Monte Carlo method (TRMC) based on the ideas of Sequential Importance Sampling (SIS) and Sequential Monte Carlo (SMC) is proposed and successfully tested with a stochastic typhoon model and the Lorenz 96 model. TRMC with SMC, which contains resampling steps, is shown to be more efficient for simulations with a larger number of time steps. A limitation of TRMC and its relation to the Bayes formula are also discussed.
• Weyl Spreading Sequence Optimizing CDMA(1602.04584)

Jan. 18, 2019 nlin.CD, cs.IT, math.IT
This paper shows an optimal spreading sequence in the Weyl sequence class, which is similar to the set of the Oppermann sequences for asynchronous CDMA systems. Sequences in Weyl sequence class have the desired property that the order of cross-correlation is low. Therefore, sequences in the Weyl sequence class are expected to minimize the inter-symbol interference. We evaluate the upper bound of cross-correlation and odd cross-correlation of spreading sequences in the Weyl sequence class and construct the optimization problem: minimize the upper bound of the absolute values of cross-correlation and odd cross-correlation. Since our optimization problem is convex, we can derive the optimal spreading sequences as the global solution of the problem. We show their signal to interference plus noise ratio (SINR) in a special case. From this result, we propose how the initial elements are assigned, that is, how spreading sequences are assigned to each users. In an asynchronous CDMA system, we also numerically compare our spreading sequences with other ones, the Gold codes, the Oppermann sequences, the optimal Chebyshev spreading sequences and the SP sequences in Bit Error Rate. Our spreading sequence, which yields the global solution, has the highest performance among the other spreading sequences tested.
• A new approach to fuzzy sets: Application to the design of nonlinear time-series, symmetry-breaking patterns, and non-sinusoidal limit-cycle oscillations(1704.00676)

Jan. 17, 2019 math-ph, math.MP, nlin.CD
It is shown that characteristic functions of sets can be made fuzzy by means of the $\mathcal{B}_{\kappa}$-function, recently introduced by the author, where the fuzziness parameter $\kappa \in \mathbb{R}$ controls how much a fuzzy set deviates from the crisp set obtained in the limit $\kappa \to 0$. As applications, we present first a general expression for a switching function that may be of interest in electrical engineering and in the design of nonlinear time-series. We then introduce another general expression that allows wallpaper and frieze patterns for every possible planar symmetry group (besides patterns typical of quasicrystals) to be designed. We show how the fuzziness parameter $\kappa$ plays an analogous role to temperature in physical applications and may be used to break the symmetry of spatial patterns. As a further, important application, we establish a theorem on the shaping of limit cycle oscillations far from bifurcations in smooth deterministic nonlinear dynamical systems governed by differential equations. Following this application, we briefly discuss a generalization of the Stuart-Landau equation to non-sinusoidal oscillators obtained as a consequence of our theorem.
• Spontaneous Symmetry Breaking, Conformal Anomaly and Incompressible Fluid Turbulence(1707.07855)

Dec. 29, 2018 hep-th, nlin.CD
We propose an effective conformal field theory (CFT) description of steady state incompressible fluid turbulence at the inertial range of scales in any number of spatial dimensions. We derive a KPZ-type equation for the anomalous scaling of the longitudinal velocity structure functions and relate the intermittency parameter to the boundary Euler (A-type) conformal anomaly coefficient. The proposed theory consists of a mean field CFT that exhibits Kolmogorov linear scaling (K41 theory) coupled to a dilaton. The dilaton is a Nambu-Goldstone gapless mode that arises from a spontaneous breaking due to the energy flux of the separate scale and time symmetries of the inviscid Navier-Stokes equations to a K41 scaling with a dynamical exponent $z=\frac{2}{3}$. The dilaton acts as a random measure that dresses the K41 theory and introduces intermittency. We discuss the two, three and large number of space dimensions cases and how entanglement entropy can be used to characterize the intermittency strength.
• Waves along fractal coastlines: From fractal arithmetic to wave equations(1707.06225)

Dec. 25, 2018 math-ph, math.MP, nlin.CD, nlin.PS, math.DS
Beginning with addition and multiplication which are intrinsic to a Koch-type curve, I formulate and solve a wave equation that describes wave propagation along a fractal coastline. As opposed to the examples known from the literature I do not replace the fractal by the continuum in which it is embedded. This seems to be the first example of a truly intrinsic description of wave propagation along a fractal curve.
• Symmetries in the time-averaged dynamics of networks: reducing unnecessary complexity through minimal network models(1710.05251)

Dec. 22, 2018 nlin.CD, physics.soc-ph
Complex networks are the subject of fundamental interest from the scientific community at large. Several metrics have been introduced to characterize the structure of these networks, such as the degree distribution, degree correlation, path length, clustering coefficient, centrality measures etc. Another important feature is the presence of network symmetries. In particular, the effect of these symmetries has been studied in the context of network synchronization, where they have been used to predict the emergence and stability of cluster synchronous states. Here we provide theoretical, numerical, and experimental evidence that network symmetries play a role in a substantially broader class of dynamical models on networks, including epidemics, game theory, communication, and coupled excitable systems. Namely, we see that in all these models, nodes that are related by a symmetry relation show the same time-averaged dynamical properties. This discovery leads us to propose reduction techniques for exact, yet minimal, simulation of complex networks dynamics, which we show are effective in order to optimize the use of computational resources, such as computation time and memory.
• Stochastic resonance in a sinusoidal potential system: An analog simulation experiment(1608.06138)

Dec. 20, 2018 nlin.CD, cond-mat.stat-mech
Recently, stochastic resonance was obtained numerically in an underdamped periodic potential system driven by a periodic force and a Gaussian white noise. In that numerical work, the occurrence of stochastic resonance was explained in terms of the existence of two dynamical states having different amplitude and phase lag. At zero temperature these two initial condition dependent dynamical states are stable. However, at elevated temperatures, these two states make transitions from one to the other at a mean rate. In the present work, we setup an analog simulation experiment to show the existence of the two dynamical states in a sinusoidal potential system as well as to verify the occurrence of stochastic resonance in the same system. The experimental procedure includes setting the initial conditions for the experiment.
• Dynamic Analysis of Digital Chaotic Maps via State-Mapping Networks(1410.7694)

Dec. 18, 2018 nlin.CD, cs.CR
Chaotic dynamics is widely used to design pseudo-random number generators and for other applications such as secure communications and encryption. This paper aims to study the dynamics of discrete-time chaotic maps in the digital (i.e., finite-precision) domain. Differing from the traditional approaches treating a digital chaotic map as a black box with different explanations according to the test results of the output, the dynamical properties of such chaotic maps are first explored with a fixed-point arithmetic, using the Logistic map and the Tent map as two representative examples, from a new perspective with the corresponding state-mapping networks (SMNs). In an SMN, every possible value in the digital domain is considered as a node and the mapping relationship between any pair of nodes is a directed edge. The scale-free properties of the Logistic map's SMN are proved. The analytic results are further extended to the scenario of floating-point arithmetic and for other chaotic maps. Understanding the network structure of a chaotic map's SMN in digital computers can facilitate counteracting the undesirable degeneration of chaotic dynamics in finite-precision domains, helping also classify and improve the randomness of pseudo-random number sequences generated by iterating chaotic maps.
• Local invariants in non-ideal flows of neutral fluids and two-fluid plasmas(1709.03356)

Local structures, beyond the well-known frozen-in' to the barotropic flows of the generalized vorticities, of the two-fluid model of plasma flows are presented. More general non-barotropic situations are also considered. A modified Euler equation [T. Tao, Ann. PDE \textbf{2}, 9 (2016)] is also accordingly analyzed and remarked from the angle of view of two-fluid model, with emphasis on the local structures. And, the local constraints of high-order differential forms such as helicity, among others, find simple formulation for possible applications in modelling the dynamics. Thus, the Cauchy invariants equation [N. Besse and U. Frisch, J. Fluid Mech. \textbf{825}, 412 (2017)] may find practical application in non-ideal flows. Some formal examples are offered to outline the calculations, and particularly interestingly the two-dimensional-three-component (2D3C) or the 2D passive scalar problem presents that a locally invariant $\Theta = 2\theta \zeta$, with $\theta$ and $\zeta$ being respectively the scalar value of the vertical velocity' (or the passive scalar) and the vertical vorticity', may be used as if it were the spatial density of the globally invariant helicity, providing a Lagrangian prescription to control the latter in some situations of studying its physical effects in rapidly rotating flows (ubiquitous in atmosphere of astrophysical objects) with marked 2D3C vortical modes or in purely 2D passive scalars.
• Diffusion on hierarchical systems of weakly-coupled networks(1303.2650)

Dec. 13, 2018 nlin.CD, physics.soc-ph
We analyse diffusion dynamics on weakly-coupled networks (interconnected networks) by means of separation of time scales. Using an adiabatic approximation we reduced the system dynamics to a Markov chain with aggregated variables and derived a transport equation that is analogous to Fick's First Law and includes a driving force. Entropy production is a sum of microscopic entropy transport, which results from the particle's migration between networks of different topologies and macroscopic entropy production of the Markov chain. Equilibrium particles partition between different sub-networks depends only on internal sub-network parameters. Our framework, confirmed by numerical simulations, is also useful for considering diffusion in nested systems corresponding to hierarchical networks with several different time scales thus it can serve to uncover hidden hierarchy levels from observations of diffusion processes.
• Asymmetric cluster and chimera dynamics in globally coupled systems(1709.08663)

Nov. 23, 2018 nlin.CD
We investigate the emergence of chimera and cluster states possessing asymmetric dynamics in globally coupled systems, where the trajectories of oscillators belonging to different subpopulations exhibit different dynamical properties. In an asymmetric chimera state, the trajectory of an element in the synchronized subset is stationary or periodic, while that of an oscillator in the desynchronized subset is chaotic. In an asymmetric cluster state, the periods of the trajectories of elements belonging to different clusters are different. We consider a network of globally coupled chaotic maps as a simple model for the occurrence of such asymmetric states in spatiotemporal systems. We employ the analogy between a single map subject to a constant drive and the effective local dynamics in the globally coupled map system to elucidate the mechanisms for the emergence of asymmetric chimera and cluster states in the latter system. By obtaining the dynamical responses of the driven map, we establish a condition for the equivalence of the dynamics of the driven map and that of the system of globally coupled maps. This condition is applied to predict parameter values and subset partitions for the formation of asymmetric cluster and chimera states in the globally coupled system.
• Dynamics of non-holonomic systems with stochastic transport(1707.04741)

This paper formulates a variational approach for treating observational uncertainty and/or computational model errors as stochastic transport in dynamical systems governed by action principles under nonholonomic constraints. For this purpose, we derive, analyze and numerically study the example of an unbalanced spherical ball rolling under gravity along a stochastic path. Our approach uses the Hamilton-Pontryagin variational principle, constrained by a stochastic rolling condition, which we show is equivalent to the corresponding stochastic Lagrange-d'Alembert principle. In the example of the rolling ball, the stochasticity represents uncertainty in the observation and/or error in the computational simulation of the angular velocity of rolling. The influence of the stochasticity on the deterministically conserved quantities is investigated both analytically and numerically. Our approach applies to a wide variety of stochastic, nonholonomically constrained systems, because it preserves the mathematical properties inherited from the variational principle. Keywords: Nonholonomic constraints, Stochastic dynamics, Transport noise.
• Stochastic Parametrization of the Richardson Triple(1708.04183)

A Richardson triple is an ideal fluid flow map $g_{t/\ep,t,\ep t} = h_{t/\ep}k_t l_{\ep t}$ composed of three smooth maps with separated time scales: slow, intermediate and fast; corresponding to the big, little, and lesser whorls in Richardson's well-known metaphor for turbulence. Under homogenisation, as $\lim \ep\to0$, the composition $h_{t/\ep}k_t$ of the fast flow and the intermediate flow is known to be describable as a single stochastic flow $\dd g$. The interaction of the homogenised stochastic flow $\dd g$ with the slow flow of the big whorl is obtained by going into its non-inertial moving reference frame, via the composition of maps $(\dd g)l_{\ep t}$. This procedure parameterises the interactions of the three flow components of the Richardson triple as a single stochastic fluid flow in a moving reference frame. The Kelvin circulation theorem for the stochastic dynamics of the Richardson triple reveals the interactions among its three components. Namely, (i) the velocity in the circulation integrand acquires is kinematically swept by the large scales; and (ii) the velocity of the material circulation loop acquires additional stochastic Lie transport by the small scales. The stochastic dynamics of the composite homogenised flow is derived from a stochastic Hamilton's principle, and then recast into Lie-Poisson bracket form with a stochastic Hamiltonian. Several examples are given, including fluid flow with stochastically advected quantities, and rigid body motion under gravity, i.e., the stochastic heavy top in a rotating frame.
• Particle-relabeling symmetry, generalized vorticity, and normal-mode expansion of ideal incompressible fluids and plasmas in three-dimensional space(1601.05477)

Oct. 20, 2018 physics.plasm-ph, nlin.CD
The Lagrangian mechanical consideration of the dynamics of ideal incompressible hydrodynamic, magnetohydrodynamic, and Hall magnetohydrodynamic media, which are formulated as dynamical systems in appropriate Lie groups equipped with Riemannian metrics, leads to the notion of generalized vorticities, as well as generalized coordinates, velocities, and momenta. The action of each system is conserved against the integral path variation in the direction of the generalized vorticity, and this invariance is associated with the particle relabeling symmetry. The generalized vorticities are formulated by the operation of integro-differential operators on the generalized velocities. The eigenfunctions of the operators provide sets of orthogonal functions, and we obtain common mathematical expressions concerning these dynamical systems using the orthogonoal functions. In particular, we find that the product of the Riemannian metric, $g_{im}$, and the structure constants of the Lie group, $C^{m}_{jk}$, is given by the product of the eigenvalue of the operator, $\Lambda(i)$, and a certain totally antisymmetric tensor, $T_{ijk}$: $g_{i\alpha}C^{\alpha}_{jk}=\Lambda(i)T_{ijk}.$ Its physical implications, including the weak interaction conjecture of MHD turbulence, are also discussed.
• On Noether's theorem for the Euler-Poincar\'e equation on the diffeomorphism group with advected quantities(1206.2976)

We show how Noether conservation laws can be obtained from the particle relabelling symmetries in the Euler-Poincar\'e theory of ideal fluids with advected quantities. All calculations can be performed without Lagrangian variables, by using the Eulerian vector fields that generate the symmetries, and we identify the time-evolution equation that these vector fields satisfy. When advected quantities (such as advected scalars or densities) are present, there is an additional constraint that the vector fields must leave the advected quantities invariant. We show that if this constraint is satisfied initially then it will be satisfied for all times. We then show how to solve these constraint equations in various examples to obtain evolution equations from the conservation laws. We also discuss some fluid conservation laws in the Euler-Poincar\'e theory that do not arise from Noether symmetries, and explain the relationship between the conservation laws obtained here, and the Kelvin-Noether theorem given in Section 4 of Holm, Marsden and Ratiu, {\it Adv. in Math.}, 1998.
• Selective decay by Casimir dissipation in fluids(1206.2607)

The problem of parameterizing the interactions of larger scales and smaller scales in fluid flows is addressed by considering a property of two-dimensional incompressible turbulence. The property we consider is selective decay, in which a Casimir of the ideal formulation (enstrophy in 2D flows, helicity in 3D flows) decays in time, while the energy stays essentially constant. This paper introduces a mechanism that produces selective decay by enforcing Casimir dissipation in fluid dynamics. This mechanism turns out to be related in certain cases to the numerical method of anticipated vorticity discussed in \cite{SaBa1981,SaBa1985}. Several examples are given and a general theory of selective decay is developed that uses the Lie-Poisson structure of the ideal theory. A scale-selection operator allows the resulting modifications of the fluid motion equations to be interpreted in several examples as parameterizing the nonlinear, dynamical interactions between disparate scales. The type of modified fluid equation systems derived here may be useful in modelling turbulent geophysical flows where it is computationally prohibitive to rely on the slower, indirect effects of a realistic viscosity, such as in large-scale, coherent, oceanic flows interacting with much smaller eddies.
• Stochastic geometric models with non-stationary spatial correlations in Lagrangian fluid flows(1703.06774)

Oct. 18, 2018 nlin.CD
Inspired by spatiotemporal observations from satellites of the trajectories of objects drifting near the surface of the ocean in the National Oceanic and Atmospheric Administration's Global Drifter Program', this paper develops data-driven stochastic models of geophysical fluid dynamics (GFD) with non-stationary spatial correlations representing the dynamical behaviour of oceanic currents. Three models are considered. Model 1 from Holm [2015] is reviewed, in which the spatial correlations are time independent. Two new models, called Model 2 and Model 3, introduce two different symmetry breaking mechanisms by which the spatial correlations may be advected by the flow. These models are derived using reduction by symmetry of stochastic variational principles, leading to stochastic Hamiltonian systems, whose momentum maps, conservation laws and Lie-Poisson bracket structures are used in developing the new stochastic Hamiltonian models of GFD.
• Exploring quantum chaos with a single nuclear spin(1703.04852)

Most classical dynamical systems are chaotic. The trajectories of two identical systems prepared in infinitesimally different initial conditions diverge exponentially with time. Quantum systems, instead, exhibit quasi-periodicity due to their discrete spectrum. Nonetheless, the dynamics of quantum systems whose classical counterparts are chaotic are expected to show some features that resemble chaotic motion. Among the many controversial aspects of the quantum-classical boundary, the emergence of chaos remains among the least experimentally verified. Time-resolved observations of quantum chaotic dynamics are particularly rare, and as yet unachieved in a single particle, where the subtle interplay between chaos and quantum measurement could be explored at its deepest levels. We present here a realistic proposal to construct a chaotic driven top from the nuclear spin of a single donor atom in silicon, in the presence of a nuclear quadrupole interaction. This system is exquisitely measurable and controllable, and possesses extremely long intrinsic quantum coherence times, allowing for the observation of subtle dynamical behavior over extended periods. We show that signatures of chaos are expected to arise for experimentally realizable parameters of the system, allowing the study of the relation between quantum decoherence and classical chaos, and the observation of dynamical tunneling.
• New elements for a network (including brain) general theory during learning period(1301.2959)

Oct. 2, 2018 nlin.CD, cs.NE, nlin.AO
This study deals with the evolution of the so called 'intelligent' networks (insect society without leader, cells of an organism, brain,...) during their learning period. First we summarize briefly the Version 2 (published in French), whose the main characteristics are: 1) A network connected to its environment is considered as immersed into an information field created by this environment which so dictates to it the learning constraints. 2) The used formalism draws one's inspiration from the one of the Quantum field theory (Principle of stationary action, gauge fields, invariance by symmetry transformations,...). 3) We obtain Lagrange equations whose solutions describe the network evolution during the whole learning period. 4) Then, while proceeding with the same formalism inspiration, we suggest other study ways capable of evolving the knowledge in the considered scope. In a second part, after a reminder of the points to be improved, we exhibit the Version 5 which brings, we think, relevant improvements. Indeed: 5) We consider the weighted averages of the variables; this introduces probabilities. 6) We define two observables (L average of information flux and A activity of the network) which could be measured and so be compared with experimental results. 7) We find that L , weighted average of information flows, is an invariant. 8) Finally, we propose two expressions for the conactance, from which we deduce the corresponding Lagrange equations which have to be solved to know the evolution of the considered weighted averages. But, at the present stage, we think that we can progress only by carrying out experiments (see projects like Human brain project) and discovering invariants, symmetries which would allow us, like in Physics, to classify networks and above all to understand better the connections between them. Indeed, and that is what we propose among the future research ways, the underlying problem is to understand how, after their learning period, several networks can connect together to produce, in the brain case for instance, what we call mental states.