We study very simple sorting algorithms based on a probabilistic comparator model. In this model, errors in comparing two elements are due to (1) the energy or effort put in the comparison and (2) the difference between the compared elements. Such algorithms repeatedly compare and swap pairs of randomly chosen elements, and they correspond to natural Markovian processes. The study of these Markov chains reveals an interesting phenomenon. Namely, in several cases, the algorithm that repeatedly compares only adjacent elements is better than the one making arbitrary comparisons: in the long-run, the former algorithm produces sequences that are “better sorted”. The analysis of the underlying Markov chain poses interesting questions as the latter algorithm yields a nonreversible chain, and therefore its stationary distribution seems difficult to calculate explicitly. We nevertheless provide bounds on the stationary distributions and on the mixing time of these processes in several restrictions.

Coupling two or more self-oscillating systems may stabilize their zero-amplitude rest state, therefore quenching their oscillation. This phenomenon is termed “amplitude death.” Well known and studied in classical self-oscillators, amplitude death was only recently investigated in quantum self-oscillators [Ishibashi and Kanamoto, Phys. Rev. E 96, 052210 (2017)]. Quantitative differences between the classical and quantum descriptions were found. Here, we demonstrate that for quantum self-oscillators with anharmonicity in their energy spectrum, multiple resonances in the mean phonon number can be observed. This is a result of the discrete energy spectrum of these oscillators, and is not present in the corresponding classical model. Experiments can be realized with current technology and would demonstrate these genuine quantum effects in the amplitude death phenomenon.

A self-consistent kinetic theory of Thomson scattering of an electromagnetic field by a nonuniform plasma is derived. We show that not only the imaginary part, but also the time and space derivatives of the real part of the dielectric susceptibility determine the amplitude and the width of the Thomson scattering spectral lines. As a result of inhomogeneity, these properties become asymmetric with respect to inversion of the sign of the frequency. Our theory provides a method of a remote probing and measurement of electron density gradients in plasma; this is based on the demonstrated asymmetry of the Thomson scattering lines.

Modeling high-temperature (tens or hundreds of eV), dense plasmas is challenging due to the multitude of non-negligible physical effects including significant partial ionization and multisite effects. These effects cause the breakdown or intractability of common methods and approximations used at low temperatures, such as pseudopotentials or plane-wave basis sets. Here we explore the Korringa-Kohn-Rostoker Green's function method at these high-temperature conditions. The method is all electron, does not rely on pseudopotentials, and uses a spherical harmonic basis set, and so avoids the aforementioned limitations. It is found to be accurate for solid density aluminum and iron plasmas when compared to a plane-wave method at low temperature, while being able to access high temperatures.

Theoretical understanding of how a deep neural network (DNN) extracts features from input images is still unclear, but it is widely believed that the extraction is performed hierarchically through a process of coarse graining. It reminds us of the basic renormalization group (RG) concept in statistical physics. In order to explore possible relations between DNN and RG, we use the restricted Boltzmann machine (RBM) applied to an Ising model and construct a flow of model parameters (in particular, temperature) generated by the RBM. We show that the unsupervised RBM trained by spin configurations at various temperatures from $T=0$ to $T=6$ generates a flow along which the temperature approaches the critical value $T_c=2.27$. This behavior is the opposite of the typical RG flow of the Ising model. By analyzing various properties of the weight matrices of the trained RBM, we discuss why it flows towards $T_c$ and how the RBM learns to extract features of spin configurations.

