We introduce a numerical method to extract the parameters of run-and-tumble dynamics from experimental measurements of the intermediate scattering function. We show that proceeding in Laplace space is unpractical and employ instead renewal processes to work directly in real time. We first validate our approach against data produced using agent-based simulations. This allows us to identify the length and time scales required for an accurate measurement of the motility parameters, including tumbling frequency and swim speed. We compare different models for the run-and-tumble dynamics by accounting for speed variability at the single-cell and population level, respectively. Finally, we apply our approach to experimental data on wild-type Escherichia coli obtained using differential dynamic microscopy.

We show that the second law of thermodynamics poses a restriction on how well we can discriminate between quantum states. By examining an ideal gas with a quantum internal degree of freedom undergoing a cycle based on a proposal by Peres, we establish a nontrivial upper bound on the attainable accuracy of quantum state discrimination. This thermodynamic bound, which relies solely on the linearity of quantum mechanics and the constraint of no work extraction, matches Holevo's bound on accessible information, but is looser than the Holevo-Helstrom bound. The result gives more evidence on the disagreement between thermodynamic entropy and von Neumann entropy and places potential limitations on proposals beyond quantum mechanics.

Liquid drainage through foam is dominated by gravity, capillary, and viscous forces. The liquid is conducted by an isotropic network of Plateau borders; however, imposed stress changes the alignment of the foam's structural elements. Previous numerical simulations predicted that a vertical drainage flow will be deflected horizontally if the foam is sheared. We investigated this phenomenon by measuring the distribution of the liquid fraction within a foam formed in a flat rectangular cell. The foam was subjected to shear stress under a forced liquid supply from the top of the cell. Neutron radiographies of unchanged and sheared foam were analyzed to extract measurements of the liquid fraction. Deflections in the distribution of the drainage liquid were detected and found to be positively correlated with increasing foam shear.

The unique ability of fire ants to form a raft to survive flooding rain has enchanted biologists as well as researchers in other disciplines. It was established during the last decade that a three-dimensional aggregation of fire ants exhibits viscoelasticity with respect to external compression and shearing among numerous unusual mechanical properties. Continuing these works, we will study the ant raft in its natural form, i.e., composing no more than two layers. This allowed us to focus on the cracks that are unique to membranes and see how their patterns are influenced by the fact that these ants are mobile and can self-repair the damage to keep their raft from disintegration. In the beginning, we show that vertical and horizontal shaking can also prompt fire ants to aggregate. The canonical view that the stability of ant raft relies on the Cheerios effect and a combination of other parameters is tested. The force-displacement experiment is performed to show that two distinct mechanical responses and fracture patterns, characteristic of ductile and brittle materials, can be elicited, depending on the magnitude of the pull speed. During the process, we counted the number of ants that actively participated in the stress-strain relation and used this information to roughly sketch out the force chain. The latter information reveals that the pull force expedites the alignment of fire ants, in analogy to the effect of an electric field on liquid crystal polymers. To highlight the self-healing nature, we employ the creep experiment to study how the length and Young's modulus of the raft change or relax with time. One major finding is that the raft can exhibit zero Poisson's ratio without resorting to specific geometry structures. This is enabled by the active recruitment of ants from the top layer to the bottom layer to keep the raft from disintegrating.

We simulate, using a particle-in-cell code, the chain of acceleration processes at work during the Compton-based interaction of a dilute electron-ion plasma with an extreme-intensity, incoherent $\gamma$-ray flux with a photon density several orders of magnitude above the particle density. The plasma electrons are initially accelerated in the radiative flux direction through Compton scattering. In turn, the charge-separation field from the induced current drives forward the plasma ions to near-relativistic speed and accelerates backwards the nonscattered electrons to energies easily exceeding those of the driving photons. The dynamics of those energized electrons is determined by the interplay of electrostatic acceleration, bulk plasma motion, inverse Compton scattering and deflections off the mobile magnetic fluctuations generated by a Weibel-type instability. The latter Fermi-like effect notably gives rise to a forward-directed suprathermal electron tail. We provide simple analytical descriptions for most of those phenomena and examine numerically their sensitivity to the parameters of the problem.

