The prevalence of wealth inequality propels us to characterize its origin and progression via empirical and theoretical studies. The yard-sale (YS) model, in which a portion of the smaller wealth is transferred between two individuals, culminates in the concentration of almost all wealth to a single individual, while distributing the rest of the wealth with a power law of exponent one. By incorporating redistribution to the model, in which the transferred wealth is proportional to the sender's wealth, we show that such extreme inequality is suppressed if the frequency ratio of redistribution to the YS-type exchange exceeds the inverse of the population size. Studying our model on a sparsely-connected population, we find that the wealth inequality ceases to grow for a period, when local rich nodes can no longer acquire wealth from their broke nearest neighbors. Subsequently, inequality resumes growth due to the redistribution effect by allowing locally amassed wealth to move and coalesce. Analyzing the Langevin equations and the coalescing random walk on complex networks, we elucidate the scaling behaviors of wealth inequality in those multiple phases. These findings reveal the influence of network structure on wealth distribution, offering a novel perspective on wealth inequality.

Combinatorial mechanical metamaterials feature spatially textured soft modes that yield exotic and useful mechanical properties. While a single soft mode often can be rationally designed by following a set of tiling rules for the building blocks of the metamaterial, it is an open question what design rules are required to realize multiple soft modes. Multimodal metamaterials would allow for advanced mechanical functionalities that can be selected on the fly. Here we introduce a transfer matrix-like framework to design multiple soft modes in combinatorial metamaterials composed of aperiodic tilings of building blocks. We use this framework to derive rules for multimodal designs for a specific family of building blocks. We show that such designs require a large number of degeneracies between constraints, and find precise rules on the real space configuration that allow such degeneracies. These rules are significantly more complex than the simple tiling rules that emerge for single-mode metamaterials. For the specific example studied here, they can be expressed as local rules for tiles composed of pairs of building blocks in combination with a nonlocal rule in the form of a global constraint on the type of tiles that are allowed to appear together anywhere in the configuration. This nonlocal rule is exclusive to multimodal metamaterials and exemplifies the complexity of rational design of multimode metamaterials. Our framework is a first step towards a systematic design strategy of multimodal metamaterials with spatially textured soft modes.

Orthotropic shell structures are ubiquitous in biology and engineering, from bacterial cell walls to reinforced domes. We present a rescaling transformation that maps an orthotropic shallow shell to an isotropic one with a different local geometry. The mapping is applicable to any shell section for which the material orthotropy directions match the principal curvature directions, assuming the commonly used Huber form for the orthotropic shear modulus. Using the rescaling transformation, we derive exact expressions for the buckling pressure as well as the linear indentation response of orthotropic cylinders and general ellipsoids of revolution, which we verify against numerical simulations. Our analysis disentangles the separate contributions of geometric and material anisotropy to shell rigidity. In particular, we identify the geometric mean of orthotropic elastic constants as the key quantifier of material stiffness, playing a role akin to the Gaussian curvature which captures the geometric stiffness contribution. Besides providing insights into the mechanical response of orthotropic shells, our work rigorously establishes the validity of isotropic approximations to orthotropic shells and also identifies situations in which these approximations might fail.

Heterogeneity is the concept we encounter in numerous research areas and everyday life. While “not mixing apples and oranges” is easy to grasp, a more quantitative approach to such segregation is not always readily available. Consider the problem from a different angle: To what extent does one have to make apples more orange and oranges more “apple-shaped” to put them into the same basket (according to their appearance alone)? This question highlights the central problem of the blurred interface between heterogeneous and homogeneous, which also depends on the metrics used for its identification. This work uncovers the physics of structural stationarity quantification, based on correlation functions (CFs) and clustering based on CFs different between image subregions. By applying the methodology to a wide variety of synthetic and real images of binary porous media, we confirmed computationally that only periodically unit-celled structures and images produced by stationary processes with resolutions close to infinity are strictly stationary. Natural structures without recurring unit cells are only weakly stationary. We established a physically meaningful definition for these stationarity types and their distinction from nonstationarity. In addition, the importance of information content of the chosen metrics is highlighted and discussed. We believe the methodology as proposed in this contribution will find its way into numerous research areas dealing with materials, structures, and measurements and modeling based on structural imaging information.

