<center><b> Welcome To Nonlinear Sciences</b></center>
# <center>Adaptation And Self-Organizing Systems</center>
<hr>
### [Domain expansion and transient scaling regimes in population networks with in-domain cyclic selection](http://arxiv.org/abs/1811.07412v1) (1811.07412v1)
<i>P. P. Avelino, J. Menezes, B. F. de Oliveira, T. A. Pereira</b>
<h10>2018-11-18</h10>
> We investigate a six-species class of May-Leonard models which leads to the formation of networks with two different types of domains (three species each). In these models every species in each domain is in direct competition with all the species in the other domain (bidirectional selection, with a fixed probability ). On the other hand, the species within each domain compete cyclically (unidirectional selection, with a fixed probability ), leading to the formation of spiral waves inside the domains. We study the evolution of such networks using stochastic numerical simulations. We show that single species domains may form during the final stages of the collapse of three-species domains and we determine the critical initial radius beyond which single-species spherically symmetric domains are expected to expand as a result of their larger density of individuals compared to that of the outer three-species domain. We further show that a transient scaling regime, with a slower average growth rate of the characteristic length scale  with time , takes place before the transition to a standard  scaling law (usually associated with curvature driven dynamics), resulting in an extended period of coexistence.
### [Complex self-sustained oscillation patterns in modular excitable networks](http://arxiv.org/abs/1807.08849v2) (1807.08849v2)
<i>Jason Danison, Miguel Perez</b>
<h10>2018-07-23</h10>
> We study the relationship between the modularity of scale-free excitable networks and their ability to support self-sustained oscillation patterns. We find that the probability for a network of given degree-distribution exponent to be able to support self-sustained oscillations is strongly affected by its modularity. In addition, both high- and low-modularity networks are more likely to exhibit long-period oscillation patterns than those with intermediate modularity, but the degrees of complexity and correlation in these two cases are different. The long-period oscillations cannot be explained by a minimum-length Winfree loop, but instead arise from the interplay between two or more propagating waves. Finally, we introduce a new method that can be used to analyze the structure of the self-sustained oscillation sources at different levels of detail and show that the period of the oscillation pattern is statistically correlated with the fraction of modules that are part of the oscillation source.
### [Multivariate Spatiotemporal Hawkes Processes and Network Reconstruction](http://arxiv.org/abs/1811.06321v1) (1811.06321v1)
<i>Baichuan Yuan, Hao Li, Andrea L. Bertozzi, P. Jeffrey Brantingham, Mason A. Porter</b>
<h10>2018-11-15</h10>
> There is often latent network structure in spatial and temporal data and the tools of network analysis can yield fascinating insights into such data. In this paper, we develop a nonparametric method for network reconstruction from spatiotemporal data sets using multivariate Hawkes processes. In contrast to prior work on network reconstruction with point-process models, which has often focused on exclusively temporal information, our approach uses both temporal and spatial information and does not assume a specific parametric form of network dynamics. This leads to an effective way of recovering an underlying network. We illustrate our approach using both synthetic networks and networks constructed from real-world data sets (a location-based social media network, a narrative of crime events, and violent gang crimes). Our results demonstrate that, in comparison to using only temporal data, our spatiotemporal approach yields improved network reconstruction, providing a basis for meaningful subsequent analysis --- such as community structure and motif analysis --- of the reconstructed networks.
### [Interacting particles with Lévy strategies: limits of transport equations for swarm robotic systems](http://arxiv.org/abs/1807.10124v3) (1807.10124v3)
<i>Gissell Estrada-Rodriguez, Heiko Gimperlein</b>
<h10>2018-07-18</h10>
> L\'{e}vy robotic systems combine superdiffusive random movement with emergent collective behaviour from local communication and alignment in order to find rare targets or track objects. In this article we derive macroscopic fractional PDE descriptions from the movement strategies of the individual robots. Starting from a kinetic equation which describes the movement of robots based on alignment, collisions and occasional long distance runs according to a L\'{e}vy distribution, we obtain a system of evolution equations for the fractional diffusion for long times. We show that the system allows efficient parameter studies for a search problem, addressing basic questions like the optimal number of robots needed to cover an area in a certain time. For shorter times, in the hyperbolic limit of the kinetic equation, the PDE model is dominated by alignment, irrespective of the long range movement. This is in agreement with previous results in swarming of self-propelled particles. The article indicates the novel and quantitative modeling opportunities which swarm robotic systems provide for the study of both emergent collective behaviour and anomalous diffusion, on the respective time scales.
### [Assessing biological models using topological data analysis](http://arxiv.org/abs/1811.04827v1) (1811.04827v1)
<i>M. Ulmer, Lori Ziegelmeier, Chad M. Topaz</b>
<h10>2018-11-12</h10>
> We use topological data analysis as a tool to analyze the fit of mathematical models to experimental data. This study is built on data obtained from motion tracking groups of aphids in [Nilsen et al., PLOS One, 2013] and two random walk models that were proposed to describe the data. One model incorporates social interactions between the insects, and the second model is a control model that excludes these interactions. We compare data from each model to data from experiment by performing statistical tests based on three different sets of measures. First, we use time series of order parameters commonly used in collective motion studies. These order parameters measure the overall polarization and angular momentum of the group, and do not rely on a priori knowledge of the models that produced the data. Second, we use order parameter time series that do rely on a priori knowledge, namely average distance to nearest neighbor and percentage of aphids moving. Third, we use computational persistent homology to calculate topological signatures of the data. Analysis of the a priori order parameters indicates that the interactive model better describes the experimental data than the control model does. The topological approach performs as well as these a priori order parameters and better than the other order parameters, suggesting the utility of the topological approach in the absence of specific knowledge of mechanisms underlying the data.
# <center>Chaotic Dynamics</center>
<hr>
### [Spectral Test of the MIXMAX Random Number Generators](http://arxiv.org/abs/1806.05243v2) (1806.05243v2)
<i>Narek Martirosyan, Konstantin Savvidy, George Savvidy</b>
<h10>2018-06-13</h10>
> An important statistical test on the pseudo-random number generators is called the spectral test. The test is aimed at answering the question of distribution of the generated pseudo-random vectors in dimensions  that are larger than the genuine dimension of a generator . In particular, the default MIXMAX generators have various dimensions:  and higher. Therefore the spectral test is important to perform in dimensions  for  generator,  for  and  for  generator. These tests have been performed by L'Ecuyer and collaborators. When  the vectors of the generated numbers fall into the parallel hyperplanes and the distances between them can be larger than the genuine "resolution" of the MIXMAX generators, which is . The aim of this article is to further study the spectral properties of the MIXMAX generators, to investigate the dependence of the spectral properties of the MIXMAX generators as a function of their internal parameters and in particular their dependence on the parameter . We found that the best spectral properties are realized when  is between  and , a range which is inclusive of the value of the  generator. We also provide the alternative parameters for the generators,  and  with  in this optimised range.
### [Symmetries in the time-averaged dynamics of networks: reducing unnecessary complexity through minimal network models](http://arxiv.org/abs/1710.05251v2) (1710.05251v2)
<i>Abu Bakar Siddique, Louis M. Pecora, Francesco Sorrentino</b>
<h10>2017-10-15</h10>
> 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.
### [Radio-wave communication with chaos](http://arxiv.org/abs/1811.07334v1) (1811.07334v1)
<i>Hai-Peng Ren, Wu-Yun Zheng, Celso Grebogi</b>
<h10>2018-11-18</h10>
> The invariance of the Lyapunov exponent of a chaotic signal as it propagates along a wireless transmission channel provides a theoretical base for the application of chaos in wireless communication. In additive Gaussian channel, the chaotic signal is proved to be the optimal coherent communication waveform in the sense of using the very simple matched filter to obtain the maximum signal-to-noise ratio. The properties of chaos can be used to reduce simply and effectively the Inter-Symbol Interference (ISI) and to achieve low bit error rate in the wireless communication system. However, chaotic signals need very wide bandwidth to be transmitted in the practical channel, which is difficult for the practical transducer or antenna to convert such a broad band signal. To solve this problem, in this work, the chaotic signal is applied to a radio-wave communication system, and the corresponding coding and decoding algorithms are proposed. A hybrid chaotic system is used as the pulse-shaping filter to obtain the baseband signal, and the corresponding matched filter is used at the receiver, instead of the conventional low-pass filter, to maximize the signal-to-noise ratio. At the same time, the symbol judgment threshold determined by the chaos property is used to reduce the Inter-Symbol Interference (ISI) effect. Simulations and virtual channel experiments show that the radio-wave communication system using chaos obtains lower bit error rate in the multi-path transmission channel compared with the traditional radio-wave communication system using Binary Phase Shift Keying (BPSK) modulation technology and channel equalization.
### [Dynamics Motivated by Sierpinski Fractals](http://arxiv.org/abs/1811.07122v1) (1811.07122v1)
<i>Marat Akhmet, Mehmet Onur Fen, Ejaily Milad Alejaily</b>
<h10>2018-11-17</h10>
> Fatou-Julia iteration (FJI) is an effective instrument to construct fractals. Famous Julia and Mandelbrot sets are strong confirmations of this. In the present study, we use the paradigm of FJI to construct and map Sierpinski fractals. The fractals can be mapped by developing a mapping iteration on the basis of FJI. Because of the close link between mappings, differential equations and dynamical systems, one can introduce dynamics for a fractal through differential equations such that it becomes points of the solution trajectory. Thus, in this paper, we consider two types of dynamics motivated by the Sierpinski fractals. The first one is the dynamics of FJI itself, and the second one is the dynamics of a fractal mapping iteration which can be performed through differential equations. The characterization of fractals as trajectory points of the dynamics can help to enhance and widen the scope of their applications in physics and engineering.
### [Resonator neuron and triggering multipulse excitability in laser with injected signal](http://arxiv.org/abs/1811.06820v1) (1811.06820v1)
<i>Axel Dolcemascolo, Bruno Garbin, Bertand Peyce, Romain Veltz, Stephane Barland</b>
<h10>2018-11-16</h10>
> Semiconductor lasers with coherent forcing are expected to behave similarly to simple neuron models in response to external perturbations, as long as the physics describing them can be approximated by that of an overdamped pendulum with fluid torque. Beyond the validity range of this approximation, more complex features can be expected. We perform experiments and numerical simulations which show that the system can display resonator and integrator features depending on parameters and that multiple pulses can be emitted in response to larger perturbations.
# <center>Cellular Automata And Lattice Gases</center>
<hr>
### [Theoretical investigation, simulation and empirical analysis of the growth pattern of traffic oscillations in the Euler coordinates](http://arxiv.org/abs/1806.04287v2) (1806.04287v2)
<i>Junfang Tian, Martin Treiber, Rui Jiang, Xiaopeng Li, Bin Jia, Ziyou Gao</b>
<h10>2018-06-12</h10>
> The formation and development of oscillations is an important traffic flow phenomenon. Recent studies found that along a vehicle platoon described in the Lagrangian specification, traffic oscillations grow in a concave way. Since stationary bottlenecks are more intuitively described in the Eulerian framework, this paper investigates whether the concave growth pattern of traffic oscillations in the Lagrangian coordinates can be transferred to the Euler coordinates (i.e. the concave increase in standard deviation is no longer measured across the vehicle indices but as a function of the road location). To this end, we theoretically unify these two ways of measuring oscillations by revealing their mapping relationship. We show that the growth pattern measured in the Lagrangian coordinates can be transferred to the Euler coordinates. We believe this finding is nontrivial since the scenarios are significantly different: while in vehicle platoons (Lagrangian view, non-penetrable moving bottleneck), the speed variance for a given vehicle is ideally constant, the drivers in the Eulerian setting (penetrable stationary bottleneck triggering the waves) experience all amplitudes, first the big ones and then the small ones. To test this proposition, we performed simulation using two different kinds of car-following models. Simulation results validate the theoretical analysis. Finally, we performed empirical analysis using the NGSIM data, which also validates the theoretical analysis.
### [One-sided scaling limit of multicolor box-ball system](http://arxiv.org/abs/1808.08074v3) (1808.08074v3)
<i>Atsuo Kuniba, Hanbaek Lyu</b>
<h10>2018-08-24</h10>
> A random box-ball system starts with occupying each of the first  boxes independently with a ball of random color from , where balls of color 0 are considered as empty boxes. The time evolution is defined by a successive application of the combinatorial , and possesses a -tuple of Young diagrams as the complete set of conserved quantities. We show that the induced carrier process is an irreducible Markov chain and that the associated energy function scaled by  converges as  at an exponential rate. This implies that each of the invariant Young diagrams converges to some limiting shape as  if we scale its rows by . By establishing a large deviations principle, we show that this limiting shape is unchanged under conditioning the initial measure on the highest states. Using these results, we give an explicit formula for the limiting shape of the invariant Young diagrams by Thermodynamic Bethe Ansatz.
### [Evacuation simulation considering action of the guard in an artificial attack](http://arxiv.org/abs/1810.05918v1) (1810.05918v1)
<i>Chang-kun Chen, Yun-he Tong</b>
<h10>2018-10-13</h10>
> To investigate the evacuation behaviors of pedestrians considering action of the guard and develop an effective evacuation strategy in the artificial attack, an extended floor field model was proposed. In this model, the assault on pedestrians, the death of pedestrians and the capture of the guard were involved simultaneously. An alternative evacuation strategy which can largely reduce the number of death was developed and effects of several key parameters such as the deterrence radius and capture distance on evacuation dynamics were studied. Results show that congestion near the exit has dual effects. More specially, the guard could catch all attackers in a short time because the attackers would have more concentrated distribution, but more casualties would happen because pedestrians are hard to escape the attack due to congestion. In contrast, when pedestrians have more preference of approaching the guard, although the guard would take more time to capture the attackers result from the dispersion of attackers, the death toll would decrease. One of the reason is the dispersal of the crowd and the decrease in congestion would be beneficial for escape. Another is the attackers would be killed before launching the attack to the people those are around the guard, in other words, the guard would protect a large number of pedestrians from being killed. Moreover, increasing capture distance of the guard can effectively reduce the casualties and the catch time. As the deterrence radius reflecting the tendency of escaping from the guard for attackers rises, it would become more difficult for the guard to catch the attackers and more casualties are caused. However, when the deterrence radius reaches a certain level, the number of deaths would be reduced because the attackers would prefer to stay as far away as possible from instead of the position where they could attack more people.
### [Von Neumann regularity, split epicness and elementary cellular automata](http://arxiv.org/abs/1804.03913v3) (1804.03913v3)
<i>Ville Salo</b>
<h10>2018-04-11</h10>
> We show that a cellular automaton on a mixing subshift of finite type is a Von Neumann regular element in the semigroup of cellular automata if and only if it is split epic onto its image in the category of sofic shifts and block maps. It follows from [S.-T\"orm\"a, 2015] that Von Neumann regularity is decidable condition, and we decide it for all elementary CA.
### [Cellular automata approach to synchronized traffic flow modelling](http://arxiv.org/abs/1805.05555v2) (1805.05555v2)
<i>Junfang Tian, Chenqiang Zhu, Rui Jiang</b>
<h10>2018-05-15</h10>
> Cellular automaton (CA) approach is an important theoretical framework for studying complex system behavior and has been widely applied in various research field. CA traffic flow models have the advantage of flexible evolution rules and high computation efficiency. Therefore, CA develops very quickly and has been widely applied in transportation field. In recent two decades, traffic flow study quickly developed, among which "synchronized flow" is perhaps one of the most important concepts and findings. Many new CA models have been proposed in this direction. This paper makes a review of development of CA models, concerning their ability to reproduce synchronized flow as well as traffic breakdown from free flow to synchronized flow. Finally, future directions have been discussed.
# <center>Pattern Formation And Solitons</center>
<hr>
### [Kink-kink and kink-antikink interactions with long-range tails](http://arxiv.org/abs/1811.07872v1) (1811.07872v1)
<i>Ivan C. Christov, Robert J. Decker, A. Demirkaya, Vakhid A. Gani, P. G. Kevrekidis, Avinash Khare, Avadh Saxena</b>
<h10>2018-11-19</h10>
> In this Letter, we address the general long-range interaction between kinks and antikinks, as well as kinks and kinks, in  field theories for . The kink-antikink interaction is generically attractive, while the kink-kink interaction is generically repulsive. We find that the force of interaction decays with the th power of separation, and we identify the general prefactor for {\it arbitrary} . Importantly, we test the resulting mathematical prediction with detailed numerical simulations of the dynamic field equation, and obtain {excellent} agreement between theory and numerics for the cases of  ( model),  ( model) and  ( model).
### [Stable and Unstable Vortex Knots in Excitable Media](http://arxiv.org/abs/1809.04567v2) (1809.04567v2)
<i>Jack Binysh, Carl A. Whitfield, Gareth P. Alexander</b>
<h10>2018-09-12</h10>
> We study the dynamics of knotted vortices in a bulk excitable medium using the FitzHugh-Nagumo model. From a systematic survey of all knots of at most eight crossings we establish that the generic behaviour is of unsteady, irregular dynamics, with prolonged periods of expansion of parts of the vortex. The mechanism for the length expansion is a long-range `wave slapping' interaction, analogous to that responsible for the annihilation of small vortex rings by larger ones. We also show that there are stable vortex geometries for certain knots; in addition to the unknot, trefoil and figure eight knots reported previously, we have found stable examples of the Whitehead link and  knot. We give a thorough characterisation of their geometry and steady state motion. For the unknot, trefoil and figure eight knots we greatly expand previous evidence that FitzHugh-Nagumo dynamics untangles initially complex geometries while preserving topology.
### [Multi-dark soliton solutions for the multi-component coupled Maccari system](http://arxiv.org/abs/1811.07504v1) (1811.07504v1)
<i>Tao Xu, Yong Chen</b>
<h10>2018-11-19</h10>
> Based on the KP hierarchy reduction method, the general multi-dark soliton solutions in Gram type determinant forms for the multi-component coupled Maccari system are constructed. Especially, the three-component coupled Maccari system comprised of two-component short waves and single-component long waves are discussed in detail. Besides, the dynamics of one and two dark-dark solitons are analyzed. It is shown that the collisions of two dark-dark solitons are elastic by asymptotic analysis. Additionally, the two dark-dark solitons bound states are studied through two different cases (stationary and moving cases). The bound states can exist up to arbitrary order in the stationary case, however, only two-soliton bound state exists in the moving case. Besides, the oblique stationary bound state can be generated for all possible combinations of nonlinearity coefficients consisted of positive, negative and mixed cases. Nevertheless, the parallel stationary and the moving bound states are only possible when nonlinearity coefficients take opposite signs.
### [Domain expansion and transient scaling regimes in population networks with in-domain cyclic selection](http://arxiv.org/abs/1811.07412v1) (1811.07412v1)
<i>P. P. Avelino, J. Menezes, B. F. de Oliveira, T. A. Pereira</b>
<h10>2018-11-18</h10>
> We investigate a six-species class of May-Leonard models which leads to the formation of networks with two different types of domains (three species each). In these models every species in each domain is in direct competition with all the species in the other domain (bidirectional selection, with a fixed probability ). On the other hand, the species within each domain compete cyclically (unidirectional selection, with a fixed probability ), leading to the formation of spiral waves inside the domains. We study the evolution of such networks using stochastic numerical simulations. We show that single species domains may form during the final stages of the collapse of three-species domains and we determine the critical initial radius beyond which single-species spherically symmetric domains are expected to expand as a result of their larger density of individuals compared to that of the outer three-species domain. We further show that a transient scaling regime, with a slower average growth rate of the characteristic length scale  with time , takes place before the transition to a standard  scaling law (usually associated with curvature driven dynamics), resulting in an extended period of coexistence.
### [Dynamical transitions between equilibria in a dissipative Klein-Gordon lattice](http://arxiv.org/abs/1809.07995v2) (1809.07995v2)
<i>D. J. Frantzeskakis, N. I. Karachalios, P. G. Kevrekidis, V. Koukouloyannis, K. Vetas</b>
<h10>2018-09-21</h10>
> We consider the energy landscape of a dissipative Klein-Gordon lattice with a  on-site potential. Our analysis is based on suitable energy arguments, combined with a discrete version of the \L{}ojasiewicz inequality, in order to justify the convergence to a single, nontrivial equilibrium for all initial configurations of the lattice. Then, global bifurcation theory is explored, to illustrate that in the discrete regime all linear states lead to nonlinear generalizations of equilibrium states. Direct numerical simulations reveal the rich structure of the equilibrium set, consisting of non-trivial topological (kink-shaped) interpolations between the adjacent minima of the on-site potential, and the wealth of dynamical convergence possibilities. These dynamical evolution results also provide insight on the potential stability of the equilibrium branches, and glimpses of the emerging global bifurcation structure, elucidating the role of the interplay between discreteness, nonlinearity and dissipation.
# <center>Exactly Solvable And Integrable Systems</center>
<hr>
### [On quantum separation of variables](http://arxiv.org/abs/1807.11572v2) (1807.11572v2)
<i>J. M. Maillet, G. Niccoli</b>
<h10>2018-07-30</h10>
> We present a new approach to construct the separate variables basis leading to the full characterization of the transfer matrix spectrum of quantum integrable lattice models. The basis is generated by the repeated action of the transfer matrix itself on a generically chosen state of the Hilbert space. The fusion relations for the transfer matrix, stemming from the Yang-Baxter algebra properties, provide the necessary closure relations to define the action of the transfer matrix on such a basis in terms of elementary local shifts, leading to a separate transfer matrix spectral problem. Hence our scheme extends to the quantum case a key feature of Liouville-Arnold classical integrability framework where the complete set of conserved charges defines both the level manifold and the flows on it leading to the construction of action-angle variables. We work in the framework of the quantum inverse scattering method. As a first example of our approach, we give the construction of such a basis for models associated to Y(gln) and argue how it extends to their trigonometric and elliptic versions. Then we show how our general scheme applies concretely to fundamental models associated to the Y(gl2) and Y(gl3) R-matrices leading to the full characterization of their spectrum. For Y(gl2) and its trigonometric deformation a particular case of our method reproduces Sklyanin's construction of separate variables. For Y(gl3) it gives new results, in particular through the proper identification of the shifts acting on the separate basis. We stress that our method also leads to the full characterization of the spectrum of other known quantum integrable lattice models, including in particular trigonometric and elliptic spin chains, open chains with general integrable boundaries, and further higher rank cases that we will describe in forthcoming publications.
### [Complete spectrum of quantum integrable lattice models associated to Y(gl(n)) by separation of variables](http://arxiv.org/abs/1810.11885v2) (1810.11885v2)
<i>J. M. Maillet, G. Niccoli</b>
<h10>2018-10-28</h10>
> We apply our new approach of quantum Separation of Variables (SoV) to the complete characterization of the transfer matrix spectrum of quantum integrable lattice models associated to gl(n)-invariant R-matrices in the fundamental representations. We consider lattices with N sites and quasi-periodic boundary conditions associated to an arbitrary twist K having simple spectrum (but not necessarily diagonalizable). In our approach the SoV basis is constructed in an universal manner starting from the direct use of the conserved charges of the models, i.e., from the commuting family of transfer matrices. Using the integrable structure of the models, incarnated in the hierarchy of transfer matrices fusion relations, we prove that our SoV basis indeed separates the spectrum of the corresponding transfer matrices. Moreover, the combined use of the fusion rules, of the known analytic properties of the transfer matrices and of the SoV basis allows us to obtain the complete characterization of the transfer matrix spectrum and to prove its simplicity. Any transfer matrix eigenvalue is completely characterized as a solution of a so-called quantum spectral curve equation that we obtain as a difference functional equation of order n. Namely, any eigenvalue satisfies this equation and any solution of this equation having prescribed properties leads to an eigenvalue. We construct the associated eigenvector, unique up to normalization, by computing its decomposition on the SoV basis that is of a factorized form written in terms of the powers of the corresponding eigenvalues. If the twist matrix K is diagonalizable with simple spectrum then the transfer matrix is also diagonalizable with simple spectrum. In that case, we give a construction of the Baxter Q-operator satisfying a T-Q equation of order n, the quantum spectral curve equation, involving the hierarchy of the fused transfer matrices.
### [Integrable coupled Linard-type systems with balanced loss and gain](http://arxiv.org/abs/1804.02366v3) (1804.02366v3)
<i>Debdeep Sinha, Pijush K. Ghosh</b>
<h10>2018-04-06</h10>
> A Hamiltonian formulation of generic many-particle systems with space-dependent balanced loss and gain coefficients is presented. It is shown that the balancing of loss and gain necessarily occurs in a pair-wise fashion. Further, using a suitable choice of co-ordinates, the Hamiltonian can always be reformulated as a many-particle system in the background of a pseudo-Euclidean metric and subjected to an analogous inhomogeneous magnetic field with a functional form that is identical with space-dependent loss/gain co-efficient.The resulting equations of motion from the Hamiltonian are a system of coupled Linard-type differential equations. Partially integrable systems are obtained for two distinct cases, namely, systems with (i) translational symmetry or (ii) rotational invariance in a pseudo-Euclidean space. A total number of  integrals of motion are constructed for a system of  particles, which are in involution, implying that two-particle systems are completely integrable. A few exact solutions for both the cases are presented for specific choices of the potential and space-dependent gain/loss co-efficients, which include periodic stable solutions. Quantization of the system is discussed with the construction of the integrals of motion for specific choices of the potential and gain-loss coefficients. A few quasi-exactly solvable models admitting bound states in appropriate Stoke wedges are presented.
### [Hypergeometric First Integrals of the Duffing and van der Pol Oscillators](http://arxiv.org/abs/1706.02506v3) (1706.02506v3)
<i>Tomasz Stachowiak</b>
<h10>2017-06-08</h10>
> The autonomous Duffing oscillator, and its van der Pol modification, are known to admit time-dependent first integrals for specific values of parameters. This corresponds to the existence of Darboux polynomials, and in fact more can be shown: that there exist Liouvillian first integrals which do not depend on time. They can be expressed in terms of the Gauss and Kummer hypergeometric functions, and are neither analytic, algebraic nor meromorphic. A criterion for this to happen in a general dynamical system is formulated as well.
### [On the Integrable Structure of Super Yang-Mills Scattering Amplitudes](http://arxiv.org/abs/1811.06324v1) (1811.06324v1)
<i>Nils Kanning</b>
<h10>2018-11-15</h10>
> The maximally supersymmetric Yang-Mills theory in four-dimensional Minkowski space is an exceptional model of mathematical physics. Even more so in the planar limit, where the theory is believed to be integrable. In particular, the tree-level scattering amplitudes were shown to be invariant under the Yangian of the superconformal algebra psu(2,2|4). This infinite-dimensional symmetry is a hallmark of integrability. In this dissertation we explore connections between these amplitudes and integrable models. Our aim is to lay foundations for an efficient integrability-based computation of amplitudes. To this end, we characterize Yangian invariants within the quantum inverse scattering method, which is an extensive toolbox for integrable spin chains. Making use of this setup, we develop methods for the construction of Yangian invariants. We show that the algebraic Bethe ansatz can be specialized to yield Yangian invariants for u(2). Our approach also allows to interpret these Yangian invariants as partition functions of vertex models. What is more, we establish a unitary Gra{\ss}mannian matrix model for the construction of u(p,q|m) Yangian invariants with oscillator representations. In a special case our formula reduces to the Brezin-Gross-Witten model. We apply an integral transformation due to Bargmann to our unitary Gra{\ss}mannian matrix model, which turns the oscillators into spinor helicity-like variables. Thereby we are led to a refined version of the Gra{\ss}mannian integral formula for certain amplitudes. The most decisive differences are that we work in Minkowski signature and that the integration contour is fixed to be a unitary group manifold. We compare Yangian invariants defined by our integral to amplitudes and recently introduced deformations thereof.
<br><hr> <center>Thank you for reading!<br> https://cdn.steemitimages.com/DQmbn3ovuKLM17k6aemZMrJj6iqKkYzXCYz5Qh1Fg7vPmRx/image.png <br> Don't forget to Follow and Resteem. @complexcity <br>Keeping everyone inform.</center>
hiveblocks