CyberSec Research

Fix a smooth Morse function $U\colon \mathbb{R}^{d}\to\mathbb{R}$ with finitely many critical points, and consider the solution of the stochastic differential equation \begin{equation*} d\bm{x}_{\epsilon}(t)=-\nabla U(\bm{x}_{\epsilon}(t))\,dt \,+\,\sqrt{2\epsilon}\, d\bm{w}_{t}\,, \end{equation*} where $(\bm{w}_{t})_{t\ge0}$ represents a $d$-dimensional Brownian motion, and $\epsilon>0$ a small parameter. Denote by $\mathcal{P}(\mathbb{R}^{d})$ the space of probability measures on $\bb R^d$, and by $\mathcal{I}_{\epsilon} \colon \mathcal{P}(\mathbb{R}^{d})\to[0,\,\infty]$ the Donsker--Varadhan level two large deviations rate functional. We express $\mc I_\epsilon$ as $\mc I_\epsilon = \epsilon^{-1} \mc J^{(-1)} + \mc J^{(0)} + \sum_{1\le p\le \mf q} (1/\theta^{(p)}_\epsilon) \, \mc J^{(p)}$, where $\mc J^{(p)}\colon \mc P(\bb R^d) \to [0,+\infty]$ stand for rate functionals independent of $\epsilon$ and $\theta^{(p)}_\epsilon$ for sequences such that $\theta^{(1)}_\epsilon \to\infty$, $\theta^{(p)}_\epsilon / \theta^{(p+1)}_\epsilon \to 0$ for $1\le p< \mf q$. The speeds $\theta^{(p)}_\epsilon$ correspond to the time-scales at which the diffusion $\bm{x}_{\epsilon}(\cdot)$ exhibits a metastable behaviour, while the functional $\mc J^{(p)}$ represent the level two, large deviations rate functionals of the finite-state, continuous-time Markov chains which describe the evolution of the diffusion $\bm{x}_{\epsilon}(\cdot)$ among the wells in the time-scale $\theta^{(p)}_\epsilon$.

Electro-viscoelastic theory for polymer melts has been extensively studied experimentally for the past century, primarily for manufacturing purposes. However, the modeling and theory for this have been minimal, leaving many questions on the mechanisms and behavior of an arbitrary flow scheme. To remedy this, previously solved overdamped Langevin equations for the Doi-Rouse model are modified to include charge and electric field potential forces. The charge sequence on the chain is hypothesized to be a cosine sequence along the chain, resembling multiple electric dipoles that conveniently correspond to a Rouse mode of the chain. These are then solved for the shear stress under homogeneous shear rates and electric fields to find directional viscosity increases depending on the shearing and electric field orientation. Using the newly derived shear stress from the Doi-Rouse approach, a continuum model is proposed that resembles a modified upper-convected Maxwell model, including polarization stresses in terms of an electric field dyadic. This new continuum model, named the upper-convected electro-Maxwell model, is verified using Kremer-Grest polymer chains simulated with molecular dynamics for multiple flow schemes and a specified charge sequence along the chain. Furthermore, the MD results verified the difference in the overall and charge sequence relaxation times through the shear and normal stress polarizations, showing the necessity for the upper-convected derivative of the electric field dyadic to correct the viscosity scaling. Finally, the dynamic properties of the polarized polymer melt are examined analytically, finding that the phase shift is unaffected by the electric field contribution.

We investigate self-organized criticality in a two-dimensional electron gas (2DEG) by introducing a lattice-based model that incorporates electron-electron interactions through the concept of coherence length. Our numerical simulations demonstrate that in the strongly interacting regime, the system exhibits a distinct set of universal critical exponents, markedly different from those observed in the weakly interacting limit. This dichotomy aligns with experimental findings on the metal-insulator transition in 2DEGs, where high interaction strength (low carrier density) leads to qualitatively different behavior. The analysis includes scaling of the average electron density with temperature, the power spectral density, and the statistics of electronic avalanches - namely their size distributions and autocorrelation functions. In all cases, the extracted exponents differ significantly between the weak and strong interaction regimes, highlighting the emergence of two universality classes governed by interaction strength. These results underscore the critical role of electron correlations in the self-organized behavior of low-dimensional electronic systems.

