In this paper, we conduct experimental investigations on the behavior of confined self-propelled particles within a circular arena, employing small commercial robots capable of locomotion, communication, and information processing. These robots execute circular trajectories, which can be clockwise or counterclockwise, based on two internal states. Using a majority-based stochastic decision algorithm, each robot can reverse its direction based on the states of two neighboring robots. By manipulating a control parameter governing the interaction, the system exhibits a transition from a state where all robots rotate randomly to one where they rotate uniformly in the same direction. Moreover, this transition significantly impacts the trajectories of the robots. To extend our findings to larger systems, we introduce a mathematical model enabling characterization of the order transition type and the resulting trajectories. Our results reveal a second-order transition from active Brownian to chiral motion.

We investigate experimentally the collective motion of polar vibrated disks in an annular geometry, varying both the packing fraction and the amplitude of the angular noise. For low enough noise and large enough density, an overall collective motion takes place along the tangential direction. The spatial organization of the flow reveals the presence of polar bands of large density, as expected from the commonly accepted picture of the transition to collective motion in systems of aligning polar active particles. However, in our case, the low density phase is also polar, consistent with what is observed when jamming takes place in a very high density flock. Interestingly, while in that case the particles in the high density bands are arrested, resulting in an upstream propagation at a constant speed, in our case the bands travel downstream with a density-dependent speed. We demonstrate from local measurements of the packing fraction, alignment, and flow speeds that the bands observed here result both from a polar ordering process and a motility induced phase separation mechanism.

Nematic tactoids are spindle-shaped droplets of a nematic phase nucleated in the co-existing isotropic phase. According to equilibrium theory, their internal structure and shape are controlled by a balance between the elastic deformation of the director field, induced by the preferred anchoring of that director field to the interface, and the interfacial free energy. Recent experiments on tactoids of chitin nanocrystals dispersed in water show that electrical fields can very strongly elongate tactoids, at least if the tactoids are sufficiently large in volume. However, this observation contradicts the predictions of equilibrium theory as well as findings from Monte Carlo simulations that do not show this kind of extreme elongation to take place at all. To explain this, we put forward a relaxational model based on the Oseen-Frank free energy of elastic deformation of a director field coupled to an anisotropic surface free energy. In our model, we use two reaction coordinates to describe the director field and the extent of elongation of the droplets and evaluate the evolution of both as a function of time following the switching on of an electric field. Depending on the relative magnitude of the fundamental relaxation rates associated with the two reaction coordinates, we find that the aspect ratio of the drops may develop a large and very long-lived overshoot before eventually relaxing to the much smaller equilibrium value. In that case, the response of the curvature of the director field lags behind, explaining the experimental observations. Our theory describes the experimental data reasonably well.

Functionalized thin elastic films and membranes frequently feature internal sites of net forces or stresses. These are, for instance, active sites of actuation, or rigid inclusions in a strained membrane that induce counterstress upon externally imposed deformations. We theoretically analyze the geometry of isotropic, flat, thin, linearly elastic films or membranes of finite thickness, laterally extended to infinity. At the mathematical core of such characterizations are the fundamental solutions for localized force and stress singularities associated with corresponding Green's functions. We derive such solutions in three dimensions and place them into the context of previous two-dimensional calculations. To this end, we consider both no-slip and stress-free conditions at the top and/or bottom surfaces. We provide an understanding for why the fully free-standing thin elastic membrane leads to diverging solutions in most geometries and compare these situations to the truly two-dimensional case. A no-slip support of at least one of the surfaces stabilizes the solution, which illustrates that the divergences in the fully free-standing case are connected to global deformations. Within the aforementioned framework, our results are important for associated theoretical characterizations of thin elastic films, whether supported or free-standing, and of membranes subject to internal or external forces or stresses.

Polymer solutions in the dilute regime play a significant role in industrial applications. Due to the intricate rheological properties of these highly viscoelastic fluids, especially in complex flow geometries, a thorough numerical analysis of their flow dynamics is imperative. In this research, we present a numerical investigation of purely elastic instability occurring in two- and three-dimensional serpentine channels under conditions where fluid inertia is negligible and across a broad spectrum of polymer relaxation times. Our findings reveal a strong qualitative agreement between the existing experimental results obtained from dilute solutions of flexible polymers in microfluidic devices and the numerical simulations conducted in two and three dimensions using the Oldroyd-B model. Spatial flow observations and statistical analysis of temporal flow features indicate that this purely elastic turbulent flow exhibits nonhomogeneous, non-Gaussian, and anisotropic characteristics across all scales. Additionally, our comparison of two- and three-dimensional simulations demonstrates that the elastic instability is primarily driven by the curvature of the streamlines induced by the flow geometry, rather than the weak secondary flow in the azimuthal direction. Therefore, our two-dimensional numerical simulations successfully replicate, at least qualitatively, the features observed in three-dimensional experiments. Furthermore, spectral analysis suggests that, in comparison to elastic turbulence in the dilute regime, the range of scales for the excited fluctuations is narrower.

