2025-12-08 | | Total: 11
Boundary-catalytic branching processes describe a broad class of natural phenomena where the population of diffusing particles grows due to their spontaneous binary branching (e.g., division, fission or splitting) on a catalytic boundary located in a complex environment. We investigate the possibility of geometric control of the population growth by compensating the proliferation of particles due to branching events by their absorptions in the bulk or on boundary absorbing regions. We identify an appropriate Steklov spectral problem to obtain the phase diagram of this out-of-equilibrium stochastic process. The principal eigenvalue determines the critical line that separates an exponential growth of the population from its extinction. In other words, we establish a powerful tool for calculating the optimal absorption rate that equilibrates the opposite effects of branching and absorption events and thus results in steady-state behavior of this diffusion-reaction system. Moreover, we show the existence of a critical catalytic rate above which no compensation is possible, so that the population cannot be controlled and keeps growing exponentially. The proposed framework opens promising perspectives for better understanding, modeling and control of various boundary-catalytic branching processes, with applications in physics and life sciences.
We investigate mixed eigenstates in systems with sharply-divided phase space, using different piecewise-linear maps whose regular-chaotic boundaries are formed by marginally unstable periodic orbits (MUPOs) or by quasi-periodic orbits. With the overlap index and the entropy localization length, we classify mixed eigenstates and show that the contribution from dynamical tunneling scales as $\sim \hbar\, \exp(-b/\hbar)$, with $b>0$ associated with the relative size of the regular region. The dominant fraction of states that remain sticky to the boundaries, referred to as sticky eigenstates, scales as $\hbar^{1/2}$ in the MUPO case and oscillates around this algebraic behavior in the quasi-periodic case. This behavior generalizes established predictions for hierarchical states in KAM systems, which scale as $\hbar^{1 - 1/γ}$, with $γ$ set by the corresponding classical stickiness reflected in the algebraic decay of cumulative RTDs $t^{-γ}$. For the piecewise-linear maps studied here, $γ= 2$. These results reveal a clear quantum signature of classical stickiness in non-KAM systems.
A dynamical transition -- seen as a sudden increase in the mean-squared displacement at a characteristic temperature that depends on the observation time -- is widely reported in neutron-scattering experiments and molecular dynamics simulations of hydrated proteins. However, its physical origin remains elusive. We show that fluctuating diffusivity in a Langevin framework leads to an observation-time-induced crossover, where the effective diffusion coefficient exhibits a temperature-dependent transition whose crossover point shifts with observation time. Analytical and numerical analyses reveal the mechanism of this crossover and delineate the conditions under which it emerges. Our findings provide a unified nonequilibrium interpretation for observation-time-induced crossover, and suggest that the protein dynamical transition can be viewed as an instance of this general crossover mechanism.
The influence of an external periodic field in a bistable system gives rise to a very interesting phenomenon of dynamic hysteresis, which has extensive applications in technologies. Environmental noise also has a crucial role in controlling this class of hysteresis. The relaxational delay of the system's response towards the periodic drive manifests in the response function-field hysteresis loops. These hysteresis loops are generally symmetric in the range of the external force and the response function. Symmetry breaking is observed under extreme conditions of the noise strength, amplitude and frequency of the periodic force. This symmetric-asymmetric dynamic transition is interpreted in terms of the relevant order parameter defined for the system under study. In our present work, we importantly identify that it is possible to induce symmetry breaking in the hysteresis loops, and consequently to observe dynamic transitions under moderate conditions as well if appropriate asymmetry is implemented in the design of the underlying potential. We develop a quantitative understanding of the variation of the hysteresis loop area and the order parameter with respect to the degree of asymmetry in the potential. Our study suggests an ingenious way to get the desired outputs in the processes related to dynamic hysteresis by effectively controlling the asymmetry of the intrinsic potential governing the dynamics. These findings not only direct a new and impactful path for technological advances involving hysteresis but also provide the basic understanding of the precise role of the underlying potential in regulating this fundamental process, which manifests frequently in natural and designed systems.
In this paper, we review and compare some frameworks for dissipation in non-equilibrium thermodynamics. We start with a brief overview of classical irreversible thermodynamics and gradient dynamics. Then we discuss several specific frameworks including Rayleigh dissipation potential and the dissipative d'Alembert framework, showing their relations with gradient dynamics. Finally, we discuss frameworks for dissipative evolution generated from Poisson brackets.
The relationship between the vanishing of the heat capacities as $T\to0^+$ and the thermal stability is examined. The heat capacities vanish as fast as or faster than $T$ as $T\to0^+$ for states at the phase space boundary ($T=0$) to sustain the standard thermal stability criterion $U_{ss}>0$. Conversely, weakly vanishing heat capacities, which signify a loss of curvature in $U(S)$ at $T=0$, are the signature of a critical condition precisely at $T=0$, as exemplified in marginal Fermi liquids. Therefore, the vanishing of the specific heat should be viewed not as a new law but as a confirmatory result of the existing framework of thermodynamics.
Dormancy is a widespread adaptive strategy that enables populations to persist in fluctuating environments, yet how its benefits depend on the temporal structure of environmental variability remains unclear. We examine how dormancy interacts with environmental correlation times using a delayed-logistic model in which dormant individuals reactivate after a fixed lag while birth rates fluctuate under temporally correlated stochasticity. Numerical simulations and analytical calculations show that the combination of demographic memory and colored multiplicative noise generates a strongly non-monotonic dependence of fitness on dormancy duration, with three distinct performance regimes. Very short dormancy maximizes linear growth but amplifies fluctuations and extinction risk. Very long dormancy buffers environmental variability, greatly increasing mean extinction times despite slower growth. Strikingly, we find a broad band of intermediate dormancy durations that is maladaptive, simultaneously reducing both growth and persistence due to a mismatch between delay times and environmental autocorrelation. An evolutionary agent-based model confirms bistability between short- and long-dormancy strategies, which avoid intermediate lag times and evolve toward stable extremes. These results show that dormancy duration is not merely a life-history parameter but an adaptive mechanism tuned to environmental timescales, and that intermediate "dangerous middle" strategies can be inherently disfavored. More broadly, this work identifies a generic mechanism by which demographic delays interacting with correlated environmental variability produce a non-monotonic fitness landscape that selects for extreme timing strategies.
We investigate uncertainty growth and chaotic dynamics in statistically steady, stably stratified three-dimensional turbulence. Using direct numerical simulations of the Boussinesq equations, we quantify the divergence of initially infinitesimal perturbations via twin simulations and decorrelator diagnostics. At short times, perturbations exhibit exponential growth, allowing us to define a (largest) Lyapunov exponent. We systematically examine how this exponent depends on stratification strength, quantified by the Brunt--Väisälä frequency and the Froude number, in a parameter regime relevant to oceanic flows. We find that increasing stratification leads to a monotonic reduction of the Lyapunov exponent, indicating suppressed chaoticity. Despite this reduction, uncertainty growth retains the universal temporal sequence observed in homogeneous isotropic turbulence -- initial decay, exponential growth, and saturation. The growth phase is characterized by self-similar decorrelator spectra, but exhibits strong anisotropy: uncertainty spreads much more slowly along the stratification direction than horizontally, with the disparity increasing with stratification strength. An analysis of the decorrelator evolution equation reveals that the suppression of chaos arises primarily from strain-mediated alignment dynamics rather than direct buoyancy coupling. Our results provide a quantitative characterization of predictability and uncertainty growth in stratified turbulence and highlight the utility of decorrelator-based methods for anisotropic geophysical flows.
We analyze a size-structured branching process in which individual cells grow exponentially according to a Feller square-root process and divide under general size-control mechanisms. We obtain exact expressions for the asymptotic population growth rate, the steady-state snapshot distribution of cell sizes, and the fluctuations of the total cell number. Our first result is that the population growth rate is exactly equal to the mean single-cell growth rate, for all noise strengths and for all division and size-regulation schemes that maintain size homeostasis. Thus square-root growth noise is neutral with respect to long-term fitness, in sharp contrast to models with size-independent stochastic growth rates. Second, we show that the steady-state population cell-size distribution is obtained from the deterministic inverse-square-law solution by a one-sided exponential convolution with kernel width set by the strength of growth fluctuations. Third, the mean-rescaled population size $N_t/\left\langle N_t\right\rangle$ converges to a stationary compound Poisson-exponential distribution that depends only on growth noise. This distribution, and hence the long-time shape of population-size fluctuations, is unchanged by division-size noise or asymmetric partitioning. These results identify Feller-type exponential growth with square-root noise as an exactly solvable benchmark for stochastic growth in size-controlled populations and provide concrete signatures that distinguish it from models with size-independent growth-rate noise.
We study the emergence of chaos in multilevel atoms with all-to-all interactions, inspired by cavity QED. Focusing on a 3-level Tavis-Cummings model in a far detuned limit, we detail its deep Hilbert space structure -- i.e. we enumerate all distinct dynamical sectors, beyond the totally symmetric subspace -- by using the Schur-Weyl duality, which is applicable thanks to the permutation symmetry in the all-to-all Hamiltonian. Strong Hilbert space fragmentation ensues from the non-abelian nature of the symmetry, with some sectors displaying regular dynamics and others being chaotic. We uncover that many permutation symmetry sectors contribute to the dynamics in the classical limit, in addition to the commonly studied totally symmetric subspace. To elucidate the dynamical responses in each of the symmetry sectors, we propose a semiclassical description in terms of spin coherent states, which is also able to explain the origin of chaotic or regular dynamics with a simple geometrical argument. Our work contributes to the study of the quantum-classical correspondence in chaotic systems, and uncovers a rich structure in multilevel all-to-all interacting models.
Gross-Neveu model in 2+1 dimensions exhibits a continuous transition from gapless Dirac semimetal to the gapped quantum anomalous Hall (QAH) insulator at a finite (attractive) coupling, at which the inversion and time-reversal symmetry become spontaneously broken, and the flavor O($M$) symmetry remains preserved. A unification of leading order parameters of 2+1 dimensional $N$ four-component Dirac fermions collects all Lorentz-singlet mass-like fermion bilinears, except the one condensing in the QAH state, into an irreducible representation of the O($M=4N$), and predicts another phase transition in the Gross-Neveu model to occur at a strong (repulsive) coupling. Here, a fermionic auxiliary-field quantum Monte Carlo algorithm is employed in order to study a lattice realization of the Gross-Neveu field theory in the repulsive regime, where the sign problem is absent. We indeed find the O($4N$) symmetry breaking transition out of Dirac semimetal to occur and to be weakly first-order for $N=2$, relevant to graphene. The size of the discontinuity and the magnitude of the critical coupling, however, both grow with $N$. Adding a finite chemical potential is found to break the symmetry and cause superconductivity. These results are in broad agreement with the predictions of the unified field theory. Our lattice model also displays an interesting exact O($2N$) symmetry, a subgroup of the low-energy O($4N$), and has the ordered ground state with the order parameter that belongs to its $N(2N-1)$-dimensional representation. Other order parameters are also examined, and a certain hierarchy among those that belong to different representations of the exact $O(2N)$ is observed.