In this work, we investigate the impact of conserved charges on the dynamics of spread complexity of quantum states. Building on the notion of symmetry-resolved Krylov complexity [1], we extend the framework to general quantum states and analyze the relation between the total spread complexity and its decomposition into fixed-charge sectors. After exploring a range of analytical examples and using orthogonal polynomial approach, we identify conditions under which spread complexity exhibits equipartition across sectors. Finally, we discuss quantum speed limits that constrain the growth of complexity in the presence of conserved charges.

For the system with inhomogeneous distribution of macroscopic parameters we obtain thermodynamic relation which depends on the spatial point (coordinate). In our approach, to obtain such a relation we use the basic ideas of the method of nonequilibrium statistical operator combined with the Hubbard-Stratonovich transformation. First of all, we define the thermodynamic relation for the system with homogeneous distribution of particles. Possible behavior peculiarities of systems with different character of interaction in nonequilibrium case are predicted. By saddle-point method we find the dominant contributions to the partition function and obtain all thermodynamic parameters of the systems with different character of interaction. The formations of saddle state in all systems of interacting particles at different temperatures and particle distributions have the same physical nature and therefore they can be described in the same way. We consider the systems with attractive and repulsive interactions as well as self-gravitating systems.

An interesting class of time dependent backgrounds in $1+1$ dimensional string theory involves worldsheet Liouville walls which move in (target space) time. When a parameter in such a background exceeds a certain critical value, the speed of the Liouville wall exceeds the speed of light, and there is no usual S-Matrix. We examine such backgrounds in the dual $c=1$ matrix model from the point of view of fluctuations of the collective field, and determine the nature of the emergent space-time perceived by these fluctuations. We show that so long as the corresponding Liouville wall remains time-like, the emergent space time is conformal to full Minkowski space with a time-like wall. However, for the cases where the Liouville wall is superluminal, the emergent space-time has a {\em space-like boundary} where the collective field couplings diverge. This appears as a space-like singularity in perturbative collective field theory. We comment on the necessity of incorporating finite $N$, as well as finite (double-scaled) coupling, effects to understand the behavior of the exact theory near this boundary.

We present a quantum algorithm for simulating open quantum systems coupled to Gaussian environments valid for any configuration and coupling strength. The algorithm is for instance applicable to problems with strongly coupled, or non-Markovian, environments; problems with multiple environments out of mutual equilibrium; and problems with time-dependent Hamiltonians. We show that the algorithm can reproduce the true dynamics of such problems at arbitrary accuracy and, for a broad range of problems, only adds a minor resource cost relative to Trotterized time evolution; the cost is polynomial in the inverse target accuracy. The algorithm is based on the insight that any Gaussian environment can be represented as a train of ancillary qubits that sequentially interact with the system through a time-local coupling, given by the convolution square root of the bath correlation function; this is a secondary result of our work. Our results open up new applications of quantum computers for efficient simulation of non-equilibrium and open quantum systems.

We apply Nambu non-equilibrium thermodynamics (NNET)-a dynamics with multiple Hamiltonians coupled to entropy-induced dissipation-to paradigmatic far-from-equilibrium systems. Concretely, we construct NNET realizations for the Belousov-Zhabotinsky (BZ) reaction (oscillatory), the Hindmarsh-Rose neuron model (spiking), and the Lorenz and Chen systems (chaotic), and analyze their dynamical and thermodynamic signatures. Across all cases the velocity field cleanly decomposes into a reversible Nambu part and an irreversible entropygradient part, anchored by a model-independent quasi-conserved quantity. This construction reproduces cycles, spikes, and strange-attractor behavior and clarifies transitions among steady, periodic, and chaotic regimes via cross-model diagnostics. These results demonstrate that NNET provides a unified, quantitatively consistent framework for oscillatory, spiking, and chaotic non-equilibrium systems, offering a systematic description beyond the scope of linear-response theories such as Onsager's relations or GENERIC.

