The autoresonant approach to generation of solitary structures in Bose-Einstein condensates by chirped frequency space-time modulation of the interaction strength is proposed. Both a spatially periodic case and a finite-size trap are studied numerically within a Gross-Pitaevskii equation. Weakly nonlinear theory of the process is developed in the spatially periodic case using Whitham's averaged variational principle. The theory also describes the threshold phenomenon setting the lowest bound on the amplitude of modulations of the interaction strength for autoresonant excitation.

The emergence of random matrix spectral correlations in interacting quantum systems is a defining feature of quantum chaos. We study such correlations in terms of the spectral form factor and its moments in interacting chaotic few- and many-body systems, modeled by suitable random-matrix ensembles. We obtain the spectral form factor exactly for large Hilbert space dimension. Extrapolating those results to finite Hilbert space dimension we find a universal transition from the noninteracting to the strongly interacting case, which can be described as a simple combination of these two limits. This transition is governed by a single scaling parameter. In the bipartite case we derive similar results also for all moments of the spectral form factor. We confirm our results by extensive numerical studies and demonstrate that they apply to more realistic systems given by a pair of quantized kicked rotors as well. Ultimately we complement our analysis by a perturbative approach covering the small-coupling regime.

From a context of evolutionary dynamics, social games can be studied as complex systems that may converge to a Nash equilibrium. Nonetheless, they can behave in an unpredictable manner when looking at the spatial patterns formed by the agents' strategies. This is known in the literature as spatial chaos. In this paper we analyze the problem for a deterministic prisoner's dilemma and a public goods game and calculate the Hamming distance that separates two solutions that start at very similar initial conditions for both cases. The rapid growth of this distance indicates the high sensitivity to initial conditions, which is a well-known indicator of chaotic dynamics.

We study the two-dimensional (2D) Ising model in a complex magnetic field in the vicinity of the Yang-Lee edge singularity. By using Baxter's variational corner transfer matrix method combined with analytic techniques, we numerically calculate the scaling function and obtain an accurate estimate of the location of the Yang-Lee singularity. The existing series expansions for susceptibility of the 2D Ising model on a triangular lattice by Chan, Guttmann, Nickel, and Perk allowed us to substantially enhance the accuracy of our calculations. Our results are in excellent agreement with the Ising field theory calculations by Fonseca, Zamolodchikov, and the recent work by Xu and Zamolodchikov. In particular, we numerically confirm an agreement between the leading singular behavior of the scaling function and the predictions of the ${\mathcal{M}}_{2/5}$ conformal field theory.

Terahertz (THz) emissions from fast electron and ion currents driven in relativistic, femtosecond laser-foil interactions are examined theoretically. We first consider the radiation from the energetic electrons exiting the backside of the target. Our kinetic model takes account of the coherent transition radiation due to these electrons crossing the plasma-vacuum interface as well as of the synchrotron radiation due to their deflection and deceleration in the sheath field they set up in vacuum. After showing that both mechanisms tend to largely compensate each other when all the electrons are pulled back into the target, we investigate the scaling of the net radiation with the sheath field strength. We then demonstrate the sensitivity of this radiation to a percent-level fraction of escaping electrons. We also study the influence of the target thickness and laser focusing. The same sheath field that confines most of the fast electrons around the target rapidly sets into motion the surface ions. We describe the THz emission from these accelerated ions and their accompanying hot electrons by means of a plasma expansion model that allows for finite foil size and multidimensional effects. Again, we explore the dependencies of this radiation mechanism on the laser-target parameters. Under conditions typical of current ultrashort laser-solid experiments, we find that the THz radiation from the expanding plasma is much less energetic—by one to three orders of magnitude—than that due to the early-time motion of the fast electrons.

