Papers

A two-component Bose–Einstein condensate can host a pseudospin Hopf texture coupled to the viscous normal flow of its thermal cloud. Motivated by OpenAI's reported finite-time Navier–Stokes singular solution, we evolve such a texture in a smooth flow reproducing the geometry and similarity scaling of the collapsing axisymmetric core. The evolution creates and annihilates preimage pairs while preserving the Hopf charge. The first one-to-three transition exhibits the square-root opening characteristic of a Thom $A_2$ fold. A subsequent three-to-five transition precedes a burst-like regime of recurrent fold catastrophes and near-caustic activity.

The logistic map is a quintessential model in the study of low-dimensional chaos. High-dimensional chaos, on the other hand, presents itself in disordered systems with many interacting and heterogeneously coupled components. Here, we formulate a system of many logistic maps, interacting through disordered couplings. Using a combination of dynamic mean-field theory, random matrix theory and numerical simulations, we show that even the smallest amount of disorder can remove the period-doubling cascade in the conventional logistic map. Instead we find a transition to high-dimensional chaos, marked by an oscillatory instability not previously reported for disordered systems. We also show that with sufficiently strong homogeneous coupling between the maps one recovers elements of the conventional period-doubling cascade. Our findings indicate that well-known phenomena in dynamical systems can be fragile in the face of disorder. At the same time, new phenomena emerge that are neither found in simple low-dimensional dynamics nor in high-dimensional disordered systems.

We study rolls and homoclinic snaking in a variation of the one-dimensional Swift-Hohenberg equation, whose standard forms are prototypical order-parameter models for pattern formation in the sciences. Motivated by classes of differential equation models that involve continuous non-smooth low order nonlinear terms, we replace the standard quadratic-cubic nonlinearity by $\nu |u|^\alpha-u^3$, $\alpha \in [1,2]$ with $\nu > 0$. In the vicinity of zero, for $\alpha<2$ this nonlinearity falls outside the scope of classical Taylor expansion and bifurcation analysis. Our partially analytical and partially numerical results highlight that the non-smooth term modifies the criticality of pattern-forming bifurcations and alters the associated branches of solutions. In particular, $\alpha\in(1,2)$ implies subcriticality of roll bifurcations for any $\nu>0$. At $\alpha=1$ differentiability is lost, which has a strong impact on the bifurcations of sign-changing rolls including the disappearance of homoclinic snaking. Homoclinic snaking thus emerges non-smoothly as $\alpha$ increases from $\alpha=1$, and persists when retaining an additional, e.g., quadratic term.

We investigate the equilibrium thermodynamics of a Dyson gas in connection with a set of integrable statistical mechanical observables satisfying the Toda Lattice hierarchy. We prove that in the thermodynamic limit, the integrable observables are state functions satisfying a set algebraic equations of state in closed form, obtained from direct integration of the Toda Lattice hierarchy in the continuum limit. We then explore the connection between regularity and critical behaviour of the state functions and the Dyson gas particle statistics via Monte Carlo simulations. We show that the properties of the integrable observables, such as regularity, multivaluedness, cusp singularities, carry information on the macroscopic particle statistics and its qualitative changes but with some limitations.

