2026-04-15 | | Total: 20
Neutrino billiards serve as a model system for the study of aspects of relativistic quantum chaos. These are relativistic quantum billiards consisting of a spin-1/2 particle which is confined to a planar domain by imposing boundary conditions on the spinor components which were proposed in [Berry and Mondragon 1987, {\it Proc. R. Soc.} A {\bf 412} 53) . We review their general features and the properties of neutrino billiards with shapes of billiards with integrable dynamics. Furthermore, we review the features of two neutrino billiards with the shapes of billiards generating a chaotic dynamics, whose nonrelativistic counterpart exhibits particular properties. Finally we briefly discuss possible experimental realizations of relativistic quantium billiards based on graphene billiards, that is, finite size sheets of graphene.
In generic Hamiltonian systems that are neither completely integrable nor fully chaotic, phase space consists of a mixture of regular and chaotic components. In classical dynamics, transitions between different invariant sets in phase space are strictly forbidden, and these sets act as dynamical barriers to one another. In quantum mechanics, in contrast, wave effects allow transitions through such dynamical barriers. This process, known as dynamical tunneling, refers to penetration through dynamical barriers in phase space and was first recognized in the early 1980s. Since then, various aspects of dynamical tunneling have been elucidated, significantly advancing our understanding of such a novel quantum phenomenon. In this article, we provide an overview of several phenomenological perspectives of dynamical tunneling, including chaos-assisted and resonance-assisted tunneling, and also introduce approaches based on classical mechanics extended into the complex domain. In particular, we seek to clarify what is meant by the common claim that "chaos leads to an enhancement of the tunneling probability", which is often made when dynamical tunneling is dressed. We discuss what regime this refers to and, if such an enhancement occurs, what its likely origin is.
We propose a theory based on dynamical systems to explain and predict the occurrence of extreme events, of which critical transitions form a subset. In fast-slow nonlinear systems, we identify a cascade of events preceding extreme events: (i) a slow regime, in which the fast covariant Lyapunov vectors (CLVs) are both tangent to the fast eigenvectors and remain transversal to the slow subspace; (ii) a transition regime, in which the fast eigenvalues become neutrally stable while the fast CLVs are no longer tangent to the fast eigenvectors; and (iii) a critical regime, in which a strong spectral gap in the eigenvalues causes both fast and slow CLVs to become tangent along the dominant fast direction, breaking the transversality between fast and slow subspaces. Building on this cascade, we propose two precursors to forewarn the occurrence of extreme events. We numerically test the theory and precursors on low- and higher-dimensional systems. The proposed precursors predict extreme events and critical transitions with 100% precision and recall. This work opens opportunities for time-forecasting extreme events using theoretically grounded precursors.
We study large-scale dynamics in the Ginzburg-Landau equation (GLE) using a reduced description derived from a WKB expansion. Rigorous mathematical results establishing that this reduced equation accurately approximates the full GLE are currently limited to the real GLE (RGLE) and exclude phase-slip dynamics. For the RGLE, we find that the reduced equation has conserved gradient form and show that, upon inclusion of a higher-order regularization, it admits exact stationary solutions. In the reduced dynamics, all nonuniform steady states are linearly unstable and among them, localized hole solutions identified through the reduced description differ from the classical hole solution of the RGLE due to Langer and Ambegaokar. In the Eckhaus-unstable regime, we derive a self-similar description of the approach to finite-time singularities in the reduced equation, with scaling exponents that agree with direct numerical simulations (DNS), and a similarity profile obtained from a nonlinear 4th-order boundary value problem. Extending the reduction to the complex GLE (CGLE) with nearly real coefficients introduces a Burgers nonlinearity that generates traveling shocks connecting two distinct plane-waves. We obtain exact expressions for the shock profile and perform extensive DNS to demonstrate convergence to the predicted profile in the appropriate large-scale, nearly real-coefficient limit of the CGLE. Away from this limit, the wave number profile loses monotonicity, which we explain in the framework of spatial dynamics. We further show that the exact shock solutions found here are qualitatively distinct from the Nozaki-Bekki solutions. Taken together, our results reveal how a single, scalar reduced equation elucidates unstable stationary states, self-similar collapse toward phase slips, and shock formation, providing an understanding large-scale phase dynamics in pattern-forming systems.
Spatiotemporal chaotic systems are difficult to characterize in a model-free manner because of their high dimensionality, strong nonlinearity, and sensitivity to initial conditions. Coupled map lattices, as a representative class of extended nonlinear systems, exhibit diverse regimes such as frozen random pattern, defect chaotic diffusion, and fully developed turbulence. In this work, we propose an ensemble version of multiplexing local reservoir computing for the data-driven characterization of spatiotemporal chaos. By constructing multiple base learners with randomized hyperparameters and combining their outputs, the method improves prediction robustness and quantifies predictive uncertainty through ensemble spread. More importantly, we show that this uncertainty contains direct dynamical information. It identifies frozen positions in frozen random pattern, supports the estimation of defect diffusion coefficients in defect chaotic diffusion, and provides an effective indicator of chaotic intensity in fully developed turbulence. Analyses of the spatial power spectrum and Lyapunov exponent spectrum further support the consistency between the uncertainty field and the intrinsic dynamical properties of the system. These results show that ensemble reservoir computing can serve not only as a prediction tool but also as a data-driven framework for the dynamical characterization of high-dimensional nonlinear systems.
We propose a framework to infer the coupling strength and the natural frequency distribution in a coupled Stuart-Landau oscillator system with a large population. The inference method uses observation of linear response of a macroscopic quantity and of an oscillator. We first solve the direct problem on the response with transforming the system into the phase-amplitude equations. Solving the inverse problem, we show that the coupling strength is inferred from observation of an oscillator and the natural frequency distribution from macroscopic responses. The proposed method requires only one-dimensional observation in the two-dimensional Stuart-Landau system. Validity of the inference theory is examined by numerical simulations.
This chapter gives an overview of transport problems where chaotic dynamics of the system plays a crucial role. We begin with single-particle transport problems and then come to conservative and then dissipative systems of identical particles, which follows the historical way of developing the theory of Quantum Chaos over the past 40 years. We also include brief descriptions of key laboratory experiments on the discussed transport problems.
Through semiclassical methods the subject of quantum chaos motivates and depends on Hamiltonian chaos research. Presented here is a selection of Hamiltonian chaos topics that in this way get directly related to any of a variety of quantum chaos research problems. The chapter begins with a description of various useful theoretical and computational tools of chaos research, e.g.~surfaces of section, paradigms of chaos, stability analysis, and symbolic dynamics... This is followed by discussions regarding the geometry of chaos, how chaotic systems respond to perturbations, and the complexification of Hamiltonian dynamics. The emphasis is on intuitive explanations and illustrations of various ideas with the references containing more mathematically rigorous expositions.
Organelle patterning and its heritability remain central mysteries in cell biology, highlighting the fundamental tension between genetic inheritance and self-assembly. Here, we explore the nonequilibrium assembly and size control of the Golgi complex and endosomes, amid a continuous flux of membrane traffic, within a stochastic framework of mechanochemical fusion-fission cycles that violate detailed balance. Using a dynamical systems approach, we identify distinct, robust regimes, ranging from fixed points to limit cycles with definite phase relations. We identify these dynamical regimes with diverse phenotypes, from stable cisternae to periodic, cell-cycle-dependent dissolution/reassembly to cisternal progression. We analyse its dynamic response to systematic perturbations or driving protocols and make definite predictions that may be tested experimentally. Our analysis reveals that the two competing models of Golgi organization-vesicular transport and cisternal progression - are, in fact, two phases of the same underlying nonequilibrium process. Finally, our framework offers a strategy for controlling cisternal chemical identity and number and by modulating the interplay between glycosylation enzymes and membrane fission-fusion dynamics.
We present an integral formulation of classical Yang-Mills theory coupled to fermionic and scalar matter fields in (1+1)-dimensional Minkowski spacetime. By reformulating the local dynamics in terms of loop-space holonomies, we demonstrate that the path independence of the holonomy eigenvalues constitutes a conservation law, yielding an infinite hierarchy of gauge-invariant, dynamically conserved charges. While a zero-curvature equation is associated with a necessary condition for this path invariance, we note that it is not strictly sufficient on its own. Employing a first-order symplectic formalism, we show that these non-abelian charges generate global symmetry transformations on the fundamental phase-space variables. We rigorously prove that these transformations preserve the physical dynamics, leaving the total Hamiltonian invariant up to first-class constraints. Furthermore, an analysis of the Poisson algebra reveals that these conserved charges are in involution, provided the boundary integration constant lies within the center of the gauge group. This exact, lower-dimensional framework provides a highly tractable setting to investigate the algebraic structures of these hidden symmetries and the meaning of the conserved charges as physical observables, establishing a classical foundation for exploring their role in the quantum regime, such as in strongly coupled lattice gauge theories.
In quantum systems with a classical limit, advanced semiclassical methods provide the crucial link between phase-space structures, reflecting the distinction between chaotic, mixed or integrable classical dynamics, and the corresponding quantum properties. Well established techniques dealing with ergodic wave interference in the usual semiclassical limit $\hbar \to 0$, where the classical limit is given by Hamiltonian mechanics of particles, constitute a now standard part of the toolkit of theoretical physics. During the last years, these ideas have been extended into the field theoretical domain of systems composed of $N$ indistinguishable particles, aka quantum fields, displaying a different type of semiclassical limit $\hbar_{\rm eff}=1/N \to 0$ and accounting for genuine many-body quantum interference. The foundational concept behind this idea of many-body interference, the many-body version of the van Vleck-Gutzwillers semiclassical propagator, is explained in detail. Based on this the corresponding semiclassical many-body theory is reviewed. It provides a unified framework for understanding a variety of quantum chaotic phenomena addressed, including random-matrix spectral correlations in many-body systems, the universal morphology of many-body eigenstates, interference effects kin to mesoscopic weak localization, and the key to the scrambling of many-body correlations characterized by out-of-time-order correlators.
Mesoscopic devices, with system sizes in the range of several to several dozens wavelengths, represent paradigmatic model systems for the observation of quantum chaotic behaviour based on semiclassical concepts. Those electronic and photonic billiard cavities are small enough for interference effects not to be ignored. Nonetheless, the classical ray or particle tracing picture can often provide a substantial understanding of the dynamics of the system along the lines of classical-quantum, or ray-wave correspondence. This well-established principle turns out to be particularly useful when applied not only in real space, but by extending it to phase space such that both location and momentum information can contribute to a deeper and more comprehensive understanding of the dynamical behaviour.
One of the fundamental manifestations of classical chaos is exponential sensitivity to initial conditions that is, two trajectories starting from nearly identical initial states diverge exponentially over time. This behavior is quantified by the Lyapunov exponents. Due to the unitary nature of quantum mechanics, such exponential divergence is elusive in quantum systems. As a result, several alternative quantities have been proposed and studied in recent years to capture analogous behavior. In this article, we present a pedagogical overview of three such quantities that have been the focus of intense research in recent years: the Loschmidt echo, out-of-time-order correlators (OTOCs), and Krylov complexity.
Out-of-Time-Order Correlators (OTOCs) quantify quantum information scrambling, but their connection to localized phase-space structures, such as chemical transition states, requires formal development. We derive a leading-order semiclassical expansion for the local microcanonical OTOC in systems with an index-1 saddle point, expressing the scrambling rate as a coherent sum over unstable periodic orbits on the Normally Hyperbolic Invariant Manifold (NHIM). Valid in the semiclassical limit and the intermediate-time regime before the Ehrenfest time, our derivation utilizes the Normal Form theory of the transition state, which transforms the Hamiltonian near the saddle into an integrable (though generally non-separable) form dependent on conserved actions. We outline the derivation of the microcanonical trace, the semiclassical propagator for integrable systems, the factorization of the stability matrix, and the Schur complement reduction of the stationary phase approximation. Our result extends periodic-orbit trace methods to scrambling observables, yielding a local instability exponent Λ(J) governing the leading semiclassical growth window. As a special case, when the observation time coincides with the intrinsic periods of the contributing orbits, the trace sum reduces to an effective 1.5Λ scaling, resulting from the competition between local hyperbolic growth and wavepacket dilution. This simplified form is conditional; the full expansion retains a coherent sum over orbit periods. Finally, we discuss how the dependence of the instability on transverse actions establishes a theoretical mechanism for mode-selective control of scrambling, and outline a numerical evaluation strategy to test these predictions.
The kicked rotor provides a simple yet powerful model for introducing many of the central concepts of classical and quantum chaos. Despite its apparent simplicity, it exhibits rich dynamical behavior and has found applications across a wide range of fields, including atomic and optical physics, condensed matter physics, and emerging quantum technologies. This chapter begins by exploring foundational ideas using the kicked rotor as a unifying framework. We first discuss the transition from regular to chaotic motion in the classical system, and then introduce key quantum phenomena such as dynamical localization and quantum resonances. Special attention is devoted to the emergence of characteristic time scales and their role in the quantum-classical correspondence. To make these ideas more concrete, we also provide a brief overview of experimental realizations of the kicked rotor and its variants, illustrating how theoretical concepts are implemented in practice. In the second part of the chapter, we guide the reader toward more recent and advanced developments. Topics include near-resonant dynamics, topological features of kicked systems, the emergence of quantum dynamical phases inferred from classical transport properties, and extensions to non-Hermitian physics. We conclude with a discussion of open problems and future perspectives, outlining directions in which the kicked rotor continues to offer valuable insights.
The notion of Shannon entropy is crucial for the theory of classical information. In quantum information theory, an analogous key role is played by the von Neumann entropy: quantum information processing is closely related to entropy dynamics. This reveals a direct link with the theory of quantum chaotic systems, which can be characterized by a positive entropy production. Furthermore, noise, which inevitably affects any quantum system, can be modeled by a random quantum operation or by coupling to an environment in a generic chaotic state. In this contribution, we emphasize the universality of quantum chaotic dynamics and discuss its implications for quantum information processing.
The Euler ensemble was recently supported by large-scale ($4096^3$) direct numerical simulations as the universal statistical attractor of decaying fluid turbulence. Previous mathematical derivations of this ensemble relied on measure-theoretic limits of discrete polygonal loop equations. In this Letter, we present a continuous algebraic derivation. By reformulating the Navier-Stokes equation as a covariant derivative operator flow in the Lagrangian frame, we analytically eliminate advection. Applying Feynman's operational calculus, the 3D non-commutative operator algebra maps to ordering discontinuities (finite-difference jumps) on a 1D momentum loop. This continuous formalism reduces to the discrete, number-theoretic geometric quantization of the Euler ensemble via roots of unity without requiring spatial lattice approximations, demonstrating that macroscopic fluid chaos is a deterministic projection of the Farey sequence.
Factor models characterize the joint behavior of large sets of financial assets through a smaller number of underlying drivers. We develop a network-based framework in which factors emerge naturally from the structure of interactions among assets rather than being imposed statistically. The market is modeled as a system of coupled iterated maps, where assets' return depends on its own past returns and those of related assets. Effectively modeling the influence of irrational traders whose decisions are based on the past movements of a collection of stocks. The interaction structure between stock returns is defined by a coupling matrix derived from an orthogonal transformation of a Laplacian matrix that gradually links initially isolated clusters into a fully connected network. Within this structure, stable patterns of co-movement arise and can be interpreted as financial factors. The relationship between the initial clustering and the number of observed factors is consistent with a center manifold reduction. We identify an optimal regime in which assets' variance is effectively explained by the set of factors produced by the network. Our framework offers a structural perspective based on interaction-based factor formation and dimension reduction in financial markets.
The 100,000-year problem concerns the dominant period of glacial-interglacial cycles over the past 800,000 years and their correlation with Earth's orbital eccentricity, despite eccentricity's weak influence on solar radiation. Two theories compete: the astronomical theory, in which orbital forcing drives the cycles with amplification from Earth system feedbacks, and the geochemical theory, in which internal dynamics dominate with orbital forcing synchronising oscillations. We investigate these theories using conceptual models. Augmentations to the Budyko energy balance model fail to reproduce the 100,000-year period, revealing formulation limitations. Linearised versions of existing non-linear ice volume models perform comparably to their full counterparts, indicating the data does not necessitate non-linear dynamics. We develop two simple linear models: a feedforward model aligned with the astronomical theory and a feedback model aligned with the geochemical theory. The feedforward model reproduces the ice volume record well and offers a novel explanation for the absence of eccentricity's 400,000-year period, arising from oceanic heat storage and tropospheric energy responding with differing phase lags. Conservative estimates show bulk ocean temperature variation can be explained by eccentricity alone, challenging the geochemical theory's core assumption. We also show that widespread use of Q65 may bias models towards geochemical explanations by underrepresenting eccentricity. The feedback model's improvement is concentrated around Marine Isotope Stage 11, suggesting this anomalous interglacial reflects Earth-based events rather than a general requirement for feedback mechanisms. We conclude that 800,000 years of glacial cycles can be largely reproduced by a linear astronomical model, emphasising the importance of parsimony when interpreting palaeoclimate data.
Self-gravitating systems of $N$ particles are chaotic. We wonder how chaotic the Galaxy is, and what the consequences are. We therefore simulate the dynamical evolution of a galaxy-scale distribution of point masses in order to measure the degree of chaos in such a system. These calculations were performed using the softened gravitational $N$-body tree-code Bonsai, with up to 40 million equal-mass particles. Smaller simulations were performed to establish the scaling of the Lyapunov time $t_L$ with $N$. We establish the relations between the degree of chaos, the number of particles, and the softening length in the gravitational force calculation of large-scale $N$-body simulations. The moment the bar forms appears insensitive to infinitesimal perturbations to the initial realisation. In contrast, the bar strength and its further evolution sensitively depend on such perturbations. Interestingly enough, the run-to-run variation in the bar strength has its maximum around the maximum bar strength, and drops to the moment the bar buckles. The galaxies we simulated are highly chaotic, but the softening in the simulations suppresses chaos. Still, our models show considerable variations in the macroscopic behaviour due to infinitesimal perturbations to the initial conditions. Real galaxies, however, should be orders of magnitude more chaotic than our simulations, and we are unable to quantify their consequences. Smooth galactic potentials to study individual stellar orbits should be handled with caution on timescales longer than the Lyapunov time. Extrapolating to the number of stars in the Galaxy, ignoring planets and other minor bodies, we conclude that the Milky Way-size galaxies are chaotic on a timescale $\lesssim 0.1$ Myr.