In the work of Colliander et al. (2020) a minimal lattice model was constructed describing the transfer of energy to high frequencies in the defocusing nonlinear Schrödinger equation. In the present ...work, we present a systematic study of the coherent structures, both standing and traveling, that arise in the context of this model. We find that the nonlinearly dispersive nature of the model is responsible for standing waves in the form of discrete compactons. On the other hand, analysis of the dynamical features of the simplest nontrivial variant of the model, namely the dimer case, yields both solutions where the intensity is trapped in a single site and solutions where the intensity moves between the two sites, which suggests the possibility of moving excitations in larger lattices. Such excitations are also suggested by the dynamical evolution associated with modulational instability. Our numerical computations confirm this expectation, and we systematically construct such traveling states as exact solutions in lattices of varying size, as well as explore their stability. A remarkable feature of these traveling lattice waves is that they are of “antidark” type, i.e., they are mounted on top of a non-vanishing background. These studies shed light on the existence, stability and dynamics of such standing and traveling states in 1+1 dimensions, and pave the way for exploration of corresponding configurations in higher dimensions.
•A minimal lattice model that describes the transfer of energy to high frequencies in the defocusing nonlinear Schrödinger equation supports both standing and traveling coherent structures.•Standing coherent structures take the form of compactons, which are solitary waves with compact support.•The dimer system can be solved exactly, and yields both solutions where the intensity is trapped in a single site and solutions where the intensity moves between the two lattice sites.•Lattice traveling waves, which take the form of a localized excitation superimposed on a nonzero background, can be constructed numerically for larger lattices.
We consider the Hamilton-Jacobi equationH(x,Du)+λ(x)u=c,x∈M, where M is a connected, closed and smooth Riemannian manifold. The functions H(x,p) and λ(x) are continuous. H(x,p) is convex, coercive ...with respect to p, and λ(x) changes the signs. The first breakthrough to this model was achieved by Jin-Yan-Zhao 11 under the Tonelli conditions. In this paper, we consider more detailed structure of the viscosity solution set and large time behavior of the viscosity solution on the Cauchy problem. To the best of our knowledge, it is the first detailed description of the large time behavior of the HJ equations with non-monotone dependence on the unknown function.
The first complete proof of Arnold diffusion—one of the most important problems in dynamical systems and mathematical physics Arnold diffusion, which concerns the appearance of chaos in classical ...mechanics, is one of the most important problems in the fields of dynamical systems and mathematical physics. Since it was discovered by Vladimir Arnold in 1963, it has attracted the efforts of some of the most prominent researchers in mathematics. The question is whether a typical perturbation of a particular system will result in chaotic or unstable dynamical phenomena. In this groundbreaking book, Vadim Kaloshin and Ke Zhang provide the first complete proof of Arnold diffusion, demonstrating that that there is topological instability for typical perturbations of five-dimensional integrable systems (two and a half degrees of freedom).This proof realizes a plan John Mather announced in 2003 but was unable to complete before his death. Kaloshin and Zhang follow Mather's strategy but emphasize a more Hamiltonian approach, tying together normal forms theory, hyperbolic theory, Mather theory, and weak KAM theory. Offering a complete, clean, and modern explanation of the steps involved in the proof, and a clear account of background material, this book is designed to be accessible to students as well as researchers. The result is a critical contribution to mathematical physics and dynamical systems, especially Hamiltonian systems.
Spin models are the prime example of simplified many-body Hamiltonians used to model complex, strongly correlated real-world materials. However, despite the simplified character of such models, their ...dynamics often cannot be simulated exactly on classical computers when the number of particles exceeds a few tens. For this reason, quantum simulation of spin Hamiltonians using the tools of atomic and molecular physics has become a very active field over the past years, using ultracold atoms or molecules in optical lattices, or trapped ions. All of these approaches have their own strengths and limitations. Here we report an alternative platform for the study of spin systems, using individual atoms trapped in tunable two-dimensional arrays of optical microtraps with arbitrary geometries, where filling fractions range from 60 to 100 per cent. When excited to high-energy Rydberg D states, the atoms undergo strong interactions whose anisotropic character opens the way to simulating exotic matter. We illustrate the versatility of our system by studying the dynamics of a quantum Ising-like spin-1/2 system in a transverse field with up to 30 spins, for a variety of geometries in one and two dimensions, and for a wide range of interaction strengths. For geometries where the anisotropy is expected to have small effects on the dynamics, we find excellent agreement with ab initio simulations of the spin-1/2 system, while for strongly anisotropic situations the multilevel structure of the D states has a measurable influence. Our findings establish arrays of single Rydberg atoms as a versatile platform for the study of quantum magnetism.
Abstract Weak form equation learning and surrogate modeling has proven to be computationally efficient and robust to measurement noise in a wide range of applications including ODE, PDE, and SDE ...discovery, as well as in coarse-graining applications, such as homogenization and mean-field descriptions of interacting particle systems. In this work we extend this coarse-graining capability to the setting of Hamiltonian dynamics which possess approximate symmetries associated with timescale separation. A smooth $$\varepsilon$$ ε -dependent Hamiltonian vector field $$X_\varepsilon$$ X ε possesses an approximate symmetry if the limiting vector field $$X_0=\lim _{\varepsilon \rightarrow 0}X_\varepsilon$$ X 0 = lim ε → 0 X ε possesses an exact symmetry. Such approximate symmetries often lead to the existence of a Hamiltonian system of reduced dimension that may be used to efficiently capture the dynamics of the symmetry-invariant dependent variables. Deriving such reduced systems, or approximating them numerically, is an ongoing challenge. We demonstrate that WSINDy can successfully identify this reduced Hamiltonian system in the presence of large perturbations imparted in the $$\varepsilon >0$$ ε > 0 regime, while remaining robust to extrinsic noise. This is significant in part due to the nontrivial means by which such systems are derived analytically. WSINDy naturally preserves the Hamiltonian structure by restricting to a trial basis of Hamiltonian vector fields. The methodology is computationally efficient, often requiring only a single trajectory to learn the global reduced Hamiltonian, and avoiding forward solves in the learning process. In this way, we argue that weak-form equation learning is particularly well-suited for Hamiltonian coarse-graining. Using nearly-periodic Hamiltonian systems as a prototypical class of systems with approximate symmetries, we show that WSINDy robustly identifies the correct leading-order system, with dimension reduced by at least two, upon observation of the relevant degrees of freedom. While our main contribution is computational, we also provide a contribution to the literature on averaging theory by proving that first-order averaging at the level of vector fields preserves Hamiltonian structure in nearly-periodic Hamiltonian systems. This provides theoretical justification for our approach as WSINDy’s computations occur at the level of Hamiltonian vector fields. We illustrate the efficacy of our proposed method using physically relevant examples, including coupled oscillator dynamics, the Hénon–Heiles system for stellar motion within a galaxy, and the dynamics of charged particles.
This article introduces a novel systematic methodology for modeling a class of multidimensional linear mechanical systems that directly allows to obtain their infinite-dimensional port-Hamiltonian ...representation. While the approach is tailored to systems governed by specific kinematic assumptions, it encompasses a wide range of models found in current literature, including ℓ-dimensional elasticity models (where ℓ = 1, 2, 3), vibrating strings, torsion in circular bars, classical beam and plate models, among others. The methodology involves formulating the displacement field using primary generalized coordinates via a linear algebraic relation. The non-zero components of the strain tensor are then calculated and expressed using secondary generalized coordinates, enabling the characterization of the skew-adjoint differential operator associated with the port-Hamiltonian representation. By applying Hamilton's principle and employing a specially developed integration by parts formula for the considered class of differential operators, the port-Hamiltonian model is directly obtained, along with the definition of boundary inputs and outputs. To illustrate the methodology, the plate modeling process based on Reddy's third-order shear deformation theory is presented as an example. To the best of our knowledge, this is the first time that a port-Hamiltonian representation of this system is presented in the literature.
•Systematic modeling of multidimensional flexible linear mechanical systems using the port-Hamiltonian framework.•Explicit definition of energy and co-energy variables for this class of port Hamiltonian systems.•Explicit definition of distributed and boundary power conjugated ports suitable for control or interconnection purposes.
We propose new symplectic networks (SympNets) for identifying Hamiltonian systems from data based on a composition of linear, activation and gradient modules. In particular, we define two classes of ...SympNets: the LA-SympNets composed of linear and activation modules, and the G-SympNets composed of gradient modules. Correspondingly, we prove two new universal approximation theorems that demonstrate that SympNets can approximate arbitrary symplectic maps based on appropriate activation functions. We then perform several experiments including the pendulum, double pendulum and three-body problems to investigate the expressivity and the generalization ability of SympNets. The simulation results show that even very small size SympNets can generalize well, and are able to handle both separable and non-separable Hamiltonian systems with data points resulting from short or long time steps. In all the test cases, SympNets outperform the baseline models, and are much faster in training and prediction. We also develop an extended version of SympNets to learn the dynamics from irregularly sampled data. This extended version of SympNets can be thought of as a universal model representing the solution to an arbitrary Hamiltonian system.
We show the Godbillon-Vey invariant arises as a 'restricted Casimir' invariant for three-dimensional ideal fluids associated to a foliation. We compare to a finite-dimensional system, the rattleback, ...where analogous phenomena occur.
Exponential/extended Runge–Kutta–Nyström (ERKN) methods have been validated by numerous theoretical and numerical results to be more efficient and suitable than classical Runge–Kutta–Nyström (RKN) ...methods in dealing with highly oscillatory Hamiltonian systems. Once numerical solutions are required to achieve a prescribed local error tolerance for practical problems, the embedded ERKN pairs with adaptive stepsize are needed. Given the higher efficiency of high-order methods than low-order methods, this paper focuses on high-order embedded ERKN pairs. Using the particular mapping from RKN methods into ERKN methods, we establish an approach to constructing high-order embedded ERKN pairs. By diagonalizing the frequency matrix, an implementation algorithm with less calculation amount than the direct calculation procedure is proposed for high-dimensional problems. In addition, the dispersion and dissipation of the ERKN pairs ERKN4(3) and ERKN8(6) are analyzed. Furthermore, the application of ERKN pairs to an important highly oscillatory problem, i.e., the wave equation prescribed with different boundary conditions is discussed. For the periodic boundary condition, a special implementation algorithm incorporated with FFT techniques is presented, which decreases the calculation amount in one time step from O(N2) to O(NlogN). Finally, numerical results with the FPU problem, the Klein–Gordon equation in the nonrelativistic limit regime, and a two-dimensional wave equation demonstrate the superiority of ERKN pairs over classical RKN pairs.
•An effective approach is presented to establish high-order embedded ERKN pairs.•A special implementation algorithm by diagonalizing the frequency matrix is proposed to reduce the calculation amount.•The application of the embedded ERKN pairs to wave equations prescribed with different boundary conditions is discussed.