A recent approach that led to compact frequency domain formulations of the cubic and quartic fluctuation-dissipation theorems (FDTs) for the classical one-component plasma (OCP) [Golden and Heath, J. Stat. Phys. 162, 199 (2016)] is generalized to accommodate binary ionic mixtures. Paralleling the procedure followed for the OCP, the basic premise underlying the present approach is that a $(\mathbf{k},\omega )$ 4-vector rotational symmetry, known to be a pivotal feature in the frequency domain architectures of the linear and quadratic fluctuation-dissipation relations for a variety of Coulomb plasmas [Golden et al., J. Stat. Phys. 6, 87 (1972); J. Stat. Phys. 29, 281 (1982); Golden, Phys. Rev. E 59, 228 (1999)], is expected to be a pivotal feature of the frequency domain architectures of the higher-order members of the FDT hierarchy. On this premise, each member, in its most tractable form, connects a single ($p+1$)-point dynamical structure function to a linear combination of ($p+1$)-order $p$ density response functions; by definition, such a combination must also remain invariant under rotation of their $({\mathbf{k}}_1,{\omega }_1),({\mathbf{k}}_2,{\omega }_2),...,({\mathbf{k}}_p,{\omega }_p)$, $({\mathbf{k}}_1+{\mathbf{k}}_2+\cdots +{\mathbf{k}}_p,{\omega }_1+{\omega }_2+\cdots +{\omega }_p)$ 4-vector arguments. Assigned to each 4-vector is a species index that corotates in lock step. Consistency is assured by matching the static limits of the resulting frequency domain cubic and quartic FDTs to their exact static counterparts independently derived in the present work via a conventional time-independent perturbation expansion of the Liouville distribution function in its macrocanonical form. The proposed procedure entirely circumvents the daunting issues of entangled Liouville space paths and nested Poisson brackets that one would encounter if one attempted to use the conventional time-dependent perturbation-theoretic Kubo approach to establish the frequency domain FDTs beyond quadratic order.

We study the ordering statistics of four random walkers on the line, obtaining a much improved estimate for the long-time decay exponent of the probability that a particle leads to time $t$, $P_{\mathrm{lead}}\left(t\right)\sim t^{-0.91287850}$, and that a particle lags to time $t$ (never assumes the lead), $P_{\mathrm{lag}}\left(t\right)\sim t^{-0.30763604}$. Exponents of several other ordering statistics for $N=4$ walkers are obtained to eight-digit accuracy as well. The subtle correlations between $n$ walkers that lag jointly, out of a field of $N$, are discussed: for $N=3$ there are no correlations and $P_{\mathrm{lead}}\left(t\right)\sim P_{\mathrm{lag}}{\left(t\right)}^2$. In contrast, our results rule out the possibility that $P_{\mathrm{lead}}\left(t\right)\sim P_{\mathrm{lag}}{\left(t\right)}^3$ for $N=4$, although the correlations in this borderline case are tiny.

We introduce the threshold $q$-voter opinion dynamics where an agent, facing a binary choice, can change its mind when at least $q_0$ among $q$ neighbors share the opposite opinion. Otherwise, the agent can still change its mind with a certain probability $\varepsilon$. This threshold dynamics contemplates the possibility of persuasion by an influence group even when there is not full agreement among its members. In fact, individuals can follow their peers not only when there is unanimity ($q_0=q$) in the lobby group, as assumed in the $q$-voter model, but also, depending on the circumstances, when there is simple majority ($q_0>q/2$), Byzantine consensus ($q_0>2q/3$), or any minimal number $q_0$ among $q$. This realistic threshold gives place to emerging collective states and phase transitions which are not observed in the standard $q$ voter. The threshold $q_0$, together with the stochasticity introduced by $\varepsilon$, yields a phenomenology that mimics as particular cases the $q$ voter with stochastic drivings such as nonconformity and independence. In particular, nonconsensus majority states are possible, as well as mixed phases. Continuous and discontinuous phase transitions can occur, but also transitions from fluctuating phases into absorbing states.

This paper proposes a phase definition, named the vector field phase, which can be defined for systems with arbitrary finite dimension and is a monotonically increasing function of time. The proposed definition can properly quantify the dynamics in the flow direction, often associated with the null Lyapunov exponent. Numerical examples that use benchmark periodic and chaotic oscillators are discussed to illustrate some of the main features of the definition, which are that (i) phase information can be obtained either from the vector field or from a time series, (ii) it permits not only detection of phase synchronization but also quantification of it, and (iii) it can be used in the phase synchronization of very different oscillators.

