Nonlinear Perturbations Research Articles

Spin cones in random-field XY models.

We determine the arrangement of spins in the ground state of the XY model with quenched, random fields, on a fully connected graph. Two types of disordered fields are considered, namely, randomly oriented magnetic fields and randomly oriented crystal fields. Orientations are chosen from a uniformly isotropic distribution, but disorder fluctuations in each realization of a finite system lead to a breaking of rotational symmetry. The result is an interesting pattern of spin orientations found by solving a system of coupled, nonlinear equationswithin perturbation theory and also by exact numerical continuation. All spins lie within a cone for small enough ratio of field to coupling strength, with an interesting distribution of spin orientations, with peaks at the cone edges. The orientation of the cone depends strongly on the realization of disorder, but the opening angle does not. In the case of random magnetic fields, the cone angle widens as the ratio increases till a critical value at which there is a first-order phase transition and the cone disappears. With random crystal fields, there is no phase transition and the cone angle approaches 180^{∘} for large values of the ratio. At finite low temperatures, Monte Carlo simulations show that the formation of a cone and its subsequent alignment along the equilibrium direction occur on two different timescales.

Active suspension hierarchical control with parameter uncertainty and external disturbance of electro-hydraulic actuators

In order to improve the ride comfort and handling stability of the electro-hydraulic active suspension system, a hierarchical control strategy is proposed. For the active suspension system with body mass uncertainty and safety constraints, an enhanced constraint adaptive backstepping controller is designed to generate the target force of body vertical motion. When dealing with constraints, nonlinear filters are introduced, backstepping technology and quadratic Lyapunov function are used to combine the control target and domain constraints, and the vertical displacement of the body and the suspension deflection control index are synthesized into a single controlled variable, so that the control variable converging to zero meets the suspension dynamic displacement limit. The ride comfort and safety of the active suspension system are improved. Secondly, stability analysis is conducted on the zero dynamic error system to ensure that all safety performance indicators are bounded. The effects of external disturbances, parameter uncertainties and modelling errors on the force tracking accuracy of asymmetric cylinder electrohydraulic systems are addressed. A backstepping dynamic surface control method based on a nonlinear perturbation observer is proposed, which estimates the composite perturbation terms such as modelling error and external perturbation, and uses the estimated value of the nonlinear observer to design a feedforward compensating control term to offset the effect of the composite perturbation on the system performance. In addition, the dynamic surface control method is used to calculate the derivative of the target force input, which avoids the derivative explosion phenomenon and reduces the computational complexity of the system. Simulation test results show that this method reduces the computational complexity of the system and improves the control performance. And under random and bumpy road conditions, the root-mean-square value of sprung mass acceleration performance index is reduced by 48.354 % and 72.677 % compared with passive suspension, which improves the ride comfort of the vehicle.

The modified cubic Benjamin–Ono equation describing internal solitary waves in the deep ocean and its related properties

In this article, we investigate the propagation of internal solitary waves in deep ocean. Based on the principles of nonlinear theory, perturbation expansion, and multi-scale analysis, a time-dependent modified cubic Benjamin–Ono (mCBO) equation is derived to describe internal solitary waves in the deep ocean with stronger nonlinearity. When the dispersive term ∂3f∂X3 vanishes, the mCBO equation transforms into the cubic BO equation. Similarly, when the dispersive term ∂3f∂X3 becomes zero and the nonlinear term ∂f3∂X degenerates into ∂f2∂X, the mCBO equation reduces to the BO equation. Furthermore, if the integral term ∂2∂X2ℵ(f) disappears, it simplifies to the mKdV equation. To gain deeper insight into the characteristics of solitary waves, conservation of mass and momentum associated with them are discussed. By employing Hirota's bilinear method, we obtain soliton solutions for the mCBO equation and subsequently investigate interactions between two solitary waves with different directions, leading to the occurrence of important events such as rogue waves and Mach reflections. Additionally, we explore how certain parameters influence Mach stem while drawing meaningful conclusions. Our discoveries reveal the complex dynamics of internal solitary waves within the deep ocean and contribute to a broader understanding of nonlinear wave phenomena.