Complicated many-body interactions between ions and surrounding particles exist in warm and hot dense plasmas. It will significantly alter the atomic structures and dynamic properties of the embedded ions. Recently, the atomic-state-dependent (ASD) screening model has been proposed and shown to be valid for investigating the screening effect in warm and hot dense plasmas over a wide range of electron densities and temperatures. By employing the ASD model, we investigate the photoionization process for the hydrogenlike carbon ion embedded in warm and hot dense plasmas with corresponding Coulomb coupling parameter ranges of 0.05 ≤ Γ ≤ 1.16, where Γ characterizes the ratio of the average potential to thermal energy. It is found that there are stronger plasma screening effects on the ionization energy and photoionization cross section due to the negative-energy electron distributions considered in the ASD model compared to those considering only free electrons. The present results from the ASD model show reasonable agreement with the classical Debye-Hückel (DH) model in weakly coupled plasmas. However, significant deviations of the ionization energy and cross section between these two models are observed in moderately and strongly coupled plasmas, due to the approximate treatment of the plasma-electron density distribution of the DH model. In the region of low photoelectron energies, the positions of the shape resonance peaks of the cross sections obtained from the ASD model differ significantly from those of the DH model due to the different screening effects.

Literature studies of the lattice Boltzmann method (LBM) demonstrate hydrodynamics beyond the continuum limit. This includes exact analytical solutions to the LBM, for the bulk velocity and shear stress of Couette flow under diffuse reflection at the walls through the solution of equivalent moment equations. We prove that the bulk velocity and shear stress of Couette flow with Maxwell-type boundary conditions at the walls, as specified by two-dimensional isothermal lattice Boltzmann models, are inherently linear in Mach number. Our finding enables a systematic variational approach to be formulated that exhibits superior computational efficiency than the previously reported moment method. Specifically, the number of partial differential equations (PDEs) in the variational method grows linearly with quadrature order while the number of moment method PDEs grows quadratically. The variational method directly yields a system of linear PDEs that provide exact analytical solutions to the LBM bulk velocity field and shear stress for Couette flow with Maxwell-type boundary conditions. It is anticipated that this variational approach will find utility in calculating analytical solutions for novel lattice Boltzmann quadrature schemes and other flows.

We study the spectral properties of the Dirichlet-to-Neumann operator and the related Steklov problem in spheroidal domains ranging from a needle to a disk. An explicit matrix representation of this operator for both interior and exterior problems is derived. We show how the anisotropy of spheroids affects the eigenvalues and eigenfunctions of the operator. As examples of physical applications, we discuss diffusion-controlled reactions on spheroidal partially reactive targets and the statistics of encounters between the diffusing particle and the spheroidal boundary.

In the early 2000s, Geniet and Leon [Phys. Rev. Lett. 89, 134102 (2002)] discovered the nonlinear supratransmission (NST) in a medium with a forbidden frequency band gap. It is a process in which nonlinear structures are created by a sinusoidal harmonic boundary condition imposed at a frequency in the band gap. The present study extends this concept and shows that an optimal shape of a periodic nonsinusoidal excitation may induce (or inhibit) energy flow through the lattice below (or above) the NST threshold, demonstrating that nonlinear supratransmission is reliant not only on the driving amplitude but also on its shape. This is evidenced through numerical simulations and mathematical calculations varying the excitation signal shape in the Fermi-Pasta-Ulam case study. Setting the shape parameter to zero recovers the results of the literature in relation to the sinusoidal signal.

In amorphous materials, plasticity is localized and occurs as shear transformations. It was recently shown by Wu et al. that these shear transformations can be predicted by applying topological defect concepts developed for liquid crystals to an analysis of vibrational eigenmodes [Z. W. Wu et al., Nat. Commun. 14, 2955 (2023)]. This study relates the $-1$ topological defects to the displacement fields expected of an Eshelby inclusion, which are characterized by an orientation and the magnitude of the eigenstrain. A corresponding orientation and magnitude can be defined for each defect using the local displacement field around each defect. These parameters characterize the plastic stress relaxation associated with the local structural rearrangement and can be extracted using the fit to either the global displacement field or the local field. Both methods provide a reasonable estimation of the molecular-dynamics-measured stress drop, confirming the localized nature of the displacements that control both long-range deformation and stress relaxation.

