Standard lattice Boltzmann methods (LBMs) are based on a symmetric discretization of the phase space, which amounts to study the evolution of particle distribution functions (PDFs) in a reference frame at rest. This choice induces a number of limitations when the simulated flow speed gets closer to the sound speed, such as velocity-dependent transport coefficients. The latter issue is usually referred to as a Galilean invariance defect. To restore the Galilean invariance of LBMs, it was proposed to study the evolution of PDFs in a comoving reference frame by relying on asymmetric shifted lattices [N. Frapolli, S. S. Chikatamarla, and I. V. Karlin, Phys. Rev. Lett. 117, 010604 (2016)]. From the numerical viewpoint, this corresponds to overcoming the rather restrictive Courant-Friedrichs-Lewy conditions on standard LBMs and modeling compressible flows while keeping memory consumption and processing costs to a minimum (therefore using the standard first-neighbor stencils). In the present work systematic physical error evaluations and stability analyses are conducted for different discrete equilibrium distribution functions (EDFs) and collision models. Thanks to them, it is possible to (1) better understand the effect of this solution on both physics and stability, (2) assess its viability as a way to extend the validity range of LBMs, and (3) quantify the importance of the reference state as compared to other parameters such as the equilibrium state and equilibration path. The results clearly show that, in theory, the concept of shifted lattices allows the scheme to deal with arbitrarily high values of the nondimensional velocity. Furthermore, just like the zero-Mach flow for the standard stencils, it is observed that setting the shift velocity to the fluid velocity results in optimal physical and numerical properties. In addition, a detailed analysis of the obtained results shows that the properties of different collision models and EDFs remain unchanged under the shift of stencil. In other words, by introducing a velocity shift in the stencil, the optimal operating point, in terms of physics and numerics, will also be shifted by the same vector regardless of the EDF or collision model considered. Eventually, while limited to the D2Q9 stencil with the nine possible first-neighbor shifts, the present study and corresponding conclusions can be extended to other stencils and velocity shifts in a straightforward manner.

Marine species reproduce and compete while being advected by turbulent flows. It is largely unknown, both theoretically and experimentally, how population dynamics and genetics are changed by the presence of fluid flows. Discrete agent-based simulations in continuous space allow for accurate treatment of advection and number fluctuations, but can be computationally expensive for even modest organism densities. In this report, we propose an algorithm to overcome some of these challenges. We first provide a thorough validation of the algorithm in one and two dimensions without flow. Next, we focus on the case of weakly compressible flows in two dimensions. This models organisms such as phytoplankton living at a specific depth in the three-dimensional, incompressible ocean experiencing upwelling and/or downwelling events. We show that organisms born at sources in a two-dimensional time-independent flow experience an increase in fixation probability.

Power-law-type distributions are extensively found when studying the behavior of many complex systems. However, due to limitations in data acquisition, empirical datasets often only cover a narrow range of observation, making it difficult to establish power-law behavior unambiguously. In this work we present a statistical procedure to merge different datasets, with two different aims. First, we obtain a broader fitting range for the statistics of different experiments or observations of the same system. Second, we establish whether two or more different systems may belong to the same universality class. By means of maximum likelihood estimation, this methodology provides rigorous statistical information to discern whether power-law exponents characterizing different datasets can be considered equal among them or not. This procedure is applied to the Gutenberg-Richter law for earthquakes and for synthetic earthquakes (acoustic emission events) generated in the laboratory: labquakes. Different earthquake catalogs have been merged finding a Gutenberg-Richter law holding for more than eight orders of magnitude in seismic moment. The value of the exponent of the energy distribution of labquakes depends on the material used in the compression experiments. By means of the procedure proposed in this manuscript, we find that the Gutenberg-Richter law for earthquakes and charcoal labquakes can be characterized by the same power-law exponent, whereas Vycor labquakes exhibit a significantly different exponent.

We study the Caldeira-Leggett model where a quantum Brownian particle interacts with an environment or a bath consisting of a collection of harmonic oscillators in the path-integral formalism. Compared to the contours that the paths take in the conventional Schwinger-Keldysh formalism, the paths in our study are deformed in the complex time plane as suggested by the recent study by C. Aron, G. Biroli, and L. F. Cugliandolo [SciPost Phys. 4, 008 (2018)]. This is done to investigate the connection between the symmetry properties in the Schwinger-Keldysh action and the equilibrium or nonequilibrium nature of the dynamics in an open quantum system. We derive the influence functional explicitly in this setting, which captures the effect of the coupling to the bath. We show that in equilibrium the action and the influence functional are invariant under a set of transformations of path-integral variables. The fluctuation-dissipation relation is obtained as a consequence of this symmetry. When the system is driven by an external time-dependent protocol, the symmetry is broken. From the terms that break the symmetry, we derive a quantum Jarzynski-like equality for a quantum mechanical worklike quantity given as a function of fluctuating quantum trajectory. In the classical limit, the transformations becomes those used in the functional integral formalism of the classical stochastic thermodynamics to derive the classical fluctuation theorem.