We study via Monte Carlo simulations the influence of quenched and mobile impurities in the contact process (CP) on two-dimensional lattice and continuum systems. In the lattice system, the effect of mobile impurity was studied for the density $n_i=0.2$ and two selected values of hopping probability for impurity particles, $w=0.5$ and 1. In the continuum system, the CP was defined by distributing spherical impurity particles of diameter ${\sigma }_i$ and number density $n_i=0.2$ and active particles of diameter unity and number density $1-n_i$ on a square substrate with periodic boundaries. In each dynamic process, a particle is selected at random; the active particle either creates with a rate $\lambda$ an offspring at a distance $r$ $(1\le r\le 1.5)$ from the active particle or annihilates with a unit rate, and the impurity particle hops a distance $r$ $(0\le r\le 1)$, both along randomly selected directions. We found that the lattice CP shows power-law behaviors with varying critical exponents depending on the values of $w$. For the continuum CP with quenched impurity, the critical behavior followed the activated scaling scenario, whereas with mobile impurity usual power-law behaviors were observed but the critical exponents varied depending on the values of ${\sigma }_i$.

We study a population of swarmalators, mobile variants of phase oscillators, which run on a ring and have both attractive and repulsive interactions. This one-dimensional (1D) swarmalator model produces several of collective states: the standard sync and async states as well as a splaylike “polarized” state and several unsteady states such as active bands or swirling. The model's simplicity allows us to describe some of the states analytically. The model can be considered as a toy model for real-world swarmalators such as vinegar eels and sperm which swarm in quasi-1D geometries.

Bifurcation phenomena are common in multidimensional multiparameter dynamical systems. Normal form theory suggests that bifurcations are driven by relatively few combinations of parameters. Models of complex systems, however, rarely appear in normal form, and bifurcations are controlled by nonlinear combinations of the bare parameters of differential equations. Discovering reparameterizations to transform complex equations into a normal form is often very difficult, and the reparameterization may not even exist in a closed form. Here we show that information geometry and sloppy model analysis using the Fisher information matrix can be used to identify the combination of parameters that control bifurcations. By considering observations on increasingly long timescales, we find those parameters that rapidly characterize the system's topological inhomogeneities, whether the system is in normal form or not. We anticipate that this novel analytical method, which we call time-widening information geometry (TWIG), will be useful in applied network analysis.

Microtubules are dynamic intracellular fibers that have been observed experimentally to undergo spontaneous self-alignment. We formulate a three-dimensional (3D) mean-field theory model to analyze the nematic phase transition of microtubules growing and interacting within a 3D space, then make a comparison with computational simulations. We identify a control parameter $G_{\text{eff}}$ and predict a unique critical value $G_{\text{eff}}=1.56$ for which a phase transition can occur. Furthermore, we show both analytically and using simulations that this predicted critical value does not depend on the presence of zippering. The mean-field theory developed here provides an analytical estimate of microtubule patterning characteristics without running time-consuming simulations and is a step towards bridging scales from microtubule behavior to multicellular simulations.

We study numerically a spherical particle settling through a density transition layer at moderate Reynolds numbers ${\text{Re}}_u=69\sim 259$ for the upper fluid. We investigate how the transition layer thickness affects the particle's bouncing behavior as it crosses the interface. The previous intuitive understanding was that the bounce occurs when the relative thickness of the transition layer, $L/D$, which is characterized by the ratio of the layer thickness $L$ to the particle diameter $D$, is small. Indeed, we report no bounce phenomenon for very thick interfaces, i.e., $L/D>10$ in the current parametric range. However, we argue that the bounce can also be inhibited when $L/D$ is too small. Upon a fixed upper layer Reynolds number ${\text{Re}}_u=207$ with varying $L/D$, we examine the flow evolution of these cases. We propose that this inhibition is attributed to two mechanisms. First, as the interface thickness decreases, the detachment of the attached lighter fluid from the upper layer occurs more rapidly, resulting in a faster decrease in buoyancy. Second, in the case of a very thin interface $(L/D=0.5–3.0)$, the residual light fluid accumulates and undergoes a secondary detachment, separating from the particle at an angle relative to the central axis. This secondary detachment reduces the drag force and effectively prevents the particle from experiencing a rebound motion.