The phase behavior and structural properties of a monolayer of hard particles is examined in such a confinement where the adsorbed particles are constrained to the surface of a narrow hard cylindrical pore. The diameter of the pore is chosen such that only first- and second-neighbor interactions occur between the hard particles. The transfer operator method of [Percus and Zhang, Mol. Phys. 69, 347 (1990)] is reformulated to obtain information about the structure of the monolayer. We have found that a true phase transition is not possible in the examined range of pore diameters. The monolayer of hard spheres undergoes a structural change from fluidlike order to a zigzaglike solid one with increasing surface density. The case of hard cylinders is different in the sense that a layering takes place continuously between a low-density one-row and a high-density two-row monolayer. Our results reveal a clear discrepancy with classical density functional theories, which do not distinguish smecticlike ordering in bulk from that in narrow periodic pores.

A series of laboratory experiments were performed to study the evolution of Peregrine breather (PB) in a wave flume in finite depth, and wave trains were initially generated in a region of quiescent water and then propagated into an adverse current region for which the current velocity strength gradually increased from zero to an approximately stable value. The PB is often considered as a prototype of oceanic freak waves that can focus wave energy into a single wave packet. In the experiment, the cases were selected with the relative water depths $k_{0}h$ ($k_0$ is the wave number in quiescent water and $h$ is the water depth) varying from 3.11 through 8.17, and the initial wave steepness $k_0{a}_0$ ($a_0$ is the background wave amplitude) ranges between 0.065 and 0.120. The experimental results show the persistence of the breather evolution dynamics even in the presence of strong opposing currents. We have shown that the characteristic spectrum of the PB persists even on strong currents, thus making it a viable characteristic for prediction of freak waves. It was also found that the adverse currents tend to shift the focusing point upstream compared to the cases without currents. Furthermore, it was found that uniform-depth adverse currents can reduce the breather extension in time domain.

Characteristic timescales and length scales for phase transitions of real materials are in ranges where a direct visualization is unfeasible. Therefore, model systems can be useful. Here, the crystallization process of a three-dimensional complex plasma under gravity conditions is considered where the system ranges up to a large extent into the bulk plasma. Time-resolved measurements exhibit the process down to a single-particle level. Primary clusters, consisting of particles in the solid state, grow vertically and, secondarily, horizontally. The box-counting method shows a fractal dimension of $d_f\approx 2.72$ for the clusters. This value gives a hint that the formation process is a combination of local epitaxial and diffusion-limited growth. The particle density and the interparticle distance to the nearest neighbor remain constant within the clusters during crystallization. All results are in good agreement with former observations of a single-particle layer.

We present experiments and numerical simulations which demonstrate that fully ionized, low-density plasma channels could be formed by hydrodynamic expansion of plasma columns produced by optical field ionization. Simulations of the hydrodynamic expansion of plasma columns formed in hydrogen by an axicon lens show the generation of 200 mm long plasma channels with axial densities of order $n_{\mathrm{e}}\left(0\right)=1\times {10}^{17}{\mathrm{cm}}^{-3}$ and lowest-order modes of spot size $W_{\mathrm{M}}\approx 40\mu \mathrm{m}$. These simulations show that the laser energy required to generate the channels is modest: of order 1 mJ per centimeter of channel. The simulations are confirmed by experiments with a spherical lens which show the formation of short plasma channels with $1.5\times {10}^{17}{\mathrm{cm}}^{-3}\lesssim n_{\mathrm{e}}\left(0\right)\lesssim 1\times {10}^{18}{\mathrm{cm}}^{-3}$ and $61\mu \mathrm{m}\gtrsim W_{\mathrm{M}}\gtrsim 33\mu \mathrm{m}$. Low-density plasma channels of this type would appear to be well suited as multi-GeV laser-plasma accelerator stages capable of long-term operation at high pulse repetition rates.