Impact crater experiments in granular media traditionally involve loosely packed sand targets. However, this study investigates granular impact craters on both loosely and more tightly packed sand targets. We report experiments that consistently adhere to power-law scaling laws for diameter as a function of impacting energy, similar to those reported by other groups for their experiments utilizing both solid and granular projectiles. In contrast, we observe significant deviations in the depth versus energy power law predicted by previous models. To address this discrepancy, we introduce a physical model of uniaxial compression that explains how depth saturates in granular collisions. Furthermore, we present an energy balance alongside this model that describes the energy transfer mechanisms acting during crater formation. We found a better way to transfer vertical momentum to horizontal degrees of freedom as the impact surface compacts, resulting in shallow craters on compacted sandbox targets. Our results reveal depth-to-diameter aspect ratios from approximately 0.051 to 0.094, allowing us to interpret the shallowness of planetary craters at the light of the uniaxial compression mechanism proposed in this work.

Statistical divergences are important tools in data analysis, information theory, and statistical physics, and there exist well-known inequalities on their bounds. However, in many circumstances involving temporal evolution, one needs limitations on the rates of such quantities instead. Here, several general upper bounds on the rates of some f-divergences are derived, valid for any type of stochastic dynamics (both Markovian and non-Markovian), in terms of information-like and/or thermodynamic observables. As special cases, the analytical bounds on the rate of mutual information are obtained. The major role in all those limitations is played by temporal Fisher information, characterizing the speed of global system dynamics, and some of them contain entropy production, suggesting a link with stochastic thermodynamics. Indeed, the derived inequalities can be used for estimation of minimal dissipation and global speed in thermodynamic stochastic systems. Specific applications of these inequalities in physics and neuroscience are given, which include the bounds on the rates of free energy and work in nonequilibrium systems, limits on the speed of information gain in learning synapses, as well as the bounds on the speed of predictive inference and learning rate. Overall, the derived bounds can be applied to any complex network of interacting elements, where predictability and thermodynamics of network dynamics are of prime concern.

We show that numerical linked cluster expansions (NLCEs) based on sufficiently large building blocks allow one to obtain accurate low-temperature results for the thermodynamic properties of spin lattice models with continuous disorder distributions. Specifically, we show that such results can be obtained computing the disorder averages in the NLCE clusters before calculating their weights. We provide a proof of concept using three different NLCEs based on $\mathsf{L}$, square, and rectangle building blocks. We consider both classical (Ising) and quantum (Heisenberg) spin-$\frac{1}{2}$ models and show that convergence can be achieved down to temperatures that are up to two orders of magnitude lower than the relevant energy scale in the model. Additionally, we provide evidence that in one dimension one can obtain accurate results for observables such as the energy down to their ground-state values.

A quantum random-walk model is established on a one-dimensional periodic lattice that fluctuates between two possible states. This model is defined by Lindblad rate equations that incorporate the transition rates between the two lattice states. Leveraging the system's symmetries, the particle velocity can be described using a finite set of equations, even though the state space is of infinite dimension. These equations yield an analytical expression for the velocity in the long-time limit, which is employed to analyze the characteristics of directed motion. Notably, the velocity can exhibit multiple inversions, and to achieve directed motion, distinct, nonzero transition rates between lattice states are required.

We discuss the two-dimensional motion of a Brownian particle that is confined to a harmonic trap and driven by a shear flow. The surrounding medium induces memory effects modeled by a linear, typically nonreciprocal coupling of the particle coordinates to an auxiliary (hidden) variable. The system's behavior resulting from the microscopic Langevin equations for the three variables is analyzed by means of exact moment equations derived from the Fokker-Planck representation, and numerical Brownian dynamics simulations. Increasing the shear rate beyond a critical value we observe, for suitable coupling scenarios with nonreciprocal elements, a transition from a stationary to a nonstationary state, corresponding to an escape from the trap. We analyze this behavior, analytically and numerically, in terms of the associated moments of the probability distribution, and from the perspective of nonequilibrium thermodynamics. Intriguingly, the entropy production rate remains finite when crossing the stability threshold.