Hyers–Ulam stability of integral equations with infinite delay

Integral equations with infinite delay are considered as functional equations in a Banach space. Two types of Hyers–Ulam stability criteria are established. First, it is shown that a linear autonomous equation is Hyers–Ulam stable if and only if it has no characteristic value with zero real part. Second, it is proved that the Hyers–Ulam stability of a linear autonomous equation is preserved under sufficiently small nonlinear perturbations. The proofs are based on a recently developed decomposition theory of linear integral equations with infinite delay.

The effects of nonlinear perturbation terms on comparison principles for the p-Laplacian

In this paper, we investigate various comparison principles for quasilinear elliptic equations of [Formula: see text]-Laplace type with lower-order terms that depend on the solution and its gradient. More specifically, we study comparison principles for equations of the following form: [Formula: see text] where [Formula: see text] is the [Formula: see text]-Laplace operator with [Formula: see text], and [Formula: see text] is a continuous function that satisfies a structure condition. Many of these results lead to comparison principles for the model equations [Formula: see text] where [Formula: see text] are non-decreasing and [Formula: see text]. Our results either improve or complement those that appear in the literature.

Verification of Long‐Term Ensemble Evapotranspiration Hindcast Using a Conditional Nonlinear Optimal Parameter Perturbation Ensemble Prediction Method on the Tibetan Plateau

AbstractIn this research, a conditional nonlinear optimal perturbation related to parameters (CNOP‐P) method is employed to propose an ensemble prediction method titled as the conditional nonlinear optimal perturbation related to parameters ensemble prediction (CNOP‐PEP) method. Within the CNOP‐PEP method, all ensemble members are generated by the CNOP‐P method. To explore the operability and validity of the CNOP‐PEP method, long‐term evapotranspiration (ET) hindcast skill is evaluated at 13 sites on the Tibetan Plateau (TP) during the period of 2001–2018 with a land surface model. Two traditional ensemble prediction methods (the stochastically perturbed parameters (SPP) method and the one‐at‐a‐time (OAT) method) are also applied to assess the ET hindcast skills. The numerical results indicate that the CNOP‐PEP method shows the best hindcast skill for estimated ET at 13 stations on the TP among the three methods (CNOP‐PEP, OAT, and SPP methods) during the period of 2001–2018. This suggests that the CNOP‐PEP method is helpful and effective for long‐term and interannual hydrological predictions, such as long‐term and interannual ET predictions. The numerical results also indicate that ensemble members with certain special properties of fast‐growing types should be selected to implement ensemble prediction compared to these ordinary and traditional samples. The CNOP‐PEP method could supply special samples to improve the ensemble prediction and hindcast skill of hydrological cycle.

Deformation in the wrinkle–crease transformation

We investigate deformation in the wrinkle–crease transformation to reveal the missing link between the wrinkling and creasing instabilities, which occur on the surface of elastomers under compression. Plane-strain finite element (FE) analysis is performed by combining a nonlinear perturbation approach to find two bifurcated solutions for the flat surface in a metastable state. The first solution is the stable deformation path for crease evolution, and the second solution is the previously unknown unstable deformation path from wrinkling to creasing. The two deformation paths are not independent, and they are connected at the critical strain for creasing εC. The resulting deformation at εC acts as a precursive deformation to cause creasing. Although the unstable path is initiated by the wrinkling instability as a bifurcation, the creasing instability (i.e., εC) is interpreted as the singular point at which the unstable path switches to the stable path. A chain reaction model is developed to describe the mechanism and deformation of the wrinkle–crease transformation, and it shows good agreement with the FE prediction within an inevitable limitation of FEs.

Near-term distributed quantum computation using mean-field corrections and auxiliary qubits