We compute the effective six-wave interaction in spatial and temporal Hamiltonian Dysthe equations and compare it with the complete degree-six coefficient of the deep-water Euler Hamiltonian. On every compact subset $K$ of a generic quadratic-dispersion resonance patch, exact resonances persist under a small cubic correction to the dispersion, and \[ W_6\bigl(\Phi_{\lambda,\varepsilon}(z)\bigr) =\varepsilon\frac{4BD}{aL(z)}+O_K(\varepsilon^2) \] uniformly in the resonance point $z$ and the deformation parameter $\lambda$. Here $L$ is the spectral diameter and $B,D$ are the local cubic and mean-flow coefficients. Thus $BD\ne0$ gives nonzero scattering on relatively open resonance sets. At the Craig–Guyenne–Sulem temporal coefficients, the corresponding exact Euler resonances satisfy \[ \frac{\mu}{4}W_6^{\rm E}\bigl(1+\mu\Psi_\mu(z)\bigr) =-\frac{8}{L(z)}+O_K(\mu), \] where $\mu$ is the relative bandwidth. The leading term comes from nine contractions of quartic vertices; every other connected degree-six contribution is uniformly bounded. We classify the coefficient values for which the generic six-wave coefficient vanishes identically; in the temporal case, this also gives cancellation on every non-permutation four-wave resonance. Outside these sets, the Zakharov–Schulman condition restricts every regular $C^3$ quadratic leading symbol to $\operatorname{span}\{1,k,\omega(k)\}$, the span of the symbols of wave action, momentum, and linear energy. Consequently the physical Dysthe models cannot support an inverse-scattering hierarchy with infinitely many linearly independent regular $C^3$ quadratic leading symbols. The cancellation sets contain known nonlinear Schrödinger, mixed Chen–Lee–Liu, Calogero–Moser derivative nonlinear Schrödinger, and Hirota representatives.

We resolve an existing question concerning dynamical and structural stability of the E {\times} B staircases in magnetically confined fusion plasma. The problem has stirred considerable interest in the literature due to a unique combination of its scientific and practical aspects. Here, we report a new approach to the problem by combining the paradigmatic Hasegawa-Wakatani model of electrostatic drift-wave turbulence with an idea of interaction renormalization, imported from quantum field theory. It is shown based on the Schrödinger-Newton equation with renormalized interaction constant that the plasma staircase is robust and stable as soon as the nonadiabaticity parameter, {\delta}, goes above a certain critical value, and is unstable otherwise. This critical value is an exact result of the interaction model and is found to be {\delta}c = {\pi}. For {\delta} < {\delta}c , the unstable system expands in the radial direction, resulting in a radial second moment that grows with time t as t^(2/3). The results obtained bear important implications with regard to the implementation of arrays of semipermeable transport barriers in tokamaks, suggesting a threshold condition for their generation.

The hypersensitivity and turnpike phenomena in the optimal control of an inverted pendulum are investigated from a dynamical-systems perspective. We show that, for a fixed terminal time and a fixed terminal state optimal control problem, (1) the hypersensitivity originates from the fractal structure of the set of initial adjoint variables in the associated Hamiltonian dynamics, (2) the turnpike arises from slow dynamics in the vicinity of a degenerate center manifold, and (3) the escape channels are formed by normally hyperbolic invariant manifolds (NHIMs). As a consequence, small perturbations in the initial adjoint variables lead to qualitatively distinct extremal trajectories, resulting in severe numerical instability in trajectory optimization. Both the fractal structure and the invariant sets are characterized numerically and analytically.

Equilibria of the Hamiltonian system associated with Pontryagin's minimum principle exhibit an exact doubling of the natural spectrum and, under a simple pairing condition, a Jordan block at every simple purely imaginary eigenvalue. Consequently, the classical Lyapunov Center Theorem does not apply to the augmented system, and no periodic orbit with nonzero optimal control bifurcates at linear order. We establish this mechanism in general and show that an optimal-control-induced periodic family emerges at cubic order in a Lindstedt–Poincaré expansion. The mechanism is illustrated in closed form for the pendulum and evaluated numerically for the planar $L_2$ equilibrium of Hill's restricted three-body problem, where the third-order approximation is validated against an independently computed family of periodic orbits.