A family of super-regular (SR) breather solutions in systems with self-steepening effect and in the case of either normal or anomalous dispersion is derived analytically. Derivation is based on the Darboux transformation with a quadratic spectral parameter. In contrast to the SR breather solutions in $t$-symmetric systems such as the nonlinear Schrödinger equation, the new breathers found in the present work evolve asymmetrically even if started from symmetric initial conditions. The initial stage of this process is modulation instability. Numerical simulations confirm the excitation of the SR breathers when started from the approximate initial conditions leading at first to modulation instability. Our results offer the possibility of experimental observations of SR breather dynamics in systems with self-steepening effects, such as optical frequency-doubling crystals or magnetized plasmas.

The breaking of a plane self-excited dust-acoustic soliton in a dust cloud formed in stratified dc glow discharge plasma is studied. Both macroscopic and kinetic parameters of the dust component near the soliton are experimentally obtained. It is shown that the breaking of a soliton can accelerate charged particles to supersonic speeds. The theoretical interpretation of the experimental results is performed in the framework of the hydrodynamic plasma approach, as well as the single-particle approximation. Both dissipative and nondissipative cases are considered.

Much work has studied effective interactions between micron-sized particles carrying linkers forming reversible, interparticle linkages. These studies allowed understanding the equilibrium properties of colloids interacting through ligand-receptor interactions. Nevertheless, understanding the kinetics of multivalent interactions remains an open problem. Here, we study how molecular details of the linkers, such as the reaction rates at which interparticle linkages form or break, affect the relative dynamics of pairs of cross-linked colloids. Using a simulation method tracking single binding and unbinding events between complementary linkers, we rationalize recent experiments and prove that particles' interfaces can move across each other while being cross-linked. We clarify how, starting from diffusing colloids, the dynamics become arrested when increasing the number of interparticle linkages or decreasing the reaction rates. Before getting arrested, particles diffuse through rolling motion. The ability to detect rolling motion will be useful to shed new light on host-pathogen interactions.

The minimum vertex cover (Min-VC) problem is a well-known NP-hard problem. Earlier studies illustrate that the problem defined over the Erdös-Rényi random graph with a mean degree $c$ exhibits computational difficulty in searching the Min-VC set above a critical point $c=e=2.718...$. Here, we address how this difficulty is influenced by the mesoscopic structures of graphs. For this, we evaluate the critical condition of difficulty for the stochastic block model. We perform a detailed examination of the specific cases of two equal-size communities characterized by in and out degrees, which are denoted by $c_{\mathrm{in}}$ and $c_{\mathrm{out}}$, respectively. Our analysis based on the cavity method indicates that the solution search once becomes difficult when $c_{\mathrm{in}}+c_{\mathrm{out}}$ exceeds $e$ from below, but becomes easy again when $c_{\text{out}}$ is sufficiently larger than $c_{\mathrm{in}}$ in the region $c_{\mathrm{out}}>e$. Experiments based on various search algorithms support the theoretical prediction.

We study the non-Markovian Brownian motion of an electrically charged harmonic oscillator through the action of both a constant magnetic field and time-dependent force fields. The generalized Langevin equation with a friction memory kernel is used to derive the generalized phase-space Fokker-Planck equation for the harmonic oscillator in the absence and in the presence of time-dependent force fields. To achieve our goal, the characteristic function method is applied to obtain, in an accurate way, the theoretical description of the problem. We explicitly calculate the correlation and cross-correlation functions for the position and velocity vectors. We show that the relevant physics behind the theory is contained in the generalized diffusion coefficient, which accounts for the natural coupling between both the harmonic oscillator and magnetic field. Our theoretical results are compared with those previously reported in the literature.

In amorphous solids at finite temperatures the particles follow chaotic trajectories which, at temperatures sufficiently lower than the glass transition, are trapped in “cages.” Averaging their positions for times shorter than the diffusion time, one can define a time-averaged configuration. Under strain or stress, these average configurations undergo sharp plastic instabilities. In athermal glasses the understanding of plastic instabilities is furnished by the Hessian matrix and its eigenvalues and eigenfunctions. Here we propose an uplifting of Hessian methods to thermal glasses, with the aim of understanding the plastic responses in the time-averaged configuration. We discuss a number of potential definitions of Hessians and identify which of these can provide eigenvalues and eigenfunctions which can explain and predict the instabilities of the time-averaged configurations. The conclusion is that the nonaffine changes in the average configurations during an instability is accurately predicted by the eigenfunctions of the low-lying eigenvalues of the chosen Hessian.