We construct a simple field theory in which a sphaleron, i.e., a saddle-point particle-like solution, forms a semi-BPS state with a background defect that is an impurity. This means that there is no static force between the sphaleron and the impurity. Therefore, such a sphaleron-impurity system is very much like usual BPS multisolitons, however, still possessing an unstable direction allowing for its decay. We study dynamics of the sphaleron in such a system.

We use linear stability analysis and hybrid lattice Boltzmann simulations to study the dynamical behavior of an active nematic confined in a channel made of viscoelastic material. We find that the quiescent, ordered active nematic is unstable above a critical activity. The transition is to a steady flow state for high elasticity of the channel surroundings. However, below a threshold elastic modulus, the system produces spontaneous oscillations with periodic flow reversals. We provide a phase diagram that highlights the region where time-periodic oscillations are observed and explain how they are produced by the interplay of activity and viscoelasticity. Our results suggest experiments to study the role of viscoelastic confinement in the spatiotemporal organization and control of active matter.

We propose generalized equilibria of a three-dimensional color-gradient lattice Boltzmann model for two-component two-phase flows using higher-order Hermite polynomials. Although the resulting equilibrium distribution function, which includes a sixth-order term on the velocity, is computationally cumbersome, its equilibrium central moments (CMs) are velocity-independent and have a simplified form. Numerical experiments show that our approach, as in Wen et al. [Phys. Rev. E 100, 023301 (2019)] who consider terms up to third order, improves the Galilean invariance compared to that of the conventional approach. Dynamic problems can be solved with high accuracy at a density ratio of 10; however, the accuracy is still limited to a density ratio of 1000. For lower density ratios, the generalized equilibria benefit from the CM-based multiple-relaxation-time model, especially at very high Reynolds numbers, significantly improving the numerical stability.

We study time-reversal symmetry breaking in non-Hermitian fluctuating field theories with conserved dynamics, comprising the mesoscopic descriptions of a wide range of nonequilibrium phenomena. They exhibit continuous parity-time ($\mathcal{PT}$) symmetry-breaking phase transitions to dynamical phases. For two concrete transition scenarios, exclusive to non-Hermitian dynamics, namely, oscillatory instabilities and critical exceptional points, a low-noise expansion exposes a pretransitional surge of the mesoscale (informatic) entropy production rate, inside the static phases. Its scaling in the susceptibility contrasts conventional critical points (such as second-order phase transitions), where the susceptibility also diverges, but the entropy production generally remains finite. The difference can be attributed to active fluctuations in the wavelengths that become unstable. For critical exceptional points, we identify the coupling of eigenmodes as the entropy-generating mechanism, causing a drastic noise amplification in the Goldstone mode.

We consider a system of globally coupled phase-only oscillators with distributed intrinsic frequencies and evolving in the presence of distributed Gaussian white noise, namely, a Gaussian white noise whose strength for every oscillator is a specified function of its intrinsic frequency. In the absence of noise, the model reduces to the celebrated Kuramoto model of spontaneous synchronization. For two specific forms of the mentioned functional dependence and for a symmetric and unimodal distribution of the intrinsic frequencies, we unveil the rich long-time behavior that the system exhibits, which stands in stark contrast to the case in which the noise strength is the same for all the oscillators, namely, in the studied dynamics, the system may exist in either a synchronized, or an incoherent, or a time-periodic state; interestingly, all these states also appear as long-time solutions of the Kuramoto dynamics for the case of bimodal frequency distributions, but in the absence of any noise in the dynamics.

We investigate the dynamics of a quantum system subjected to a time-dependent and conditional resetting protocol. Namely, we ask what happens when the unitary evolution of the system is repeatedly interrupted at random time instants with an instantaneous reset to a specified set of reset configurations taking place with a probability that depends on the current configuration of the system at the instant of reset? Analyzing the protocol in the framework of the so-called tight-binding model describing the hopping of a quantum particle to nearest-neighbor sites in a one-dimensional open lattice, we obtain analytical results for the probability of finding the particle on the different sites of the lattice. We explore a variety of dynamical scenarios, including the one in which the resetting time intervals are sampled from an exponential as well as from a power-law distribution, and a setup that includes a Floquet-type Hamiltonian involving an external periodic forcing. Under exponential resetting, and in both the presence and absence of the external forcing, the system relaxes to a stationary state characterized by localization of the particle around the reset sites. The choice of the reset sites plays a defining role in dictating the relative probability of finding the particle at the reset sites as well as in determining the overall spatial profile of the site-occupation probability. Indeed, a simple choice can be engineered that makes the spatial profile highly asymmetric even when the bare dynamics does not involve the effect of any bias. Furthermore, analyzing the case of power-law resetting serves to demonstrate that the attainment of the stationary state in this quantum problem is not always evident and depends crucially on whether the distribution of reset time intervals has a finite or an infinite mean.