We relate the model of vorticity wave interaction to Krein collision, $\mathcal{PT}$-symmetry breaking, and the formation of exceptional points in shear flow instabilities. We show that the dynamical system of coupled vorticity waves is a pseudo-Hermitian system with nonreciprocal coupling terms. Krein signatures of the eigenvalues are illustrated as the signs of the action of the vorticity waves. Interaction between positive-action and negative-action vorticity waves then corresponds to the Krein collision between eigenvalues with opposite Krein signatures, the spontaneous breaking of $\mathcal{PT}$ symmetry, and the formation of exceptional points. The control parameter of the $\mathcal{PT}$-symmetry-breaking bifurcation is the ratio between frequency detuning and coupling strength of the vorticity waves. The critical behavior near the exceptional points is described as a transition between phase-locking and phase-slip dynamics of the vorticity waves. The phase-slip dynamics correspond to nonmodal, transient growth of perturbations in the regime of unbroken $\mathcal{PT}$ symmetry, and the phase-slip frequency $\mathrm{\Omega }\propto |k-k_c{|}^{1/2}$ shares the same critical exponent with the phase rigidity of system eigenvectors.

We characterize thermalization slowing down of Josephson junction networks in one, two, and three spatial dimensions for systems with hundreds of sites by computing their entire Lyapunov spectra. The ratio of Josephson coupling $E_J$ to energy density $h$ controls two different universality classes of thermalization slowing down, namely, the weak-coupling regime, $E_J/h\to 0$, and the strong-coupling regime, $E_J/h\to \infty$. We analyze the Lyapunov spectrum by measuring the largest Lyapunov exponent and by fitting the rescaled spectrum with a general ansatz. We then extract two scales: the Lyapunov time (inverse of the largest exponent) and the exponent for the decay of the rescaled spectrum. The two universality classes, which exist irrespective of network dimension, are characterized by different ways the extracted scales diverge. The universality class corresponding to the weak-coupling regime allows for the coexistence of chaos with a large number of near-conserved quantities and is shown to be characterized by universal critical exponents, in contrast with the strong-coupling regime. We expect our findings, which we explain using perturbation theory arguments, to be a general feature of diverse Hamiltonian systems.

Chemical reactions involve the movement of charges, and this paper presents a mathematical model for describing chemical reactions in electrolytes. The model is developed using an energy variational method that aligns with classical thermodynamics principles. It encompasses both electrostatics and chemical reactions within consistently defined energetic and dissipative functionals. Furthermore, the energy variation method is extended to account for open systems that involve the input and output of charge and mass. Such open systems have the capability to convert one form of input energy into another form of output energy. In particular, a two-domain model is developed to study a reaction system with self-regulation and internal switching, which plays a vital role in the electron transport chain of mitochondria responsible for ATP generation—a crucial process for sustaining life. Simulations are conducted to explore the influence of electric potential on reaction rates and switching dynamics within the two-domain system. It shows that the electric potential inhibits the oxidation reaction while accelerating the reduction reaction.

Chromonic nematics are lyotropic liquid crystals that have already been known for half a century, but have only recently raised interest for their potential applications in life sciences. Determining elastic constants and anchoring strengths for rigid substrates has thus become a priority in the characterization of these materials. Here we present a method to determine chromonics' planar anchoring strength. We call it geometric as it is based on recognition and fitting of the stable equilibrium shapes of droplets surrounded by the isotropic phase in a thin cell with plates enforcing parallel alignments of the nematic director. We apply our method to shapes observed in experiments; they resemble elongated rods with round ends, which are called bâtonnets. Our theory also predicts other droplets' equilibrium shapes, which are either slender and round, called discoids, or slender and pointed, called tactoids. In particular, sufficiently small droplets are expected to display shape bistability, with two equilibrium shapes, one tactoid and one discoid, exchanging roles as stable and metastable shapes upon varying their common area.

We report experimental results on breath figures (BFs) observed on substrates with quenched disorder. The evolution of BFs is found to be primarily influenced by global parameters associated with boundary conditions. We investigate classical statistical measures and explore topological properties using persistent homology techniques based on a modified Vietoris-Rips complex. Our findings reveal that the evolution of the number surface density of condensed droplets plays a crucial role in determining various condensation stages previously considered distinct. This evolution is significantly influenced by the distribution of nucleation sites and the individual growth law governing water droplets when coalescence does not occur. Ultimately, we demonstrate the capability to predict coalescence events based on the topological characteristics of BFs at a given point in time.