At an $O(2)$-equivariant Hopf bifurcation, the real parts $\zeta_R$ and $\xi_R$ of the two cubic normal-form coefficients govern selection between traveling and standing waves. In a fourth-order regularized $p$-system, $\xi_R$ was observed to vanish and the cubic nonlinearity to shift only frequencies. We show that both observations are structural: variational nonlinearities enforce them even when the linear part is dissipative and non-Hamiltonian. We consider a class of two-component PDEs on the circle whose quadratic and cubic nonlinearities are generated by reflection-invariant local Hamiltonians of the fields and finitely many of their spatial derivatives, through a common constant-coefficient Poisson operator. Then $\xi_R=0$ for every admissible linear part, while $\zeta_R$ depends only on the quadratic nonlinearity. For variational nonlinearities the direct cubic terms are purely imaginary. The mixed cascade term loses its real part because its second harmonic has zero temporal frequency and the resolvent keeps its parity, while the self-cascade at the doubled critical frequency keeps a real part. When $\zeta_R\ne0$, the bifurcating traveling waves are saddles, while the standing waves are supercritical and orbitally asymptotically stable for $\zeta_R<0$ and subcritical and unstable for $\zeta_R>0$. When the nonlinearities are the Poisson operator applied to polynomials in the fields alone, and its symbol does not vanish at the critical and second-harmonic wavenumbers, cancellation for every admissible linear part conversely characterizes variationality. In a dissipatively regularized Boussinesq family, as dissipation tends to zero, $\zeta_R$ vanishes linearly off resonance and diverges inversely at the 2:1 resonance. The Hamiltonian limit of the cubic coefficient is therefore singular.

We investigate the discrete-time version of the Kuramoto model with phase lag, which comprises globally-coupled phase oscillators of distributed frequencies that are evolving under a nonlinear map. In the continuum limit of an infinite number of oscillators ($N\to \infty$), we derive the exact Frobenius-Perron equation for the time evolution of the single-oscillator probability density, and study linear stability of the incoherent state. Instability signals onset of synchronization. The corresponding synchronization threshold is obtained analytically for the case of a Lorentzian distribution of the oscillator frequencies. The threshold differs from that of the continuous-time Kuramoto model, reflecting the fundamentally different stability conditions for discrete-time maps and continuous-time flows. Beyond synchronization threshold, we observe several interesting nonlinear phenomena: Unlike the classical Kuramoto model, the discrete-time version exhibits periodic and chaotic states. Numerical simulations of the finite-$N$ system confirm the analytical prediction for the synchronization threshold, while highlighting breakdown of the celebrated Ott-Antonsen ansatz invoked to conveniently study the continuous-time Kuramoto model in terms of a low-dimensional description.

Empirical regularities described by parametrized families of probability distributions can shape how complex systems are interpreted and explained. A recurring problem is to determine the extent to which these descriptions are supported by observations, and how they can be used to characterize the underlying system. Here, we introduce a general framework for constructing constrained surrogates capable of representing arbitrary families of continuous distributions. We demonstrate its utility by providing explicit algorithms for several canonical distributions in statistical physics, including continuous powerlaw, lognormal, exponential, and truncated Gaussian distributions. These constrained surrogates enable unbiased estimation and accurate tests of model validity based on arbitrary (physically meaningful) sample statistics, without requiring specification or estimation of model parameters. We show how these theoretical advantages translate into more informative and nuanced conclusions in applications ranging from bushfire magnitudes and word usage in texts to brain connectomes and city populations.

Which cellular-automaton diagonal families form bases in every finite window? For canonical polynomial lifts of elementary rules, two truth-table bits determine triangularity, and units on the matrix diagonal determine invertibility. Exactly 24 rules give universal binary bases; all remain universal over every modulus. Among triangular binary coordinate maps, the Pascal transform is uniquely characterized by converting OR convolution into pointwise multiplication, while increment becomes strict prefix summation. Explicit inverses and coordinate comparisons distinguish sparsity from evaluation cost. A Rule 30 polynomial construction gives Fibonacci bounds on interpolation order and prime-modulus periods. Exact additive periods anchor a finite census modulo two and three. These results separate all-window basis classification from optimization of a representation and from period patterns observed in finite windows.