Lattice Boltzmann (LB) models used for the computation of fluid flows represented by the Navier-Stokes (NS) equations on standard lattices can lead to non-Galilean-invariant (GI) viscous stress involving cubic velocity errors. This arises from the dependence of their third-order diagonal moments on the first-order moments for standard lattices, and strategies have recently been introduced to restore Galilean invariance without such errors using a modified collision operator involving corrections to either the relaxation times or the moment equilibria. Convergence acceleration in the simulation of steady flows can be achieved by solving the preconditioned NS equations, which contain a preconditioning parameter that can be used to tune the effective sound speed, and thereby alleviating the numerical stiffness. In the present paper, we present a GI formulation of the preconditioned cascaded central-moment LB method used to solve the preconditioned NS equations, which is free of cubic velocity errors on a standard lattice, for steady flows. A Chapman-Enskog analysis reveals the structure of the spurious non-GI defect terms and it is demonstrated that the anisotropy of the resulting viscous stress is dependent on the preconditioning parameter, in addition to the fluid velocity. It is shown that partial correction to eliminate the cubic velocity defects is achieved by scaling the cubic velocity terms in the off-diagonal third-order moment equilibria with the square of the preconditioning parameter. Furthermore, we develop additional corrections based on the extended moment equilibria involving gradient terms with coefficients dependent locally on the fluid velocity and the preconditioning parameter. Such parameter dependent corrections eliminate the remaining truncation errors arising from the degeneracy of the diagonal third-order moments and fully restore Galilean invariance without cubic defects for the preconditioned LB scheme on a standard lattice. Several conclusions are drawn from the analysis of the structure of the non-GI errors and the associated corrections, with particular emphasis on their dependence on the preconditioning parameter. The GI preconditioned central-moment LB method is validated for a number of complex flow benchmark problems and its effectiveness to achieve convergence acceleration and improvement in accuracy is demonstrated.

A fluctuation relation is derived to extract the order parameter function $q\left(x\right)$ in weakly ergodic systems. The relation is based on measuring and classifying entropy production fluctuations according to the value of the overlap $q$ between configurations. For a fixed value of $q$, entropy production fluctuations are Gaussian distributed allowing us to derive the quasi-FDT so characteristic of aging systems. The theory is validated by extracting the $q\left(x\right)$ in various types of glassy models. It might be generally applicable to other nonequilibrium systems and experimental small systems.

We study the influence of thermal fluctuations on the buckling behavior of thin elastic capsules with spherical rest shape. Above a critical uniform pressure, an elastic capsule becomes mechanically unstable and spontaneously buckles into a shape with an axisymmetric dimple. Thermal fluctuations affect the buckling instability by two mechanisms. On the one hand, thermal fluctuations can renormalize the capsule's elastic properties and its pressure because of anharmonic couplings between normal displacement modes of different wavelengths. This effectively lowers its critical buckling pressure [Košmrlj and Nelson, Phys. Rev. X 7, 011002 (2017)]. On the other hand, buckled shapes are energetically favorable already at pressures below the classical buckling pressure. At these pressures, however, buckling requires to overcome an energy barrier, which only vanishes at the critical buckling pressure. In the presence of thermal fluctuations, the capsule can spontaneously overcome an energy barrier of the order of the thermal energy by thermal activation already at pressures below the critical buckling pressure. We revisit parameter renormalization by thermal fluctuations and formulate a buckling criterion based on scale-dependent renormalized parameters to obtain a temperature-dependent critical buckling pressure. Then we quantify the pressure-dependent energy barrier for buckling below the critical buckling pressure using numerical energy minimization and analytical arguments. This allows us to obtain the temperature-dependent critical pressure for buckling by thermal activation over this energy barrier. Remarkably, both parameter renormalization and thermal activation lead to the same parameter dependence of the critical buckling pressure on temperature, capsule radius and thickness, and Young's modulus. Finally, we study the combined effect of parameter renormalization and thermal activation by using renormalized parameters for the energy barrier in thermal activation to obtain our final result for the temperature-dependent critical pressure, which is significantly below the results if only parameter renormalization or only thermal activation is considered.