Diffusive motion is a fundamental transport mechanism in physical and biological systems, governing dynamics across a wide range of scales -- from molecular transport to animal foraging. In many complex systems, however, diffusion deviates from classical Brownian behaviour, exhibiting striking phenomena such as Brownian yet non-Gaussian diffusion (BYNGD) and anomalous diffusion. BYNGD describes a frequently observed statistical feature characterised by the coexistence of linear mean-square displacement (MSD) and non-Gaussian displacement distributions. Anomalous diffusion, in contrast, involves a nonlinear time dependence of the MSD and often reflects mechanisms such as trapping, viscoelasticity, heterogeneity, or active processes. Both phenomena challenge the conventional framework based on constant diffusivity and Gaussian statistics. This review focuses the theoretical modelling of such behaviour via the Langevin equation with fluctuating diffusivity (LEFD) -- a flexible stochastic framework that captures essential features of diffusion in heterogeneous media. LEFD not only accounts for BYNGD but also naturally encompasses a wide range of anomalous transport phenomena, including subdiffusion, ageing, and weak ergodicity breaking. Ergodicity is discussed in terms of the correspondence between time and ensemble averages, as well as the trajectory-to-trajectory variability of time-averaged observables. The review further highlights the empirical relevance of LEFD and related models in explaining diverse experimental observations and underscores their value to uncovering the physical mechanisms governing transport in complex systems.

The Hopfield model describes a neural network that stores memories using all-to-all-coupled spins. Memory patterns are recalled under equilibrium dynamics. Storing too many patterns breaks the associative recall process because frustration causes an exponential number of spurious patterns to arise as the network becomes a spin glass. Despite this, memory recall in a spin glass can be restored, and even enhanced, under quantum-optical nonequilibrium dynamics because spurious patterns can now serve as reliable memories. We experimentally observe associative memory with high storage capacity in a driven-dissipative spin glass made of atoms and photons. The capacity surpasses the Hopfield limit by up to seven-fold in a sixteen-spin network. Atomic motion boosts capacity by dynamically modifying connectivity akin to short-term synaptic plasticity in neural networks, realizing a precursor to learning in a quantum-optical system.

This article presents a theoretical study of the scaling properties of the kinetic energy spectrum in compressible turbulence. From the fundamental symmetries and linear transformations of the microscopic action, we derive exact relations between the correlation functions and their generators. These relations put strong constraints on the possible scaling relations in the system as a function of scale. One of the main results of this study is that the action can be split between an incompressible part, which is the same as in the usual stochastic Navier-Stokes theory whatever the value of the Mach number, and a longitudinal part, whose behavior is to be compared to the three dimensional Burgers equation, which presents a much richer phase diagram as its usually discussed one dimensionalcounterpart.

Understanding and predicting the glassy dynamics of small organic molecules is critical for applications ranging from pharmaceuticals to energy and food preservation. In this work, we present a theoretical framework that combines molecular dynamics simulations and Elastically Collective Nonlinear Langevin Equation (ECNLE) theory to predict the structural relaxation behavior of small organic glass-formers. By using propanol, glucose, fructose, and trehalose as model systems, we estimate the glass transition temperature (Tg) from stepwise cooling simulations and volume-temperature analysis. These computed Tg values are then inserted into the ECNLE theory to calculate temperature-dependent relaxation times and diffusion coefficients. Numerical results agree well with experimental data in prior works. This approach provides a predictive and experimentally-independent route for characterizing glassy dynamics in molecular materials.

We explore the evolution of a strongly interacting dissipative quantum Ising spin chain that is driven by a slowly varying time-dependent transverse field. This system possesses an extensive number of instantaneous (adiabatic) stationary states which are coupled through non-adiabatic transitions. We analytically calculate the generator of the ensuing slow dynamics and analyze the creation of coherences through non-adiabatic processes. For a certain choice of the transverse field shape, we show that the system solely undergoes transitions among classical basis states after each pulse. The concatenation of many of such pulses leads to an evolution of the spin chain under a many-body dynamics that features kinetic constraints. Our setting not only allows for a quantitative investigation of adiabatic theorems and non-adiabatic corrections in a many-body scenario. It also directly connects to many-body systems in the focus of current research, such as ensembles of interacting Rydberg atoms which are resonantly excited by a slowly varying laser pulse and subject to dephasing noise.