By producing localized wave regions at the ends of an open-field-line magnetic confinement system, ponderomotive walls can be used to differentially confine different species in the plasma. Furthermore, if the plasma is rotating, this wall can be magnetostatic in the laboratory frame, resulting in simpler engineering and better power flow. However, recent work on such magnetostatic walls has shown qualitatively different potentials than those found in the earlier, nonrotating theory. Here, using a simple slab model of a ponderomotive wall, we resolve this discrepancy. We show that the form of the ponderomotive potential in the comoving plasma frame depends on the assumption made about the electrostatic potential in the laboratory frame. If the laboratory-frame potential is unperturbed by the magnetic oscillation, one finds a parallel-polarized wave in the comoving frame, while if each field line remains equipotential throughout the perturbation region, one finds a perpendicularly polarized wave. This in turn dramatically changes the averaged ponderomotive force experienced by a charged particle along the field line, not only its scaling, but also its direction.

We present a simple two-dimensional model for a phase transition, then study its predictions, in particular the memory properties. The direct transformation is modeled by randomly placing small squares, “nuclei”, on an initially empty surface. Then, the nuclei expand (“grow”) up to finite final sizes which are randomly chosen in a given range, while keeping their square shape. An important issue is the “interaction” which forces some squares to remain at smaller sizes if the surrounding squares get in the way of their growth. Interestingly, this naturally leads to quasiequal total area covered by the squares of each size after a complete direct transformation. Next, it is shown that the system “remembers” incomplete (“arrested”) reverse transformations taking place in reversed order of the squares sizes. The memory is “encrypted” in the distribution of the squares sizes after a next direct transformation and manifests as a significant imbalance between the areas covered by the “big” and “small” (relative to the arrest size) squares. We are able to also reproduce the so-called “hammer effect” and the memorizing of multiple arrest points. Our model is particularly relevant for the thermal memory effect in shape memory alloys, and we actually borrowed many features from existing thermodynamic models addressing this effect. However, here we eliminate the explicit thermodynamics and end up with a statistical geometry model, presumably easier to reproduce.

We present the chiral knife edge rattleback, an alternative version of previously presented systems that exhibit spin inversion. We offer a full treatment of the model using qualitative arguments, analytical solutions as well as numerical results. We treat a reduced, one–mode problem which not only contains the essence of the physics of spin inversion, but that also exhibits an unexpected connection to the Chaplygin sleigh, providing insight into the nonholonomic structure of the problem. We also present exact results for the full problem together with estimates of the time between inversions that agree with previous results in the literature.

The brain is a nonequilibrium system whose dynamics change in different brain states, such as wakefulness and deep sleep. Thermodynamics provides the tools for revealing these nonequilibrium dynamics. We used violations of the fluctuation-dissipation theorem to describe the hierarchy of nonequilibrium dynamics associated with different brain states. Together with a whole-brain model fitted to empirical human neuroimaging data, and deriving the appropriate analytical expressions, we were able to capture the deviation from equilibrium in different brain states that arises from asymmetric interactions and hierarchical organization.

Eukaryotic cells can polarize and migrate in response to electric fields via “galvanotaxis,” which aids wound healing. Experimental evidence suggests cells sense electric fields via molecules on the cell's surface redistributing via electrophoresis and electroosmosis, though the sensing species has not yet been conclusively identified. We develop a model that links sensor redistribution and galvanotaxis using maximum likelihood estimation. Our model predicts a single universal curve for how galvanotactic directionality depends on field strength. We can collapse measurements of galvanotaxis in keratocytes, neural crest cells, and granulocytes to this curve, suggesting that stochasticity due to the finite number of sensors may limit galvanotactic accuracy. We find cells can achieve experimentally observed directionalities with either a few ($\sim 100$) highly polarized sensors or many ($\sim {10}^4$) sensors with an $\sim 6–10\%$ change in concentration across the cell. We also identify additional signatures of galvanotaxis via sensor redistribution, including the presence of a tradeoff between accuracy and variance in cells being controlled by rapidly switching fields. Our approach shows how the physics of noise at the molecular scale can limit cell-scale galvanotaxis, providing important constraints on sensor properties and allowing for new tests to determine the specific molecules underlying galvanotaxis.