The $\pm J$ Ising model is a simple frustrated spin model, where the exchange couplings independently take the discrete value $-J$ with probability $p$ and $+J$ with probability $1-p$. It is especially appealing due to its connection to quantum error correcting codes. Here, we investigate the nonequilibrium critical behavior of the two-dimensional $\pm J$ Ising model, after a quench from different initial conditions to a critical point $T_c\left(p\right)$ on the paramagnetic-ferromagnetic (PF) transition line, especially above, below, and at the multicritical Nishimori point (NP). The dynamical critical exponent $z_c$ seems to exhibit nonuniversal behavior for quenches above and below the NP, which is identified as a preasymptotic feature due to the repulsive fixed point at the NP, whereas for a quench directly to the NP, the dynamics reaches the asymptotic regime with $z_c\simeq 6.02\left(6\right)$. We also consider the geometrical spin clusters (of like spin signs) during the critical dynamics. Each universality class on the PF line is uniquely characterized by the stochastic Loewner evolution with corresponding parameter $\kappa$. Moreover, for the critical quenches from the paramagnetic phase, the model, irrespective of the frustration, exhibits an emergent critical percolation topology at the large length scales.

We introduce modifications to Monte Carlo simulations of the Feynman path integral that improve sampling of localized interactions. The algorithms generate trajectories in simple background potentials designed to concentrate them around the interaction region, reminiscent of importance sampling. This improves statistical sampling of the system and overcomes a long-time undersampling problem caused by the spatial diffusion inherent in Brownian motion. We prove the validity of our approach using previous analytic work on the distribution of values of the Wilson line over path integral trajectories and illustrate the improvements on some simple quantum mechanical systems.

Continuous attractor neural networks (CANN) form an appealing conceptual model for the storage of information in the brain. However a drawback of CANN is that they require finely tuned interactions. We here study the effect of quenched noise in the interactions on the coding of positional information within CANN. Using the replica method we compute the Fisher information for a network with position-dependent input and recurrent connections composed of a short-range (in space) and a disordered component. We find that the loss in positional information is small for not too large disorder strength, indicating that CANN have a regime in which the advantageous effects of local connectivity on information storage outweigh the detrimental ones. Furthermore, a substantial part of this information can be extracted with a simple linear readout.

We study the ergodic properties of one-dimensional Brownian motion with resetting. Using generic classes of statistics of times between resets, we find respectively for thin- or fat-tailed distributions the normalized or non-normalized invariant density of this process. The former case corresponds to known results in the resetting literature and the latter to infinite ergodic theory. Two types of ergodic transitions are found in this system. The first is when the mean waiting time between resets diverges, when standard ergodic theory switches to infinite ergodic theory. The second is when the mean of the square root of time between resets diverges and the properties of the invariant density are drastically modified. We then find a fractional integral equation describing the density of particles. This finite time tool is particularly useful close to the ergodic transition where convergence to asymptotic limits is logarithmically slow. Our study implies rich ergodic behaviors for this nonequilibrium process which should hold far beyond the case of Brownian motion analyzed here.

Craters formed by the impact of agglomerated materials are commonly observed in nature, such as asteroids colliding with planets and moons. In this paper, we investigate how the projectile spin and cohesion lead to different crater shapes. For that, we carried out discrete element method computations of spinning granular projectiles impacting onto cohesionless grains for different bonding stresses, initial spins, and initial heights. We found that, as the bonding stresses decrease and the initial spin increases, the projectile's grains spread farther from the collision point, and in consequence, the crater shape becomes flatter, with peaks around the rim and in the center of the crater. Our results shed light on the dispersion of the projectile's material and the different shapes of craters found on Earth and other planetary environments.