Distributed quantum computation is often proposed to increase the scalability of quantum hardware, as it reduces cooperative noise and requisite connectivity by sharing quantum information between distant quantum devices. However, such exchange of quantum information itself poses unique engineering challenges, requiring high gate fidelity and costly non-local operations. To mitigate this, we propose near-term distributed quantum computing, focusing on approximate approaches that involve limited information transfer and conservative entanglement production. We first devise an approximate distributed computing scheme for the time evolution of quantum systems split across any combination of classical and quantum devices. Our procedure harnesses mean-field corrections and auxiliary qubits to link two or more devices classically, optimally encoding the auxiliary qubits to both minimize short-time evolution error and extend the approximate scheme’s performance to longer evolution times. We then expand the scheme to include limited quantum information transfer through selective qubit shuffling or teleportation, broadening our method’s applicability and boosting its performance. Finally, we build upon these concepts to produce an approximate circuit-cutting technique for the fragmented pre-training of variational quantum algorithms. To characterize our technique, we introduce a non-linear perturbation theory that discerns the critical role of our mean-field corrections in optimization and may be suitable for analyzing other non-linear quantum techniques. This fragmented pre-training is remarkably successful, reducing algorithmic error by orders of magnitude while requiring fewer iterations.

Existence of solutions for anisotropic parabolic Ni-Serrin type equations originated from a capillary phenomena with nonstandard growth nonlinearity

We consider an initial boundary value problem for a class of anisotropic parabolic Ni-Serrin type equations with nonstandard nonlinearity in a bounded smooth domain with homogeneous Dirichlet boundary condition. Because the nonlinear perturbation leads to difficulties (it does not have a definite sign) in obtaining a priori estimates in the energy method, we had to modify the Tartar method significantly. Under suitable assumptions, we obtain the global existence, decay and extinction of solutions.

Oscillons in gapless theories

We show that large scale oscillons, i.e., quasiperiodic, long-living particlelike solutions, may exist in massless theories, too. Their existence is explained using an effective (smeared) mass threshold which takes into account nonlinear (finite) perturbations. Published by the American Physical Society 2024

Finite‐time mixed ℋ∞$$ {\mathscr{H}}_{\infty } $$ and passive filtering for nonlinear singular systems under fading channels

AbstractThis article focuses on the finite‐time mixed and passive filtering problem for discrete‐time singular systems with randomly occurring nonlinear perturbations. The main purpose is to design a filter, which ensures the singular error system is stochastically finite‐time bounded and satisfies mixed and passive performance. By constructing an augmented matrix, the properties of matrix determinant and rank are used to prove that the singular system is regular and causal, Finsler's lemma and Projection lemma are used in the design process to acquire additional slack variable matrices, thereby enhancing the solution space with extra degrees of freedom. Then, the unknown parameters of the filter are obtained by solving the less conservative linear matrix inequalities. Finally, the feasibility of the proposed method is verified through a common numerical example and by controlling a DC motor for an inverted pendulum.

Systematic search for islets of stability in the standard map for large parameter values

In the seminal paper, Chirikov (Phys Rep 52:263–379, 1979) showed that the standard map does not exhibit a boundary to chaos, but rather that there are small islands (“islets”) of stability for arbitrarily large values of the nonlinear perturbation. In this context, he established that the area of the islets in the phase space and the range of parameter values where they exist should decay following power laws with exponents -2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$-2$$\\end{document} and -1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$-1$$\\end{document}, respectively. In this paper, we carry out a systematic numerical search for islets of stability and we show that the power laws predicted by Chirikov hold. Furthermore, we use high-resolution 3D islets to reveal that the islets’ volume decays following a similar power law with exponent -3\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$-3$$\\end{document}.

Modified Nested Saturated Control for Uncertain Multiple Integrators With High-Order Nonlinear Perturbation.