The single-seed Rule 30 cone has unbounded least diagonal periods in both directions. On the finite-support side, first contact fixes the least Newton index and an integer polynomial lift gives a Fibonacci ceiling. For reflected Rule 86, a backward map on periodic tail profiles retains the zero boundary and bounds periods independently of transients. At depth $m$, the Rule 30 period is at least $\lfloor m/2\rfloor+1$; for natural Rule 86 cuts through depth $m$, the largest period is at least $\tfrac12\log_2(m+2)$. This proves infinitely many reset/integration period doublings. The physical Rule 135 orbit is a translated complement of Rule 30 after its first step, yielding finite or cofinite supports. A separately initialized auxiliary recurrence has exact index-by-index stabilization. Residue tables and block formulas describe the remaining structure and permit evaluation from supplied representations. Construction and minimization remain separate costs.

A wide range of physical and biological systems are adaptive networks, in which the dynamics of the nodes and of the edges connecting them co-evolve. Mean-field reductions of such systems typically track only the average coupling strength, and therefore cannot determine when the coupling stays effectively homogeneous and when structured connectivity emerges. Here, we present a second-order moment closure that allows us to derive mean-field equations for the coupling-weight variance in networks of heterogeneous phase oscillators with adaptive coupling, starting from uniform coupling weights. In agreement with network simulations, we find a nonlinear, non-monotonic dependence of the relative weight variance on the oscillator heterogeneity that is mediated by the phase coherence. Moreover, we find that deviations from global coupling strongly depend on an interaction between the oscillator heterogeneity and the adaptation rule. Whereas symmetric adaptation causes a strongly coupled core of coherent oscillators to emerge and creates a bistable regime that is absent without adaptation, antisymmetric adaptation leads to antisymmetric coupling within the same core, thereby destabilizing it. Our equations therefore delineate the regimes in which adaptive networks behave like globally coupled systems from those in which more complex coupling patterns form.

The hat and spectre monotiles, discovered in 2023, tile the plane only aperiodically; no cellular automaton dynamics on these tilings has previously been reported. Cellular automata are studied here on patches generated by finite-state transducers, so that every experiment regenerates deterministically from a small record. Within edge-adjacency semi-totalistic rules, exhaustive and evolutionary searches find only mortal travelers: gliders are absent. Guided by a reproduction of the known Penrose-tiling glider, the rule space is extended to vertex neighborhoods and to priority-table rules whose non-quiescent states are visible to neighbors. Evolutionary search then discovers gliders on both monotilings; tracked by a sliding window that regenerates the patch along the flight, they travel one million rings at constant speed and heading. All headings are quantized, to millidegrees, onto a six-spoke compass - the fast axes of the tiling's graph metric. An ablation shows both rule-space extensions are individually necessary. All results replay exactly in an accompanying interactive essay.

We present an equation-free method for maintaining a dynamical system operating in a desired regime, even as system parameters are disturbed in such a way that qualitatively different dynamics would emerge. This method allows for real-time identification of changes made to a system's evolution equations induced by the parameter changes. If the equations governing the original system's dynamics are known, our method can allow one to determine the values of the parameters in real time. We illustrate our method with numerical examples using a reservoir computer control implementation.

Estimating the dynamic aperture of circular particle accelerators typically requires expensive long-term tracking. Using classical homology, we construct a two-dimensional theoretical classification of regular and chaotic orbits and use persistent homology for its numerical implementation on short orbit segments. In four dimensions, we introduce topological chaos indicators based on persistent homology and use normal-form theory to apply the classification to two symplectic-plane projections. The methods reproduce the expected organization of regular and chaotic motion in several canonical Hamiltonian systems. Applied to symplectified one-turn maps of a nonlinear lattice model of the Integrable Optics Test Accelerator, the results obtained from the first 200 turns identify the central connected stability region and its surrounding transition to chaos. They are consistent with the short-term dynamic aperture of the electron configuration obtained from $10^4$ turns and the long-term dynamic aperture of the proton configuration obtained from $10^8$ turns. For the configurations studied, persistent homology of short orbit segments predicts the dynamic aperture obtained by substantially longer tracking, potentially reducing the need for exhaustive simulations in studies of existing accelerators and the design of future facilities. The methods could also complement classical chaos indicators and modern approaches based on artificial intelligence.