Conical surfaces pose an interesting challenge to crystal growth: A crystal growing on a cone can wrap around and meet itself at different radii. We use a disk-packing algorithm to investigate how this closure constraint can geometrically frustrate the growth of single crystals on cones with small opening angles. By varying the crystal seed orientation and cone angle, we find that—except at special commensurate cone angles—crystals typically form a seam that runs along the axial direction of the cone, while near the tip, a disordered particle packing forms. We show that the onset of disorder results from a finite-size effect that depends strongly on the circumference and not on the seed orientation or cone angle. This finite-size effect occurs also on cylinders, and we present evidence that on both cylinders and cones, the defect density increases exponentially as circumference decreases. We introduce a simple model for particle attachment at the seam that explains the dependence on the circumference. Our findings suggest that the growth of single crystals can become frustrated even very far from the tip when the cone has a small opening angle. These results may provide insights into the observed geometry of conical crystals in biological and materials applications.

We study the effect of spatially varying potential and diffusivity on the dispersion of a tracer particle in single-file diffusion. Noninteracting particles in such a system exhibit normal diffusion at late times, which is characterized by an effective diffusion constant $D_{\mathrm{eff}}$. Here we demonstrate the physically appealing result that the dispersion of single-file tracers in this system has the same long-time behavior as that for Brownian particles in a spatially homogeneous system with constant diffusivity $D_{\mathrm{eff}}$. Our results are based on a late-time analysis of the Fokker-Planck equation, motivated by the mathematical theory of homogenization. The findings are confirmed by numerical simulations for both annealed and quenched initial conditions.

Through numerous experiments that analyzed rare event statistics in heterogeneous media, it was discovered that in many cases the probability density function for particle position, $P(X,t)$, exhibits a slower decay rate than the Gaussian function. Typically, the decay behavior is exponential, referred to as Laplace tails. However, many systems exhibit an even slower decay rate, such as power-law, log-normal, or stretched exponential. In this study, we utilize the continuous-time random walk method to investigate the rare events in particle hopping dynamics and find that the properties of the hop size distribution induce a critical transition between the Laplace universality of rare events and a more specific, slower decay of $P(X,t)$. Specifically, when the hop size distribution decays slower than exponential, such as $e^{-{\left|x\right|}^{\beta }}$ ($\beta >1$), the Laplace universality no longer applies, and the decay is specific, influenced by a few large events, rather than by the accumulation of many smaller events that give rise to Laplace tails.

Time-dependent protocols that perform irreversible logical operations, such as memory erasure, cost work and produce heat, placing bounds on the efficiency of computers. Here we use a prototypical computer model of a physical memory to show that it is possible to learn feedback-control protocols to do fast memory erasure without input of work or production of heat. These protocols, which are enacted by a neural-network “demon,” do not violate the second law of thermodynamics because the demon generates more heat than the memory absorbs. The result is a form of nonlocal heat exchange in which one computation is rendered energetically favorable while a compensating one produces heat elsewhere, a tactic that could be used to rationally design the flow of energy within a computer.

We investigate the influence of long-range (LR) interactions on the phase ordering dynamics of the one-dimensional random-field Ising model (RFIM). Unlike the usual RFIM, a spin interacts with all other spins through a ferromagnetic coupling that decays as $r^{-(1+\sigma )}$, where $r$ is the distance between two spins. In the absence of LR interactions, the size of coarsening domains $R\left(t\right)$ exhibits a crossover from pure system behavior $R\left(t\right)\sim t^{1/2}$ to an asymptotic regime characterized by logarithmic growth: $R\left(t\right)\sim {(lnt)}^2$. The LR interactions affect the preasymptotic regime, which now exhibits ballistic growth $R\left(t\right)\sim t$, followed by $\sigma$-dependent growth $R\left(t\right)\sim t^{1/(1+\sigma )}$. Additionally, the LR interactions also affect the asymptotic logarithmic growth, which becomes $R\left(t\right)\sim {(lnt)}^{\alpha \left(\sigma \right)}$ with $\alpha \left(\sigma \right)<2$. Thus, LR interactions lead to faster growth than for the nearest-neighbor system at short times. Unexpectedly, this driving force causes a slowing down of the dynamics ($\alpha <2$) in the asymptotic logarithmic regime. This is explained in terms of a nontrivial competition between the pinning force caused by the random field and the driving force introduced by LR interactions. We also study the spatial correlation function and the autocorrelation function of the magnetization field. The former exhibits superuniversality for all $\sigma$, i.e., a scaling function that is independent of the disorder strength. The same holds for the autocorrelation function when $\sigma <1$, whereas a signature of the violation of superuniversality is seen for $\sigma >1$.