The energy and spin diffusion behaviors in the one-dimensional classical Heisenberg spin chain have been systematically investigated using the equilibrium diffusion method. The spatiotemporal autocorrelation functions for energy and spin are calculated at finite and infinite temperatures. As conserved quantities, the spreading of excess energy and spin can be used to determine their actual diffusion behaviors. At low temperatures, the energy diffusion shows almost ballistic behavior, and spin shows superdiffusion behavior for finite chain size. For energy diffusion, normal diffusion behavior can be obtained when the temperature is higher than 0.75. For spin diffusion, normal diffusion behavior is observed at infinite temperature.

We have studied the distribution of traffic flow $q$ for the Nagel-Schreckenberg model by computer simulations. We applied a large-deviation approach, which allowed us to obtain the distribution $P\left(q\right)$ over more than one hundred decades in probability, down to probabilities like ${10}^{-140}$. This allowed us to characterize the flow distribution over a large range of the support and identify the characteristics of rare and even very rare traffic situations. We observe a change of the distribution shape when increasing the density of cars from the free flow to the congestion phase. Furthermore, we characterize typical and rare traffic situations by measuring correlations of $q$ to other quantities like density of standing cars or number and size of traffic jams.

Social discrimination seems to be a persistent phenomenon in many cultures. It is important to understand the mechanisms that lead people to judge others by the group to which they belong rather than individual qualities. It was recently shown that evolutionary (imitation) dynamics can lead to a hierarchical discrimination between agents marked with observable, but otherwise meaningless, labels. These findings suggest that it can give useful insight to describe the phenomenon of social discrimination in terms of spontaneous symmetry breaking. The investigations so far have, however, only considered binary labels. In this contribution we extend the investigations to models with up to seven different labels. We find the features known from the binary label model remain remarkably robust when the number of labels is increased. We also discover a new feature, namely that it is more likely for neighbors to have strategies which are similar, in the sense that they agree on how to act toward a subset of the labels.

While the timing of intracellular events is essential for many cellular processes, gene expression inside a single cell can exhibit substantial cell-to-cell variability, raising the question of how cells ensure precision in event timing despite such stochasticity. We address this question by analyzing a biologically reasonable model of gene expression in the context of first passage time (FPT), focusing on two experimentally measurable statistics: mean FPT (MFPT) and timing variability (TV). We show that (1) transcriptional burst size (BS) and burst frequency (BF) can minimize the TV; (2) translational BS monotonically reduces the MFPT to a nonzero low bound; (3) the timescale of promoter kinetics can minimize both the MFPT and the TV, depending on the ratio of the on-switching rate over the off-switching rate; and (4) positive feedback regulation of any form can all minimize the TV, whereas negative feedback regulation of transcriptional BF or BS always enhances the TV. These control strategies can have broad implications for diverse cellular processes relying on precise temporal triggering of events.

We present the analysis of the impedance spectra for a binary electrolyte confined between blocking electrodes with dielectric layers. An expression for the impedance is derived from Poisson-Nernst-Planck equations in the linear approximation taking into account the voltage drop on the dielectric layer. The analysis shows that characteristic features of the frequency dependence of the impedance are determined by the ratio of the Debye length and the effective thickness of the dielectric layer. The impact of the dielectric layer is especially strong in the case of high concentrated electrolytes, where the Debye length is small and thus comparable to the effective thickness of the dielectric layer. To verify the model, measurements of the impedance spectra and transient currents in a liquid crystal 4-n-pentyl-$4^{\prime }$-cyanobiphenyl (5CB) confined between polymer-coated electrodes in cells of different thicknesses are performed. The estimates for the diffusion coefficient and ion concentration in 5CB obtained from the analysis of the impedance spectra and the transient currents are consistent and agree with previously reported data. We demonstrate that calculations of the ion parameters from the impedance spectra without taking into account the dielectric layer contribution lead in most cases to incorrect results. Application of the model to analyze violations of the low-frequency impedance scaling and contradictions in the estimates of the ion parameters recently found in some ionic electrolytes are discussed.