In this work, we develop a generalized Hamiltonian formalism for spatially nonlocal field theories whose Lagrangian densities depend explicitly on both the local field and its spatially reflected counterpart. Starting from a generalized variational principle, we derive generalized Euler-Lagrange equations and introduce a generalized functional derivative that consistently accounts for reflected-field contributions. The proposed formalism is applied to three spatially nonlocal nonlinear Schrödinger equations. For the Ablowitz-Musslimani equation, we construct, to the best of our knowledge, the first standard Lagrangian density formulated directly in terms of the complex fields and derive its complete Hamiltonian formulation. The same framework is subsequently applied to the two Lagrangian nonlocal nonlinear Schrödinger equations introduced by Velasco-Juan and Fujioka, yielding consistent Hamiltonian formulations that reproduce the corresponding generalized Euler-Lagrange equations. These results establish a unified Hamiltonian framework for a broad class of spatially nonlocal nonlinear field theories.

The Lorenz system is the canonical low-order model of convective instability, yet its dissipative and driving terms have never been checked against, nor constructed from, an explicit thermodynamic bookkeeping. We derive a modification that satisfies the local-detailed-balance condition for macroscopic relaxation toward nonequilibrium steady states, thereby identifying the thermodynamic force, entropy-production rate and frenesy of the resulting flow. The resulting model produces a transition from a quiescent fixed point to a robust, large-amplitude relaxation oscillation, closely analogous to recharge-discharge oscillator paradigms used for the El Nino-Southern Oscillation. The system alternates between a long, nearly reversible recharge phase and a brief, violently frenetic discharge burst, during which essentially all of the cycle's activity and entropy production is concentrated.

Optical frequency combs (OFCs) are revolutionary light sources characterized by discrete and equally spaced spectral lines, and they have found widespread applications in metrology, spectroscopy, and communications. In the early stages, OFCs were realized using mode-locked lasers. With advancements in the fabrication of high-quality factor ($Q$) microresonators and the increasing demand for miniaturized and integrable photonic chips, microresonator-based OFCs, commonly referred to as microcombs, have been developed. Although early studies of microcombs primarily focused on single microresonators or a few resonators, a significant breakthrough occurred in 2021 when it was theoretically predicted that light propagating in the topological edge channel of an array of ring resonators could generate nested frequency combs known as topological OFCs. Since then, the field of topological OFCs has progressed rapidly, with experimental observations made in 2024. This Perspective will introduce the history of OFCs, placing particular emphasis on the emergence and development of topological OFCs, as well as exploring the research challenges and opportunities associated with them.

Dendritic silhouettes are cheap. Pond ice neurons, diffusion-limited aggregates, wiring-minimizing arbors, and biological neurons can share a branching statistic because they share a growth instability, not because they are functionally similar. We take that observation as a modeling constraint. The ice neuron only structurally resembles a biological neuron. A biological neuron carries an internal state and a field that it helps generate, letting geometry, dynamics, and field co-sculpt one another. In our construction, the ice neuron is the zero-gain limit of an interfacial system in which a FitzHugh–Nagumo state lives on the moving front and writes back into the field that drives it, with dimensionless gain $\gamma$. At $\gamma=0$, morphology and internal dynamics decouple into a Mullins–Sekerka growth band and a damped oscillator. Finite $\gamma$ deforms the spectrum, shifts the selected instability, and opens an oscillatory region below the isolated Hopf threshold. On a dendrite grown at $\gamma=0$ and then frozen, field closure increases distal persistence without changing the geometry; its first-order effect remains positive across repeated realizations of the same growth rule and localizes to a small set of non-cable field shortcuts. Nonlinear dynamics remain finite beyond the linear instability: below threshold, closure is almost the transmission; at spike amplitude, it is a correction to a regenerative cable event. A dendritic tree therefore is not intrinsically a structural spandrel. Shape alone is morphology without an internal state or a source term. That is the dendritic silhouette of the ice neuron. Excitability adds an active cable. Field closure is the biophysical mélange in which the neuron helps write the field that acts back on it, sculpting its dynamics and, during growth, its structure. That is a biological neuron.