The spreading of epidemics is very much determined by the structure of the contact network, which may be impacted by the mobility dynamics of the individuals themselves. In confined scenarios where a small, closed population spends most of its time in localized environments and has easily identifiable mobility patterns—such as workplaces, university campuses, or schools—it is of critical importance to identify the factors controlling the rate of disease spread. Here, we present a discrete-time, metapopulation-based model to describe the transmission of susceptible-infected-susceptible-like diseases that take place in confined scenarios where the mobilities of the individuals are not random but, rather, follow clear recurrent travel patterns. This model allows analytical determination of the onset of epidemics, as well as the ability to discern which contact structures are most suited to prevent the infection to spread. It thereby determines whether common prevention mechanisms, as isolation, are worth implementing in such a scenario and their expected impact.

Living cells sense the mechanical properties of their surrounding environment and respond accordingly. Crawling cells detect the rigidity of their substratum and migrate in certain directions. They can be classified into two categories: slow-moving and fast-moving cell types. Slow-moving cell types, such as fibroblasts, smooth muscle cells, mesenchymal stem cells, etc., move toward rigid areas on the substratum in response to a rigidity gradient. However, there is not much information on rigidity sensing in fast-moving cell types whose size is ∼10 $\mu \mathrm{m}$ and migration velocity is ∼10 $\mu \mathrm{m}/\min$. In this study, we used both isotropic substrata with different rigidities and an anisotropic substratum that is rigid on the $x$ axis but soft on the $y$ axis to demonstrate rigidity sensing by fast-moving Dictyostelium cells and neutrophil-like differentiated HL-60 cells. Dictyostelium cells exerted larger traction forces on a more rigid isotropic substratum. Dictyostelium cells and HL-60 cells migrated in the “soft” direction on the anisotropic substratum, although myosin II–null Dictyostelium cells migrated in random directions, indicating that rigidity sensing of fast-moving cell types differs from that of slow types and is induced by a myosin II–related process.

A theoretical analysis of active motion on curved surfaces is presented in terms of a generalization of the telegrapher equation. Such a generalized equation is explicitly derived as the polar approximation of the hierarchy of equations obtained from the corresponding Fokker-Planck equation of active particles diffusing on curved surfaces. The general solution to the generalized telegrapher equation is given for a pulse with vanishing current as initial data. Expressions for the probability density and the mean squared geodesic displacement are given in the limit of weak curvature. As an explicit example of the formulated theory, the case of active motion on the sphere is presented, where oscillations observed in the mean squared geodesic displacement are explained.

We examine the connection between the singularities or quasisingularities in the solutions of the incompressible Navier-Stokes equation (INSE) and the local energy transfer and dissipation, in order to explore in detail how the former contributes to the phenomenon of intermittency. We do so by analyzing the velocity fields (a) measured in the experiments on the turbulent von Kármán swirling flow at high Reynolds numbers and (b) obtained from the direct numerical simulations of the INSE at a moderate resolution. To compute the local interscale energy transfer and viscous dissipation in experimental and supporting numerical data, we use the weak solution formulation generalization of the Kármán-Howarth-Monin equation. In the presence of a singularity in the velocity field, this formulation yields a nonzero dissipation (inertial dissipation) in the limit of an infinite resolution. Moreover, at finite resolutions, it provides an expression for local interscale energy transfers down to the scale where the energy is dissipated by viscosity. In the presence of a quasisingularity that is regularized by viscosity, the formulation provides the contribution to the viscous dissipation due to the presence of the quasisingularity. Therefore, our formulation provides a concrete support to the general multifractal description of the intermittency. We present the maps and statistics of the interscale energy transfer and show that the extreme events of this transfer govern the intermittency corrections and are compatible with a refined similarity hypothesis based on this transfer. We characterize the probability distribution functions of these extreme events via generalized Pareto distribution analysis and find that the widths of the tails are compatible with a similarity of the second kind. Finally, we make a connection between the topological and the statistical properties of the extreme events of the interscale energy transfer field and its multifractal properties.