We study the critical behavior of three-dimensional (3D) lattice Abelian Higgs (AH) gauge models with noncompact gauge variables and multicomponent complex scalar fields, along the transition line between the Coulomb and Higgs phases. Previous works that focused on gauge-invariant correlations provided evidence that, for a sufficiently large number of scalar components, these transitions are continuous and associated with the stable charged fixed point of the renormalization-group flow of the 3D AH field theory (scalar electrodynamics), in which charged scalar matter is minimally coupled with an electromagnetic field. Here we extend these studies by considering gauge-dependent correlations of the gauge and matter fields, in the presence of two different gauge fixings, the Lorenz and the axial gauge fixing. Our results for $N=25$ are definitely consistent with the predictions of the AH field theory and therefore provide additional evidence for the characterization of the 3D AH transitions along the Coulomb-Higgs line as charged transitions in the AH field-theory universality class. Moreover, our results give additional insights on the role of the gauge fixing at charged transitions. In particular, we show that scalar correlations are critical only if a hard Lorenz gauge fixing is imposed.

The limits on a system's response to external perturbations inform our understanding of how physical properties can be shaped by microscopic characteristics. Here, we derive constraints on the steady-state nonequilibrium response of physical observables in terms of the topology of the microscopic state space and the strength of thermodynamic driving. Notably, evaluation of these limits requires no kinetic information beyond the state-space structure. When applied to models of receptor binding, we find that sensitivity is bounded by the steepness of a Hill function with a Hill coefficient enhanced by the chemical driving beyond the structural equilibrium limit.

We generate nonlattice packings of spheres in up to 22 dimensions using the geometrical constraint satisfaction algorithm RRR. Our aggregated data suggest that it is easy to double the density of Ball's lower bound and, more tentatively, that the exponential decay rate of the density can be improved relative to Minkowski's longstanding 1/2.

We investigate the magnetic fluctuations in a mesoscopic critical region formed at the interface due to smooth time-independent spatial variations of a control parameter around its critical value. In the proximity of the spatial critical point, the order parameter fluctuations exhibit a mesoscopic nature, characterized by their significant size compared to the lattice constant, while gradually decaying away from the critical region. To explain this phenomenon, we present a minimal model that effectively captures this behavior and demonstrates its connection to the integrable Painlevé-II equation governing the local order parameter. By leveraging the well-established mathematical properties of this equation, we gain valuable insights into the nonlinear susceptibilities exhibited within this region.

In the past decade, great strides have been made to quantify the dynamics of single-cell growth and division in microbes. In order to make sense of the evolutionary history of these organisms, we must understand how features of single-cell growth and division influence evolutionary dynamics. This requires us to connect processes on the single-cell scale to population dynamics. Here, we consider a model of microbial growth in finite populations which explicitly incorporates the single-cell dynamics. We study the behavior of a mutant population in such a model and ask: can the evolutionary dynamics be coarse-grained so that the forces of natural selection and genetic drift can be expressed in terms of the long-term fitness? We show that it is in fact not possible, as there is no way to define a single fitness parameter (or reproductive rate) that defines the fate of an organism even in a constant environment. This is due to fluctuations in the population averaged division rate. As a result, various details of the single-cell dynamics affect the fate of a new mutant independently from how they affect the long-term growth rate of the mutant population. In particular, we show that in the case of neutral mutations, variability in generation times increases the rate of genetic drift, and in the case of beneficial mutations, variability decreases its fixation probability. Furthermore, we explain the source of the persistent division rate fluctuations and provide analytic solutions for the fixation probability as a multispecies generalization of the Euler-Lotka equation.