In this paper, we explore a variable-coefficient Fokas-Lenells system with gain/loss. Applying a nonlocal symmetry reduction, the system leads to a PT-invariant variable-coefficient Fokas-Lenells equation (vcFLE) with specific constraints on dispersion and nonlinearity coefficients. We prove the Lax integrability of the given system and construct one- and two-soliton solutions using a nonstandard Hirota bilinearization method. We show that, depending on the soliton parameters, the solutions can retain or violate the underlying symmetry. With appropriate choice of dispersion and nonlinearity, the system shows periodic-wave and soliton structures.

We study the numerical evaluation of contour integral representations for solutions of the linearized second and fourth Painlevé equations. The construction combines a nonlinear background solution, continuation of canonical Lax-pair columns, normalization matching, and quadrature over cycles with decaying branches in distinct sectors. Accuracy is assessed using a reference fundamental matrix, differential-equation residuals, and variations of monodromy data computed independently from the spectral problem. An exact identity describes the propagation of initial-basis errors. A stagewise a posteriori refinement criterion uses increments in monodromy variations to allocate numerical accuracy. The contour algorithm is compared with the Dormand–Prince method on a regular background, in a region of rapid solution growth, and through hard loss of stability followed by approximately regular oscillations. Accumulated offsets in monodromy variations and subsequent drift are measured separately. The experiments show how drift depends on spectral conditioning, quadrature accuracy, and error allocation between initial and subsequent parts of the computation. Five additional initial-data and parameter cases demonstrate the dependence of comparative performance on the nonlinear background.

Multiple pendulum with branching is proposed as a convenient platform to study energy transfer between different modes. The antiphase oscillation mode is localized to the "child" links, which makes it easy to prepare initial conditions without exciting the antiphase mode. A manageable expression for the energy transfer is derived theoretically and evaluated with experimental data.

In 2020 M. Mucciconi and L. Petrov introduced a long-range deformation of the quantum open non-relativistic Toda system. We prove the integrability of the deformed Toda system by constructing a $2 \times 2$ Lax operator, which produces the commutative family of differential operators containing the Hamiltonian of the deformed Toda system. Moreover, we show that the same integrable deformation exists on both classical and quantum levels and can be applied to both non-relativistic and relativistic Toda systems. For the open non-relativistic deformed Toda systems we also present an $n \times n$ Lax matrix and prove that it produces the same family of Hamiltonians. We also show how to obtain the van Diejen-type deformed Toda system. Lastly, we show that on the quantum level the algebraic Bethe ansatz technique can be applied to the deformed Toda system.

Phase reduction is a powerful technique to obtain phase models from highly dimensional oscillatory systems. Starting from pairwise interactions, a first order approximation yields Kuramoto and Winfree phase models, while nonpairwise interactions emerge at the second order. Recently, phase models have been extended to account for nonpairwise interactions, notably, the higher-order Kuramoto model. Is there an intrinsic difference between physical and emergent nonpairwise interactions? And can we make use of the former to even out the effects of the latter? In this work, we exploit a recently developed parametrization method to compute the phase reduction and answer these questions. After a comparison of the new method with the classic Kuramoto-style phase reduction, we solve the network motifs of the emergent nonpairwise interactions and compare them with the physical once. Lastly, we exploit the motifs to adopt a coupling design approach. Our framework paves the way for further exploitations of physical nonpairwise interactions for applications in synchronization engineering.

