Non-KAM (Kolmogorov-Arnold-Moser) systems, when perturbed by weak time-dependent fields, offer a fast route to the classical chaos through an abrupt breaking of the invariant phase space tori. However, such behavior is not ubiquitous but rather contingent on whether the total system is in resonance. The resonances are usually determined by the ratios of characteristic frequencies associated with the system and the perturbation. Under the resonance condition, the classical dynamics are highly susceptible to variations in the system parameters. In this work, we employ out-of-time-order correlators (OTOCs) to study the dynamical sensitivity of a perturbed non-KAM system in the quantum limit as the parameter that characterizes the resonances and non-resonances is slowly varied. For this purpose, we consider a quantized kicked harmonic oscillator (KHO) model with the kick being the external time-dependent perturbation. Although the Lyapunov exponent of the KHO at resonances remains close to zero in the weak perturbative regime, making the system weakly chaotic in the conventional sense, the classical phase space undergoes significant structural changes. Motivated by this, we study the OTOCs when the system is in resonance and contrast the results with the non-resonant case. At resonances, we observe that the asymptotic dynamics of the OTOCs are sensitive to these structural changes, where they grow quadratically as opposed to linear or stagnant growth at non-resonances. On the other hand, our findings suggest that the short-time dynamics remain relatively more stable to the variations in the parameter. We will back our results by providing analytical expressions for the OTOCs for a few special cases. We will then extend our findings concerning the non-resonant cases to a broad class of KAM systems.

The predictability of weather and climate is strongly state-dependent: special and extremely relevant atmospheric states like blockings are associated with anomalous instability. Indeed, typically, the instability of a chaotic dynamical system can vary considerably across its attractor. Such an attractor is in general densely populated by unstable periodic orbits that can be used to approximate any forward trajectory through the so-called shadowing. Dynamical heterogeneity can lead to the presence of unstable periodic orbits with different number of unstable dimensions. This phenomenon - unstable dimensions variability - implies a serious breakdown of hyperbolicity and has considerable implications in terms of the structural stability of the system and of the possibility to describe accurately its behaviour through numerical models. As a step in the direction of better understanding the properties of high-dimensional chaotic systems, we provide here an extensive numerical study of the dynamical heterogeneity of the Lorenz '96 model in a parametric configuration leading to chaotic dynamics. We show that the detected variability in the number of unstable dimensions is associated with the presence of many finite-time Lyapunov exponents that fluctuate about zero also when very long averaging times are considered. The transition between regions of the attractor with different degrees of instability comes with a significant drop of the quality of the shadowing. By performing a coarse graining based on the shadowing unstable periodic orbits, we can characterize the slow fluctuations of the system between regions featuring, on the average, anomalously high and anomalously low instability. In turn, such regions are associated, respectively, with states of anomalously high and low energy, thus providing a clear link between the microscopic and thermodynamical properties of the system.

We investigate the dynamical evolution of globally connected Stuart-Landau oscillators coupled through conjugate or dis-similar variables on simplicial complexes. We report a first-order explosive phase transition from oscillatory state to death state, with 2-simplex (triadic) interactions, as opposed to the second-order transition with only 1-simplex (dyadic) interactions. Moreover, the system displays four distinct homogeneous steady states in the presence of triadic interactions, in contrast to the two homogeneous steady states observed with dyadic interactions. We calculate the backward transition point analytically, confirming the numerical results and providing the origin of the dynamical states in the transition region. The study will be useful in understanding complex systems, such as ecological and epidemiological, having higher-order interactions and coupling through conjugate variables.

Time-delayed optical feedback is known to trigger a wide variety of complex dynamical behavior in semiconductor lasers. Adding a second optical feedback loop is naturally expected to further increase the complexity of the system and its dynamics, but due to interference between the two feedback arms it was also quickly identified as a way to improve the laser stability. While these two aspects have already been investigated, the influence of the feedback phases, i.e. sub-wavelength changes in the mirror positions, on the laser behavior still remains to be thoroughly studied, despite indications that this parameter could have a significant impact. Here, we analyze the effect of the feedback phase on the laser stability in a dual-feedback configuration. We show an increased sensitivity of the laser system to feedback phase changes when two feedback loops are present, and clarify the interplay between the frequency shift induced by the feedback and the interferometric effect between the two feedback arms.

Chaotic dynamical systems (DS) are ubiquitous in nature and society. Often we are interested in reconstructing such systems from observed time series for prediction or mechanistic insight, where by reconstruction we mean learning geometrical and invariant temporal properties of the system in question (like attractors). However, training reconstruction algorithms like recurrent neural networks (RNNs) on such systems by gradient-descent based techniques faces severe challenges. This is mainly due to exploding gradients caused by the exponential divergence of trajectories in chaotic systems. Moreover, for (scientific) interpretability we wish to have as low dimensional reconstructions as possible, preferably in a model which is mathematically tractable. Here we report that a surprisingly simple modification of teacher forcing leads to provably strictly all-time bounded gradients in training on chaotic systems, and, when paired with a simple architectural rearrangement of a tractable RNN design, piecewise-linear RNNs (PLRNNs), allows for faithful reconstruction in spaces of at most the dimensionality of the observed system. We show on several DS that with these amendments we can reconstruct DS better than current SOTA algorithms, in much lower dimensions. Performance differences were particularly compelling on real world data with which most other methods severely struggled. This work thus led to a simple yet powerful DS reconstruction algorithm which is highly interpretable at the same time.

This work presents the building-blocks of an integrability-based representation for multi-point Fishnet Feynman integrals with any number of loops. Such representation relies on the quantum separation of variables (SoV) of a non-compact spin-chain with symmetry $SO(1,5)$ explained in the first paper of this series. The building-blocks of the SoV representation are overlaps of the wave-functions of the spin-chain excitations inserted along the edges of a triangular tile of Fishnet lattice. The zoology of overlaps is analyzed along with various worked out instances in order to achieve compact formulae for the generic triangular tile. The procedure of assembling the tiles into a Fishnet integral is presented exhaustively. The present analysis describes multi-point correlators with disk topology in the bi-scalar limit of planar $\gamma$-deformed $\mathcal{N}=4$ SYM theory, and it verifies some conjectural formulae for hexagonalization of Fishnets CFTs present in the literature. The findings of this work are suitable of generalization to a wider class of Feynman diagrams.