The present work is concerned with the uncertainty propagation of the wave turbulent system. In particular, we study the temporal development and long-term behavior of the probability with respect to the amplitude and phase of complex-valued waves constituting the generic four-wave system of turbulence. Our approach to approximating the target distribution function is via the three steps: (i) to grasp the physical process described by the true turbulence model as random process, (ii) to determine the stochastic differential equation whose solution exhibits statistically similar behavior with the underlying turbulent signal, and (iii) to solve the corresponding Kolmogorov forward equation. Our implementation of the methodology is distinguished by employing a number of simplified stochastic models and applying one of them in the adaptive fashion which varies subject to the different parameter regime of the true dynamical system model. Accordingly, we become able to demonstrate the effectiveness of this reduced-order modeling framework for the analysis of the turbulent system characterized by not only weak but strong interactions among the nonlinear waves. We numerically corroborate our theoretical predictions in the context of the generalized Majda-Mclaughlin-Tabak wave turbulence prototype.

We consider a generic four-stroke quantum Otto engine consisting of two unitary and two thermalization strokes with an arbitrary many-body working medium. Using the Schwinger-Keldysh nonequilibrium Green's function formalism, we provide an analytical expression for the cumulant generating function corresponding to the joint probability distribution of nonequilibrium work and heat. The obtained result is valid up to the second order of the external driving amplitude. We then focus on the linear response limit and obtained Onsager's transport coefficients for the generic Otto cycle and show that the traditional fluctuation-dissipation relation for the total work is violated in the quantum domain, whereas for heat it is preserved. This leads to remarkable consequences in obtaining universal constraints on heat and work fluctuations for engine and refrigerator regimes of the Otto cycle and further allows us to make connections to the thermodynamic uncertainty relations. These findings are illustrated using a paradigmatic model that can be feasibly implemented in experiments.

We study the universal nonstationary evolution of wave turbulence (WT) in Bose-Einstein condensates (BECs). Their temporal evolution can exhibit different kinds of self-similar behavior corresponding to a large-time asymptotic of the system or to a finite-time blowup. We identify self-similar regimes in BECs by numerically simulating the forced and unforced Gross-Pitaevskii equation (GPE) and the associated wave kinetic equation (WKE) for the direct and inverse cascades, respectively. In both the GPE and the WKE simulations for the direct cascade, we observe the first-kind self-similarity that is fully determined by energy conservation. For the inverse cascade evolution, we verify the existence of a self-similar evolution of the second kind describing self-accelerating dynamics of the spectrum leading to blowup at the zero mode (condensate) at a finite time. We believe that the universal self-similar spectra found in the present paper are as important and relevant for understanding the BEC turbulence in past and future experiments as the commonly studied stationary Kolmogorov-Zakharov (KZ) spectra.

Different types of cells exhibit a universal power-law rheology, but the mechanism underneath is still unclear. Based on the exponential distribution of actin filament length, we treat the cell cortex as a collection of chains of crosslinkers with exponentially distributed binding energy, and show that the power-law exponent of its stress relaxation should scale with the chain length. Through this model, we are able to explain how the exponent can be regulated by the crosslinker number and imposed strain during cortex relaxation. Network statistics show that the average length of filament-crosslinker chains decreases with the crosslinker number, which endows a denser network with lower exponent. Due to gradual molecular alignment with the stretch direction, the number of effectively stretched crosslinkers in the network is found to increase with the imposed strain. This effective growth in network density diminishes the exponent under large strain. By incorporating the inclined angle of crosslinkers into the model without in-series structure, we show that the exponent cannot be altered by crosslinker rotation directly, refining our previous conjectures. This work may help to understand cellular mechanics from the molecular perspective.