The great transoceanic voyages of the Age of Discovery imposed severe human costs on their crews. Long travel durations, deteriorating sanitary conditions, and finite food and water supplies made disease, exhaustion, and mortality central features of maritime expansion. We propose a minimal compartmental model describing the voyage itself as a dynamical system coupling crew health and resource depletion. The population is divided into healthy, ill, and dead individuals, while an additional variable represents the effective level of provisions on board. Illness incidence depends nonlinearly on resource availability, capturing the increasing health risk associated with deteriorating food, water, and sanitary conditions. The coupled ordinary differential equations are analyzed in the transient regime relevant to finite-duration voyages. Rather than reconstructing specific historical expeditions, the model identifies generic mechanisms linking provisioning, health deterioration, mortality, and voyage duration. Numerical results show that resource limitations can substantially reduce the healthy and operational fraction of the crew even when overall survival remains relatively high, revealing a distinction between demographic survival and operational viability. Increasing initial provisioning delays the crossover at which ill individuals become as numerous as healthy individuals and extends the period over which the crew remains predominantly operational. We also introduce a critical provisioning time, defined as the time required for resources to reach a prescribed stress threshold, and find an approximately linear dependence on the initial resource buffer. The framework provides a quantitative perspective on the internal dynamics of early transoceanic navigation and illustrates how dynamical-systems approaches can contribute to the study of historical processes.

We study the off-diagonal matrix elements of real-space observables in time-reversal-invariant single-particle chaotic systems. By analyzing the semiclassical expression for the off-diagonal variance derived from Berry's conjecture, we show that the banded structure of the observable matrix emerges naturally. For local observables, we identify a characteristic bandwidth associated with a late-time timescale inversely proportional to the particle velocity. We further show that, for systems with steep-wall confinement, the predicted magnitude follows the entropy scaling of the eigenstate thermalization hypothesis (ETH), multiplied by an additional kinetic-energy-dependent factor that is independent of spatial dimension and is not captured by conventional many-body ETH. We illustrate these results through a case study of quantum billiards and verify the semiclassical scaling numerically in a generalized quarter-Sinai billiard. Our results elucidate the dynamical implications of Berry's conjecture and provide a comparison between single-particle eigenstate thermalization and many-body ETH.

In this article, we study the dipolar Bosonic quantum fluid. The fluid experiences mean-field, beyond mean-field, and three body interactions. We investigate their competition with dipolar interaction and fragmentation as a result of this competition. We further investigate the elementary excitations and note two distinct dispersion regimes, namely roton-mode and modulational instability. We support our observation by calculating the superfluid fraction and the condensate fraction.

We introduce a trigonometric version of the Nosé-Hoover oscillator in which the quadratic mechanical terms and unbounded thermostat coupling are replaced by bounded trigonometric functions. This formulation replaces the harmonic potential by a pendulum-type potential and confines the thermostat interaction to a bounded periodic form. The resulting two-parameter system is naturally defined on the three-dimensional torus and reduces near the origin, to leading order, to the classical polynomial Nosé-Hoover model. We investigate its global dynamics using Poincaré sections, bifurcation diagrams, Lyapunov spectra, Kaplan-Yorke dimensions, and the Lyapunov Integrability Test (LIT). The numerical results reveal the coexistence of regular and chaotic dynamics and characterize changes in dissipative behavior across the parameter plane. We then analyze the two limiting cases associated with the parameter axes. For $a=0$, we construct two functionally independent first integrals on regular domains, whereas for $b=0$ the dynamics reduces to a family of two-dimensional systems on invariant tori, which are analyzed using Darboux polynomials and exponential factors. First-order averaging near the intersection of these integrable limits yields periodic solutions bifurcating from unperturbed periodic orbits and an obstruction to regular $C^1$ first integrals in their neighborhoods. Independently, differential Galois theory applied to the normal variational equation, together with the Ayoul-Zung and Li-Shi criteria, excludes meromorphic $B$-integrability and non-constant meromorphic first integrals near a particular non-equilibrium phase curve for $ab\neq0$. Thus, despite retaining the local structure of the classical Nosé-Hoover oscillator, its trigonometric counterpart exhibits markedly different global dynamics and integrability.