We find an instability resulting in generation of large-scale vorticity in a fast-rotating small-scale turbulence or turbulent convection with inhomogeneous fluid density along the rotational axis in anelastic approximation. The large-scale instability causes excitation of two modes: (i) the mode with dominant vertical vorticity and with the mean velocity being independent of the vertical coordinate; (ii) the mode with dominant horizontal vorticity and with the mean momentum being independent of the vertical coordinate. The mode with the dominant vertical vorticity can be excited in a fast-rotating density-stratified hydrodynamic turbulence or turbulent convection. For this mode, the mean entropy is depleted inside the cyclonic vortices, while it is enhanced inside the anticyclonic vortices. The mode with the dominant horizontal vorticity can be excited only in a fast-rotating density-stratified turbulent convection. The developed theory may be relevant for explanation of an origin of large spots observed as immense storms in great planets, e.g., the Great Red Spot in Jupiter and large spots in Saturn. It may be also useful for explanation of an origin of high-latitude spots in rapidly rotating late-type stars.

A kagome-white-eye-honeycomb hexagonal superlattice pattern is observed in dielectric barrier discharge and studied by two photomultipliers, an intensified charge-coupled device, and a spectrograph. It consists of four sublattices, which are kagome structure, halos of white-eye hexagon, honeycomb frame, and central spots hexagon. As an unconventional symmetry sublattice, the kagome structure results in amazing sublattice discharge order. Based on the consideration that any discharge is related to all previous discharges, a mapping matrix which governs self-organization of the pattern is proposed. These results will promote the development of the pattern dynamics.

The nonstationary diffusion of particles in a medium with static random traps or sinks is considered. The question of the self-averaging of the diffusion coefficient (or, equivalently, of the mean-square displacement) is addressed for the fluctuation regime in the long-time limit. The property of self-averaging is needed for the result of a single measurement to be representative and reproducible. It is demonstrated that the diffusion coefficient of the surviving particles is a strongly non-self-averaging quantity: In a $d$-dimensional system its reciprocal standard deviation grows with time exponentially $\approx exp\left[{\mathrm{const}}_{d,1}{t}^{d/\left(d+2\right)}\right]$. The same result is reproduced in the “normalized” formulation “per one survivor on average.” The case when all the particles, both the survivors and the trapped ones, are contributing to the diffusion coefficient and its variance is considered also. Non-self-averaging is demonstrated for this case as well, the fluctuations of the diffusion coefficient being of the same order as its average value. The critical dimension, above which the mean-field result becomes exact, is infinite—due to the drastic difference between the classes of trajectories, upon which the corresponding results are being built. In high dimensions the strong non-self-averaging of survivors is preserved. For the case of all the particles taken into account, the nonstrong non-self-averaging is retained for any finite dimension. However, for $d\to \infty$ the limiting value of the reciprocal standard deviation, calculated for all the particles, decreases to zero. This signifies restoration of the self-averaging in some sense. In all the cases, the time evolution of the average characteristics and of their variances is governed by the decaying concentration of the survivors in fluctuational cavities.

We investigate ionic criticality on the basis of a specially devised spherical model that accounts both for Coulomb and nonionic forces in binary systems. We show in detail here the consequences of the entanglement of density and charge correlation functions $G_{NN}$ and $G_{ZZ}$ on criticality and screening. We also show on this soluble model how, because of electroneutrality, the long-range Coulomb interactions do not change the universality class of criticality in the model driven primarily by sufficiently attractive nonionic interactions. Near criticality, $G_{NN}$ and $G_{ZZ}$ are fully decoupled in charge symmetric systems. However, in more realistic nonsymmetric models, charge and density fluctuations couple in leading order so that the charge and density correlation lengths diverge asymptotically in a similar way. Similarly, the Stillinger-Lovett sum rule, which characterizes a conducting fluid, is violated at criticality in nonsymmetric models when the critical-point density-decay exponent $\eta$ vanishes. In addition, if quantum effects are accounted for semiclassically by incorporating algebraically decaying interactions, $G_{ZZ}$ decays only as a power law in the whole phase space, contrary to the usually expected exponential Debye screening. We expect these results on this soluble toy model to be general and to reveal general mechanisms ruling ionic criticality.

The critical point is a fixed point in finite-size scaling. To quantify the behavior of such a fixed point, we define, at a given temperature and scaling exponent ratio, the width of scaled observables for different sizes. The minimum of the width reveals the position of the fixed point, its corresponding phase transition temperature, and scaling exponent ratio. The value of this ratio tells the nature of the fixed point, which can be a critical point, a point of the first-order phase transition line, or a point of the crossover region. To demonstrate the effectiveness of this method, we apply it to three typical samples produced by the three-dimensional three-state Potts model. Results show the method to be more precise and effective than conventional methods. Finally, we discuss a possible application at the Beam Energy Scan program of the Relativistic Heavy Ion Collider.