2025-12-30 | | Total: 15
This paper investigates the dynamic behavior of a simplified single reed instrument model subject to a stochastic forcing of white noise type when one of its bifurcation parameters (the dimensionless blowing pressure) increases linearly over time and crosses the Hopf bifurcation point of its trivial equilibrium position. The stochastic slow dynamics of the model is first obtained by means of the stochastic averaging method. The resulting averaged system reduces to a non-autonomous one-dimensional Itô stochastic differential equation governing the time evolution of the mouthpiece pressure amplitude. Under relevant approximations the latter is solved analytically treating separately cases where noise can be ignored and cases where it cannot. From that, two analytical expressions of the bifurcation parameter value for which the mouthpiece pressure amplitude gets its initial value back are deduced. These special values of the bifurcation parameter characterize the effective appearance of sound in the instrument and are called deterministic dynamic bifurcation point if the noise can be neglected and stochastic dynamic bifurcation point otherwise. Finally, for illustration and validation purposes, the analytical results are compared with direct numerical integration of the model in both deterministic and stochastic situations. In each considered case, a good agreement is observed between theoretical results and numerical simulations, which validates the proposed analysis.
Covariant Lyapunov vectors (CLVs) are useful in multiple applications, but the optimal time windows needed to accurately compute these vectors are yet unclear. To remedy this, we investigate two methods for determining when to safely terminate the forward and backward transient phases of the CLV computation algorithm by Ginelli et al.~\cite{GinelliEtAl2007} when applied to chaotic orbits of conservative Hamiltonian systems. We perform this investigation for two prototypical Hamiltonian systems, namely the well-known Hénon-Heiles system of two degrees of freedom and a system of three nonlinearly coupled harmonic oscillators having three degrees of freedom, finding very similar results for the two methods and thus recommending the more efficient one. We find that the accuracy of two-dimesnional center subspace computations is significantly reduced when the backward evolution stages of the algorithm are performed over long time intervals. We explain this observation by examining the tangent dynamics of the center subspace wherein CLVs tend to align/anti-align, and we propose an adaptation of the algorithm that improves the accuracy of such computations over long times by preventing this alignment/anti-alignment of CLVs in the center subspace.
Integrable discretizations of the sine-Gordon equation in characteristic (or light-cone) coordinates have been extensively studied after the seminal works of Hirota and Orfanidis in the late 1970s. In contrast, integrable discretizations of the sine-Gordon equation in non-characteristic coordinates have been scarcely studied except the lattice sine-Gordon model proposed by Izergin and Korepin in the early 1980s. In this paper, using the zero-curvature representation, we propose integrable space discretizations of the sine-Gordon equation in three distinct cases of non-characteristic coordinates. For the most interesting case of the sine-Gordon equation in laboratory coordinates, the integrable space discretization is unwieldy; as a remedy, we rewrite the sine-Gordon equation as a two-component evolutionary system and present an aesthetically acceptable space discretization.
In the present work, we review two well-established quasi-continuum models of a Fermi-Pasta-Ulam lattice with Hertzian type potentials, and utilize these two models to approximate the discrete dispersive shock waves (DDSWs) which are numerically observed in the simulation of the lattice. To perform analysis on the various characteristics of the DDSW, we analytically derive the Whitham modulation equations of the two quasi-continuum models, which govern the slowly varying spatial and temporal dynamics of distinct parameters of the periodic solutions. We then perform a very useful reduction of the Whitham modulation system to gain a system of initial-value problems whose solutions can provide important insights on edge features of the DSWs such as their edge speeds. In addition, we also study the numerical rarefaction waves (RWs) of the lattice based on the two quasi-continuum models. In particular, we analytically compute and compare their self-similar solutions with the numerical discrete RW of the lattice. These comparisons made for both DSWs and RWs reveal to be reasonably good, which suggest the impressive performance of both quasi-continuum models.
In this work, we investigate a simple nonequilibrium system with many interconnected, open subsystems, each exchanging a globally conserved resource with an external reserve. The system is represented by a random graph, where nodes represent the subsystems connected through edges. At each time step, a randomly selected node gains a token (i.e, a resource) from the reserve with probability (1-p) or loses a token to the reserve with probability p. When a node loses a token, its neighbors also lose a token each. This asymmetric token exchange breaks the detailed balance. We investigate the steady state behavior of our model for different types of random graphs: graphs without edges, regular graphs, Erdős-Rényi, and Barabási-Albert graphs. In all cases, the system exhibits a sharp, switch-like transition between a token-saturated state and an empty state. When the control parameter p is below a critical threshold, almost all tokens accumulate on the graph. Furthermore, in a non-regular graph, most tokens accumulate or condense on nodes of minimum degree. A slight increase in p beyond the threshold drains almost all the tokens from the graph. This switching transition results from the interplay between drift and the conservation of tokens. However, the position of the critical threshold and the behavior at the transition zone depend on graph topology.
This study introduces amangkurat, an open-source Python library designed for the robust numerical simulation of relativistic scalar field dynamics governed by the nonlinear Klein-Gordon equation in $(1+1)$D spacetime. The software implements a hybrid computational strategy that couples Fourier pseudo-spectral spatial discretization with a symplectic Størmer-Verlet temporal integrator, ensuring both exponential spatial convergence for smooth solutions and long-term preservation of Hamiltonian structure. To optimize performance, the solver incorporates adaptive timestepping based on Courant-Friedrichs-Lewy (CFL) stability criteria and utilizes Just-In-Time (JIT) compilation for parallelized force computation. The library's capabilities are validated across four canonical physical regimes: dispersive linear wave propagation, static topological kink preservation in phi-fourth theory, integrable breather dynamics in the sine-Gordon model, and non-integrable kink-antikink collisions. Beyond standard numerical validation, this work establishes a multi-faceted analysis framework employing information-theoretic entropy metrics (Shannon, Rényi, and Tsallis), kernel density estimation, and phase space reconstruction to quantify the distinct phenomenological signatures of these regimes. Statistical hypothesis testing confirms that these scenarios represent statistically distinguishable dynamical populations. Benchmarks on standard workstation hardware demonstrate that the implementation achieves high computational efficiency, making it a viable platform for exploratory research and education in nonlinear field theory.
Classical chaos theory rests on the notion of universality, whereby disparate dynamical systems share identical scaling laws. Existing universality classes, however, implicitly assume Markovian dynamics. Here, a logistic map endowed with power law memory is used to show that Feigenbaum universality breaks down when temporal correlations decay sufficiently slowly. A critical memory exponent is identified that separates perturbative and memory dominated regimes, demonstrating that long range memory acts as a relevant renormalisation operator and generates a new universality class of chaotic dynamics. The onset of chaos is accompanied by fractional scaling of Lyapunov exponents, in quantitative agreement with analytical predictions. These results establish temporal correlations as a previously unexplored axis of universality in chaotic systems, with implications for physical, biological and geophysical settings where memory effects are intrinsic.
We consider the dispersionless limit of the recently introduced multi-component Pfaff-Toda hierarchy. Its dispersionless version is a set of nonlinear differential equations for the dispersionless limit of logarithm of the tau-function (the F-function). They are obtained as limiting cases of bilinear equations of the Hirota-Miwa type. The analysis of the Pfaff-Toda hierarchy is substantially simplified by using the observation that the full (not only dispersionless) N-component Pfaff-Toda hierarchy is actually equivalent to the 2N-component DKP hierarchy. In the dispersionless limit, there is an elliptic curve built in the structure of the hierarchy, with the elliptic modular parameter being a dynamical variable. This curve can be uniformized by elliptic functions, and in the elliptic parametrization the hierarchy acquires a compact and especially nice form.
We study the $φ^{6}$ model and derive two broad classes of lattice discretizations that admit static, translationally invariant kinks; that is, stationary kink profiles that can be centered at an arbitrary position relative to the lattice. These discretizations are constructed using a one-dimensional map, $φ_{n+1}=F(φ_{n})$, which provides a direct and systematic algorithm for generating such models. Numerical computations for two representative cases show that the discrete kinks do not possess internal modes, consistent with the continuum theory, although an additional high-frequency mode may appear above the phonon band. We also show that generic discretizations of the $φ^{6}$ model do not support static kink solutions. Instead, the resulting dynamics produce auto-traveling and self-accelerating kinks that propagate at the maximal group velocity while continuously emitting radiation.
Topological solitons and tunneling configurations (monopoles, vortices, Skyrmions, sphalerons, and bounces) are characterized by profile functions that encode their physical properties. Derrick's scaling relation provides a global integral constraint on these profiles, but reveals nothing about the local balance between kinetic and potential contributions in different regions. We derive a continuous family of virial identities for O($n$) symmetric configurations, parameterized by an exponent $α$ that controls the radial weighting: negative $α$ emphasizes the core where topological boundary conditions are imposed, large positive $α$ emphasizes the asymptotic tail, and $α= 1$ recovers the classical Derrick relation. The family provides a systematic decomposition of the global constraint into radially-resolved components, with special $α$ values isolating specific mechanisms. For BPS configurations, where the Bogomolny equations imply pointwise equality between kinetic and potential densities, the virial identity is satisfied for all valid $α$. We verify the formalism analytically for the Fubini-Lipatov instanton, BPS monopole, and BPST instanton. Numerical tests on the Coleman bounce and Nielsen-Olesen vortex illustrate how the $α$-dependence of errors distinguishes core from tail inaccuracies: the vortex shows errors growing at negative $α$ (core), while the bounce shows errors growing at positive $α$ (tail). Applications to the electroweak sphaleron, where the Higgs mass explicitly breaks scale invariance, and the hedgehog Skyrmion illustrate the formalism in systems with multiple competing length scales.
We introduce the eigen microstate entropy ($S_{\text{EM}}$), a novel metric of complexity derived from the probabilities of statistically independent eigen microstates. After establishing its scaling behavior in equilibrium systems and demonstrating its utility in critical phenomena (mean spherical, Ising, and Potts models), we apply $S_{\text{EM}}$ to non-equilibrium complex systems. Our analysis reveals a consistent precursor signal: a significant increase in $S_{\text{EM}}$ precedes major phase transitions. Specifically, we observe this entropy rise before biomolecular condensate formation in liquid-liquid phase separation in living cells and months ahead of El Niño events. These findings position $S_{\text{EM}}$ as a general framework for detecting and interpreting phase transitions in non-equilibrium systems.
We report on the presence of families of exact solutions for a complex scalar field that behaves according to the rules of discrete $Z_N$ symmetry. Since the family of models is exactly solved, the results appear to be of interest to integrability, to build junctions and networks of localized structures and to describe scalar dark matter in high energy physics.
Westerly wind bursts (WWBs) have long been known to have a major impact on the development of El Niño events. In particular, they amplify these events, with stronger events associated with a higher number of WWBs. We further find indications that WWBs lead to a more monotonically increasing evolution of warming events. We consider here a noise-driven recharge oscillator model of ENSO. Commonly, WWBs are represented by a state-dependent Gaussian noise which naturally reproduces the amplification of warm events. However, we show that many properties of WWBs and their effects on sea surface temperature (SST) are not well captured by such Gaussian noise. Instead, we show that conditional additive and multiplicative (CAM) noise presents a promising alternative. In addition to recovering the sporadic nature of WWBs, CAM noise leads to an asymmetry between El Niño and La Niña events without the need for deterministic nonlinearities. Furthermore, CAM noise generates a more monotonic increase of extreme warming events with a higher frequency of WWBs accompanying the largest events. This suggests that extreme warm events are better modelled by CAM noise. To cover the full spectrum of warm events we propose a conditional noise model in which the wind stress is modelled by additive Gaussian noise for sufficiently small SSTs and by additive CAM noise once the SST exceeds a certain threshold. We show that this conditional noise model captures the observed properties of WWBs reasonably well.
The quantum internet is a rapidly developing technological reality, yet, it remains unclear what kind of quantum network structures might emerge. Since indirect quantum communication is already feasible and preserves absolute security of the communication channel, a new node joining the quantum network does not need to connect directly to its desired target. Instead, in our proposed quantum preferential attachment model, it uniformly randomly connects to any node within the proximity of the target, including, but not restricted to, the target itself. This local flexibility is found to qualitatively change the global network behavior, leading to two distinct classes of complex network architectures, both of which are small-world, but neither of which is scale-free. Our numerical findings are supported by rigorous analytic results, in a framework that incorporates quantum and classical variants of preferential attachment in a unified phase diagram. Besides quantum networks, we expect that our results will have broad implications for classical scenarios where there is flexibility in establishing new connections.
We revisit the canonical quantization to assess the spectrum of the modified Emden equation $\ddot{x} + kx\dot{x} + ω^2 x + \frac{k^2}{9}x^3 = 0$, which is an isochronous case of the Liénard-Kukles equation. While its classical isochronicity and canonical quantization, leading to polynomial solutions with an exactly-equispaced spectrum have been discussed earlier, including in the recent paper [Int. J. Theor. Phys. 64, 212 (2025)], the present study focuses on the quantization of its branched Hamiltonians. For small $k$, we show numerically that the resulting energy spectrum is no longer perfectly harmonic but only approximately equispaced, exhibiting quasi-harmonic behavior characterized by deviations from uniform spacing. Our numerical results are precisely validated by analytical calculations based on perturbation theory.