This article studies the global asymptotic stability problem of multiple integrators subject to uncertain parameters and high-order nonlinear perturbation. By introducing a modified saturation function with positive tunable parameters, three types of modified nested saturated controllers with different parameter assignments are proposed to handle the uncertain integrators with high-order nonlinear perturbation. The specific mechanism of the modified nested saturated controllers with different parameter assignments, in which control parameters appear before states, saturation functions, or saturation levels, is comprehensively analyzed. With the help of a modified saturation function, the global convergence domain of each state can be kept invariable in the saturation reduction analysis. Meanwhile, a large saturation level is obtained by adjusting positive tunable parameters. Combining with derivatives calculation and M -matrix-based comparison theory, proposed saturated control schemes are utilized to handle the considered system. Finally, an active magnetic bearing (AMB) system is revisited as an uncertain mechanical system and used to verify the feasibility of the proposed control designs.

On nonlocal Fokker–Planck equations with nonlinear force fields and perturbations

On nonlocal Fokker–Planck equations with nonlinear force fields and perturbations

The sampling method for optimal precursors of El Niño–Southern Oscillation events

Abstract. The El Niño–Southern Oscillation (ENSO) is a significant climate phenomenon that appears periodically in the tropical Pacific. The intermediate coupled ocean–atmosphere Zebiak–Cane (ZC) model is the first and classical one designed to numerically forecast the ENSO events. Traditionally, the conditional nonlinear optimal perturbation (CNOP) approach has been used to capture optimal precursors in practice. In this paper, based on state-of-the-art statistical machine learning techniques1, we investigate the sampling algorithm proposed in Shi and Sun (2023) to obtain optimal precursors via the CNOP approach in the ZC model. For the ZC model, or more generally, the numerical models with a large number O(104−105) of degrees of freedom, the numerical performance, regardless of the statically spatial patterns and the dynamical nonlinear time evolution behaviors as well as the corresponding quantities and indices, shows the high efficiency of the sampling method compared to the traditional adjoint method. The sampling algorithm does not only reduce the gradient (first-order information) to the objective function value (zeroth-order information) but also avoids the use of the adjoint model, which is hard to develop in the coupled ocean–atmosphere models and the parameterization models. In addition, based on the key characteristic that the samples are independently and identically distributed, we can implement the sampling algorithm by parallel computation to shorten the computation time. Meanwhile, we also show in the numerical experiments that the important features of optimal precursors can still be captured even when the number of samples is reduced sharply.

Dual sliding mode control of linear maglev synchronous motor based on novel extended state observer

In order to improve the convergence speed and disturbance immunity of the linear maglev synchronous motor (LMLSM) maglev system, this paper proposes a dual sliding mode control strategy based on a novel expanded state observer (NESO). Firstly, a dual terminal sliding mode (DTSM) surface is designed, which exhibits the desired dynamic and steady state performance. The problem of the TSM singularity is solved by the dual sliding mode technique, which turns the control signal into a continuous switching control through the filtering function and effectively reduces the jitter problem. Since the maglev system is susceptible to uncertainties such as nonlinearity, strong coupling and external perturbations, which affect its tracking accuracy, a continuous smooth ESO is designed to estimate and compensate for these uncertainty perturbations, thus solving the problem that the traditional ESO is continuous but not derivable at the switching point. Then, it is proved that the DTSM control strategy based on NESO converges to the equilibrium point in finite time by Lyapunov function. Finally, the effectiveness and superiority of the proposed strategy are verified on the LMLSM maglev system platform.

Linear versus nonlinear modeling of black hole ringdowns