Using the recently introduced parsimonious Metropolis Monte Carlo algorithm, bead-stick polymers both with infinite-range Lennard-Jones interaction and with truncation are simulated. The focus lies on determining the Boyle temperature for long chains with thousands of repeat units and on testing for theoretically predicted logarithmic corrections. Subsequently the behavior at the infinite-chain transition temperature, i.e., the $\mathrm{\Theta }$ temperature, is studied for chains with up to $N=32768$ repeat units by investigation of the scaling of the end-to-end distance, the radius of gyration, the specific heat, and their derivatives with $N$.

We study the nonreciprocal Cahn-Hilliard model with thermal noise as a prototypical example of a generic class of non-Hermitian stochastic field theories, analyzed in two companion papers [Suchanek, Kroy, and Loos, Phys. Rev. Lett. 131, 258302 (2023); Phys. Rev. E 108, 064123 (2023)]. Due to the nonreciprocal coupling between two field components, the model is inherently out of equilibrium and can be regarded as an active field theory. Beyond the conventional homogeneous and static-demixed phases, it exhibits a traveling-wave phase, which can be entered via either an oscillatory instability or a critical exceptional point. By means of a Fourier decomposition of the entropy production rate, we quantify the associated scale-resolved time-reversal symmetry breaking, in all phases and across the transitions, in the low-noise regime. Our perturbative calculation reveals its dependence on the strength of the nonreciprocal coupling. Surging entropy production near the static-dynamic transitions can be attributed to entropy-generating fluctuations in the longest wavelength Fourier mode and heralds the emerging traveling wave. Its translational dynamics can be mapped on the dissipative ballistic motion of an active (quasi)particle.

When objects are forced to flow through constrictions their transport can be frustrated temporarily or permanently due to the formation of arches in the region of the bottleneck. While such systems have been intensively studied in the case of solid particles in a gas phase being forced by gravitational forces, the case of solid particles suspended in a liquid phase, forced by the liquid itself, has received much less attention. In this case, the influence of the liquid flow on the transport efficiency is not well understood yet, leading to several apparently trivial but yet unanswered questions, e.g., would an increase of the liquid flow improve the transport of particles or worsen it? Although some experimental data are already available, they lack enough detail to give a complete answer to such a question. Numerical models would be needed to scrutinize the system deeper. In this paper, we study this system making use of an advanced discrete particle solver (mercurydpm) and an approximated numerical model for the liquid drag and compare the results with experimental data.

The compression of brittle porous media can lead to the propagation of compaction bands. Although such localization phenomena have been observed in different geometries, including cuboidal and axisymmetric uniaxial compression, the role of boundary geometry on compaction features has yet to be explored, despite its relevance in geological conditions and industrial processes. To this end, we investigate the influence of shaped boundaries and inhomogeneous inclusions in a model brittle material made of puffed rice cereal. Using a variety of geometries, we show that compaction bands assume the shape of nearby boundaries, but return to a default planar form a distance away from them. Remarkably, the band aligns parallel to characteristic lines of minor principal stress obtained from a simple linear elastic model. The compelling correlation between the rotation of the principal stress directions and compaction band orientation holds implications for the geological interpretation of localized patterns in rocks and for comprehending the formation of weak planes in pharmaceutical tablets.

Coupled limit cycle oscillators with pairwise interactions are known to depict phase transitions from an oscillatory state to amplitude or oscillation death. This Research Letter introduces a scheme to incorporate higher-order interactions which cannot be decomposed into pairwise interactions and investigates the dynamical evolution of Stuart-Landau oscillators under the impression of such a coupling. We discover an oscillator death state through a first-order (explosive) phase transition in which a single, coupling-dependent stable death state away from the origin exists in isolation without being accompanied by any other stable state usually existing for pairwise couplings. We call such a state a solitary death state. Contrary to widespread subcritical Hopf bifurcation, here we report homoclinic bifurcation as an origin of the explosive death state. Moreover, this explosive transition to the death state is preceded by a surge in amplitude and followed by a revival of the oscillations. The analytical value of the critical coupling strength for the solitary death state agrees with the simulation results. Finally, we point out the resemblance of the results with different dynamical states associated with epileptic seizures.