Genomic regions can acquire heritable epigenetic states through unique histone modifications, which lead to stable gene expression patterns without altering the underlying DNA sequence. However, the relationship between chromatin conformational dynamics and epigenetic stability is poorly understood. In this paper, we propose kinetic models to investigate the dynamic fluctuations of histone modifications and the spatial interactions between nucleosomes. Our model explicitly incorporates the influence of chemical modifications on the structural stability of chromatin and the contribution of chromatin contacts to the cooperative nature of chemical reactions. Through stochastic simulations and analytical theory, we have discovered distinct steady-state outcomes in different kinetic regimes, resembling a dynamical phase transition. Importantly, we have validated that the emergence of this transition, which occurs on biologically relevant timescales, is robust against variations in model design and parameters. Our findings suggest that the viscoelastic properties of chromatin and the timescale at which it transitions from a gel-like to a liquidlike state significantly impact dynamic processes that occur along the one-dimensional DNA sequence.

Monolayers of growing bacteria, confined within channel geometries, exhibit self-organization into a highly aligned laminar state along the axis of the channel. Although this phenomenon has been observed in experiments and simulations under various boundary conditions, the underlying physical mechanism driving this alignment remains unclear. In this study, we conduct simulations of growing bacteria in two-dimensional channel geometries perturbed by fixed obstacles, either circular or arc shaped, placed at the channel's center. Our findings reveal that even sizable obstacles cause only short-ranged disruptions to the baseline laminar state. These disruptions arise from a competition between local planar anchoring and bulk laminar alignment. At smaller obstacle sizes, bulk alignment fully dominates, while at larger sizes planar anchoring induces increasing local disruptions. Furthermore, our analysis indicates that the resulting configurations of the bacterial system display a striking resemblance to the arrangement of hard-rod smectic liquid crystals around circular obstacles. This suggests that modeling hard-rod bacterial monolayers as smectic, rather than nematic, liquid crystals may yield successful outcomes. The insights gained from our study contribute to the expanding body of research on bacterial growth in channels. Our work provides perspectives on the stability of the laminar state and extends our understanding to encompass more intricate confinement schemes.

There have been some interesting recent advances in understanding the notion of mechanical disorder in structural glasses and the statistical mechanics of these systems' low-energy excitations. Here we contribute to these advances by studying a minimal model for structural glasses' elasticity in which the degree of mechanical disorder—as characterized by recently introduced dimensionless quantifiers—is readily tunable over a very large range. We comprehensively investigate a number of scaling laws observed for various macro, meso and microscopic elastic properties, and rationalize them using scaling arguments. Interestingly, we demonstrate that the model features the universal quartic glassy vibrational density of states as seen in many atomistic and molecular models of structural glasses formed by cooling a melt. The emergence of this universal glassy spectrum highlights the role of self-organization (toward mechanical equilibrium) in its formation, and elucidates why models featuring structural frustration alone do not feature the same universal glassy spectrum. Finally, we discuss relations to existing work in the context of strain stiffening of elastic networks and of low-energy excitations in structural glasses, in addition to future research directions.

The self-propulsion (translational instability) of a gas bubble in a liquid undergoing parametrically induced axisymmetric shape distortion due to being forced by a temporally sinusoidal, spatially constant acoustic field is investigated. Employing a model which accounts for the nonlinear coupling between the spherical oscillations, the axial translation and shape deformation of the bubble, the parametric excitement of two neighboring shape modes by the fundamental resonance, at the same driving frequency is studied. It is shown that provided pertinent driving pressure threshold values are exceeded, the respective shape modes are excited on different timescales. The growth of the shape mode on the faster timescale saturates giving rise to sustained constant amplitude oscillations, while the growth of the shape mode on the slower timescale is both modulated and unbounded. During the growth of the second shape mode, growing, oscillatory bubble translation is also observed.

Glass formers are, in general, classified as strong or fragile depending on whether their relaxation rates follow Arrhenius or super-Arrhenius temperature dependence. There are, however, notable exceptions, such as water, which exhibit a fragile-to-strong (FTS) transition and behave as fragile and strong, respectively, at high and low temperatures. In this work, the FTS transition is studied using a distinguishable-particle lattice model previously demonstrated to be capable of simulating both strong and fragile glasses [C.-S. Lee, M. Lulli, L.-H. Zhang, H.-Y. Deng, and C.-H. Lam, Phys. Rev. Lett. 125, 265703 (2020)]. Starting with a bimodal pair-interaction distribution appropriate for fragile glasses, we show that by narrowing down the energy dispersion in the low-energy component of the distribution, a FTS transition is observed. The transition occurs at a temperature at which the stretching exponent of the relaxation is minimized, in agreement with previous molecular dynamics simulations.