The ringdown (RD) phase of gravitational waves is of prime interest for testing general relativity (GR). The modeling of the linear quasinormal modes (QNMs) within the Kerr spectrum—or with agnostic parametrized deviations to that GR spectrum—has become ordinary; however, specific attention has recently emerged to calibrate the effects of nonlinear perturbations for the predominant quadrupolar l=2, m=2 mode. In this paper, we test the performance of a few nonlinear toy models and of the nonlinear inspiral-merger-ringdown (IMR) model IMRPhenomD for faithfully representing the RD regime and we compare them with the results obtained using linear solutions as sums of QNM tones. Using several quasicircular, nonprecessing numerical waveforms, we fit the dominant l=2, m=2 mode of the strain, and we assess the results in terms of both the Bayes factor and the inferred posterior distributions for the mass and spin of the final black hole. We find that the nonlinear models can be comparable or preferred over the linear QNM-only solutions when the analysis is performed from the peak of the strain, especially at high signal-to-noise ratios consistent with third-generation observatories. Since the calibration of the tones’ relative amplitudes and phases in high-overtone models to the progenitor parameters is still missing, or even not achievable, we consider the use of nonlinear models to be more pertinent for performing confident tests of general relativity based on the RD regime starting from early times. Published by the American Physical Society 2024

Observation of fixed lines induced by a nonlinear resonance in the CERN Super Proton Synchrotron

The motion of systems with linear restoring forces and recurring nonlinear perturbations is of central importance in physics. When a system’s natural oscillation frequencies and the frequency of the nonlinear restoring forces satisfy certain algebraic relations, the dynamics become resonant. In accelerator physics, an understanding of resonances and nonlinear dynamics is crucial for avoiding the loss of beam particles. Here we confirm the theoretical prediction of the dynamics for a single two-dimensional coupled resonance by observing so-called fixed lines. Specifically, we use the CERN Super Proton Synchrotron to measure the position of a particle beam at discrete locations around the accelerator. These measurements allow us to construct the Poincaré surface of section, which captures the main features of the dynamics in a periodic system. In our setting, any resonant particle passing through the Poincaré surface of section lies on a curve embedded in a four-dimensional phase space, the fixed line. These findings are relevant for mitigating beam degradation and thus for achieving high-intensity and high-brightness beams, as required for both current and future accelerator projects.

Cosmic-Eν: An- emulator for the non-linear neutrino power spectrum

ABSTRACT Cosmology is poised to measure the neutrino mass sum Mν and has identified several smaller-scale observables sensitive to neutrinos, necessitating accurate predictions of neutrino clustering over a wide range of length scales. The FlowsForTheMasses non-linear perturbation theory for the the massive neutrino power spectrum, $\Delta ^2_\nu (k)$, agrees with its companion N-body simulation at the $10~{{\ \rm per\ cent}}-15~{{\ \rm per\ cent}}$ level for k ≤ 1 h Mpc−1. Building upon the Mira-Titan IV emulator for the cold matter, we use FlowsForTheMasses to construct an emulator for $\Delta ^2_\nu (k)$, Cosmic-Eν, which covers a large range of cosmological parameters and neutrino fractions Ων, 0h2 ≤ 0.01 (Mν ≤ 0.93 eV). Consistent with FlowsForTheMasses at the 3.5 per cent level, it returns a power spectrum in milliseconds. Ranking the neutrinos by initial momenta, we also emulate the power spectra of momentum deciles, providing information about their perturbed distribution function. Comparing a Mν = 0.15 eV model to a wide range of N-body simulation methods, we find agreement to 3 per cent for k ≤ 3kFS = 0.17 h Mpc−1 and to 19 per cent for k ≤ 0.4 h Mpc−1. We find that the enhancement factor, the ratio of $\Delta ^2_\nu (k)$ to its linear-response equivalent, is most strongly correlated with Ων, 0h2, and also with the clustering amplitude σ8. Furthermore, non-linearities enhance the free-streaming-limit scaling $\partial \log (\Delta ^2_\nu /\Delta ^2_{\rm m}) / \partial \log (M_\nu)$ beyond its linear value of 4, increasing the Mν-sensitivity of the small-scale neutrino density.

On Convergence Rate Bounds for a Class of Nonlinear Markov Chains

A new approach is developed for evaluating the convergence rate for nonlinear Markov chains (MC) based on the recently developed spectral radius technique of Markovian coupling for linear MC and the idea of small nonlinear perturbations of linear MC. The method further enhances recent advances in the problem of convergence for such models.