The macroscopic fluctuation theory is a powerful tool to characterise the large scale dynamical properties of diffusive systems, both in- and out-of-equilibrium. It relies on an action formalism in which, at large scales, the dynamics is fully determined by the minimum of the action. Within this formalism, the analysis of the statistical properties of a given observable reduces to solving the Euler-Lagrange equations with the appropriate boundary conditions. One must then compute the action at its minimum to deduce the cumulant generating function of the observable. This typically involves computing multiple integrals of cumbersome expressions. Recently, a simple formula has been conjectured to shortcut this last step, and compute the cumulant generating function of different observables (integrated current or position of a tracer) without the need to compute any integral. In this work, we prove this simple formula, and extend it to more general observables. We then illustrate the efficiency of this approach by applying it to compute the variance of a generalised current in the semi-infinite symmetric exclusion process and the joint properties of two occupation times in any diffusive system. In the case of the integrated current, our formula can be interpreted as a generalisation of Fick's law to obtain all the cumulants of the current beyond the average value.

A key challenge in quantum complexity is how entanglement structure emerges from dynamics, highlighted by advances in simulators and information processing. The Lieb--Robinson bound sets a locality-based speed limit on information propagation, while the Small-Incremental-Entangling (SIE) theorem gives a universal constraint on entanglement growth. Yet, SIE bounds only total entanglement, leaving open the fine entanglement structure. In this work, we introduce Spectral-Entangling Strength, measuring the structural entangling power of an operator, and prove a Spectral SIE theorem: a universal limit for R\'enyi entanglement growth at $\alpha \ge 1/2$, revealing a robust $1/s^2$ tail in the entanglement spectrum. At $\alpha=1/2$ the bound is qualitatively and quantitatively optimal, identifying the universal threshold beyond which growth is unbounded. This exposes the detailed structure of Schmidt coefficients, enabling rigorous truncation-based error control and linking entanglement to computational complexity. Our framework further establishes a generalized entanglement area law under adiabatic paths, extending a central principle of many-body physics to general interactions. Practically, we show that 1D long-range interacting systems admit polynomial bond-dimension approximations for ground, time-evolved, and thermal states. This closes the quasi-polynomial gap and proves such systems are simulable with polynomial complexity comparable to short-range models. By controlling R\'enyi entanglement, we also derive the first rigorous precision-guarantee bound for the time-dependent density-matrix-renormalization-group algorithm. Overall, our results extend SIE and provide a unified framework that reveals the detailed structure of quantum complexity.

We consider the mean first passage time (MFPT) for a diffusive particle in a potential landscape with the extra condition that the particle is reset to its original position with some rate r. We study non-smooth and non-convex potentials and focus on the case where the restart rate depends on the space coordinate. There, we show that it is beneficial to restart at a lower rate once you are closer to your intended target.

We continue the work of Belliard, Pimenta and Slavnov (2024) studying the modified rational six vertex model. We find another formula of the partition function for the inhomogeneous model, in terms of a determinant that mix the modified Izergin one and a Vandermonde one. This expression enables us to compute the partition function in the homogeneous limit for the rectangular lattice, and then to study the thermodynamic limit. It leads to a new result, we obtain the first order of free energy with boundary effects in the thermodynamic limit.

The periodic Temperley-Lieb category consists of connectivity diagrams drawn on a ring with $N$ and $N'$ nodes on the outer and inner boundary, respectively. We consider families of modules, namely sequences of modules $\mathsf{M}(N)$ over the enlarged periodic Temperley-Lieb algebra for varying values of $N$, endowed with an action $\mathsf{M}(N') \to \mathsf{M}(N)$ of the diagrams. Examples of modules that can be organised into families are those arising in the RSOS model and in the XXZ spin-$\frac12$ chain, as well as several others constructed from link states. We construct a fusion product which outputs a family of modules from any pair of families. Its definition is inspired from connectivity diagrams drawn on a disc with two holes. It is thus defined in a way to describe intermediate states in lattice correlation functions. We prove that this fusion product is a bifunctor, and that it is distributive, commutative, and associative.