
A theoretical study of delayed feedback in spintorque nanooscillators is
presented. A macrospin geometry is considered, where selfsustained
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 selfinjected 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.

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.

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 2dimensional
systems, investigates their properties, shows their application to several
example systems, and presents their extension to higher dimensions.

We derive Edgeworth expansions that describe corrections to the Gaussian
limiting behaviour of slowfast systems. The Edgeworth expansion is achieved
using a semigroup formalism for the transfer operator, where a DuhamelDyson
series is used to asymptotically determine the corrections at any desired order
of the time scale parameter $\varepsilon$. The corrections involve integrals
over higherorder autocorrelation 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.

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 twomode Hamiltonian is developed for a
shallow arch/buckled beam. The resulting form of the potential energytwo
stable wells connected by rank1 saddle pointsshows 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.

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}, 438441.]. 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}, 114130] are thoroughly truncated,
leaving the system with positive definite helicity and energy, negative
temperature and the corresponding largescale concentration of vortical modes
are not allowed, unlike the incompressible case.

The "backward simulation" of a stochastic process is defined as the
stochastic dynamics that trace a timereversed 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.

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 crosscorrelation is low. Therefore, sequences in the Weyl sequence
class are expected to minimize the intersymbol interference. We evaluate the
upper bound of crosscorrelation and odd crosscorrelation of spreading
sequences in the Weyl sequence class and construct the optimization problem:
minimize the upper bound of the absolute values of crosscorrelation and odd
crosscorrelation. 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.

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 timeseries. 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 StuartLandau equation
to nonsinusoidal oscillators obtained as a consequence of our theorem.

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 KPZtype equation for the anomalous
scaling of the longitudinal velocity structure functions and relate the
intermittency parameter to the boundary Euler (Atype) 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
NambuGoldstone gapless mode that arises from a spontaneous breaking due to the
energy flux of the separate scale and time symmetries of the inviscid
NavierStokes 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.

Beginning with addition and multiplication which are intrinsic to a Kochtype
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.

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 timeaveraged 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.

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.

Chaotic dynamics is widely used to design pseudorandom number generators and
for other applications such as secure communications and encryption. This paper
aims to study the dynamics of discretetime chaotic maps in the digital (i.e.,
finiteprecision) 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 fixedpoint arithmetic, using the Logistic map and the
Tent map as two representative examples, from a new perspective with the
corresponding statemapping 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 scalefree properties of the Logistic
map's SMN are proved. The analytic results are further extended to the scenario
of floatingpoint 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
finiteprecision domains, helping also classify and improve the randomness of
pseudorandom number sequences generated by iterating chaotic maps.

Local structures, beyond the wellknown `frozenin' to the barotropic flows
of the generalized vorticities, of the twofluid model of plasma flows are
presented. More general nonbarotropic 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 twofluid model,
with emphasis on the local structures. And, the local constraints of highorder
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 nonideal flows. Some formal examples are offered
to outline the calculations, and particularly interestingly the
twodimensionalthreecomponent (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.

We analyse diffusion dynamics on weaklycoupled 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 subnetworks depends only on
internal subnetwork 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.

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.

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 HamiltonPontryagin variational principle, constrained by a
stochastic rolling condition, which we show is equivalent to the corresponding
stochastic Lagranged'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.

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 wellknown 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
noninertial 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 LiePoisson 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.

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 integrodifferential 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.

We show how Noether conservation laws can be obtained from the particle
relabelling symmetries in the EulerPoincar\'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 timeevolution 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
EulerPoincar\'e theory that do not arise from Noether symmetries, and explain
the relationship between the conservation laws obtained here, and the
KelvinNoether theorem given in Section 4 of Holm, Marsden and Ratiu, {\it Adv.
in Math.}, 1998.

The problem of parameterizing the interactions of larger scales and smaller
scales in fluid flows is addressed by considering a property of twodimensional
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 LiePoisson structure of the ideal
theory. A scaleselection 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 largescale,
coherent, oceanic flows interacting with much smaller eddies.

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
datadriven stochastic models of geophysical fluid dynamics (GFD) with
nonstationary 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 LiePoisson bracket structures are used in developing the
new stochastic Hamiltonian models of GFD.

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
quasiperiodicity 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 quantumclassical boundary, the emergence of chaos remains among
the least experimentally verified. Timeresolved 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.

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.