The Biham-Middleton-Levine (BML) traffic model, a cellular automaton with eastbound and northbound cars moving by turns on a square lattice, has been an underpinning model in the study of collective behavior by cars, pedestrians, and even internet packages. Contrary to initial beliefs that the model exhibits a sharp phase transition from freely flowing to fully jammed, it has been reported that it shows intermediate stable phases, where jams and freely flowing traffic coexist, but there is no clear understanding of their origin. Here, we analyze the model as an anisotropic system with a preferred fluid direction (northeast) and find that it exhibits two differentiated phase transitions: the system is either longer in the flow direction (longitudinal) or perpendicular to it (transversal). The critical densities where these transitions occur enclose the density interval of intermediate states and can be approximated by mean-field analysis, all derived from the anisotropic exponent relating the longitudinal and transversal correlation lengths. Thus, we arrive at the interesting result that the puzzling intermediate states in the original model are just a superposition of these two different behaviors of the phase transition, solving by the way most mysteries behind the BML model, which turns out to be a paradigmatic example of such anisotropic critical systems.

The long-time effect of weak rotation on an internal solitary wave is the decay into inertia-gravity waves and the eventual emergence of a coherent, steadily propagating, nonlinear wave packet. There is currently no entirely satisfactory explanation as to why these wave packets form. Here the initial value problem is considered within the context of the Gardner-Ostrovsky, or rotation-modified extended Korteweg–de Vries, equation. The linear Gardner-Ostrovsky equation has maximum group velocity at a critical wave number, often called the zero-dispersion point. It is found here that a nonlinear splitting of the wave-number spectrum at the zero-dispersion point, where energy is shifted into the modulationally unstable regime of the Gardner-Ostrovsky equation, is responsible for the wave-packet formation. Numerical comparisons of the decay of a solitary wave in the Gardner-Ostrovsky equation and a derived nonlinear Schrödinger equation at the zero-dispersion point are used to confirm the spectral splitting.

We study the effect of pulsed driving and kicked driving of the interaction term on the nonequilibrium phase transition in the Dicke model. Within the framework of Floquet theory, we observe the emergence of new nontrivial phases on impingement by such periodic pulses. Notably, our study reveals that a greater control over the dynamical quantum criticality is possible through the variation of multiple parameters related to the pulse, as opposed to a single parameter control in a monochromatic drive. Furthermore, the probability of the system remaining trapped in a metastable state during the observed first order transition from the superradiant to normal phase is found to be higher for small number of kicks (or pulses) in comparison to the sinusoidal perturbation.

We analyze dynamic behavior of chemically reactive species in a porous medium, subject to anomalous transport. In this context, we present transport experiments in a refraction-index-matched, three-dimensional, water-saturated porous medium. A pH indicator (Congo red) was used as either a conservative or a reactive tracer, depending on the tracer solution pH relative to that of the background solution. The porous medium consisted of an acrylic polymer material formed as spherical beads that have pH-buffering capacity. The magnitude of reaction during transport through the porous medium was related to the color change of the Congo red, via image analysis. Here, we focused on point injection of the tracer into a macroscopically uniform flow field containing water at a pH different from that of the injected tracer. The setup yielded measurements of the temporally evolving spatial (local-in-space) concentration field. Parallel experiments with the same tracer, but without reactions (no changes in pH), enabled identification of the transport itself to be anomalous (non-Fickian); this was quantified by a continuous time random walk (CTRW) formulation. A CTRW particle tracking model was then used to quantify the spatial and temporal migration of both the conservative and reactive tracer plumes. Model parameters related to the anomalous transport were determined from the conservative tracer experiments. An additional term accounting for chemical reaction was established solely from analysis of the reactant concentrations, and significantly, no other fitting parameters were required. The measurements and analysis emphasized the localized nature of reaction, caused by small-scale concentration fluctuations and preferential pathways. In addition, a threshold radius for pH-controlled reactive transport processes was defined under buffering conditions, which delineated the region in which reactions occurred rapidly.

A stochastic process, when subject to resetting to its initial condition at a constant rate, generically reaches a nonequilibrium steady state. We study analytically how the steady state is approached in time and find an unusual relaxation mechanism in these systems. We show that as time progresses an inner core region around the resetting point reaches the steady state, while the region outside the core is still transient. The boundaries of the core region grow with time as power laws at late times with new exponents. Alternatively, at a fixed spatial point, the system undergoes a dynamical transition from the transient to the steady state at a characteristic space-dependent timescale $t^{*}\left(x\right)$. We calculate analytically in several examples the large deviation function associated with this spatiotemporal fluctuation and show that, generically, it has a second-order discontinuity at a pair of critical points characterizing the edges of the inner core. These singularities act as separatrices between typical and atypical trajectories. Our results are verified in the numerical simulations of several models, such as simple diffusion and fluctuating one-dimensional interfaces.

In this work we consider a simple random walk embedded in a generic branched structure and we find a close-form formula to calculate the hitting time $H\left(i,f\right)$ between two arbitrary nodes $i$ and $j$. We then use this formula to obtain the set of hitting times $\left\{H\left(i,f\right)\right\}$ for combs and their expectation values, namely, the mean first-passage time, where the average is performed over the initial node while the final node $f$ is given, and the global mean first-passage time, where the average is performed over both the initial and the final node. Finally, we discuss applications in the context of reaction-diffusion problems.

The flow of sand on a rough inclined plane is investigated experimentally. We directly show that a jammed layer of grains spontaneously forms below the avalanche. Its properties and its relation with the rheology of the flowing layer of grains are presented and discussed. In a second part, we study the dynamics of erosion and deposition solitary waves in the domain where they are transversally stable. We characterize their shapes and velocity profiles. We relate their translational velocity to the stopping height and to the mass trapped in the avalanche. Finally, we use the velocity profile to get insight into the rheology very close to the jamming limit.

The problem of estimating the size of the cooperatively rearranging regions (CRRs) in supercooled polymeric melts from an analysis of the $\alpha$-process in ordinary relaxation experiments is addressed. The mechanism whereby a CRR changes its configuration is viewed as consisting of two distinct steps: a reduced number of monomers reaches initially an activated state, allowing for some local rearrangement; then, the subsequent regression of the energy fluctuation may take place through the configurational degrees of freedom, thus allowing for further rearrangements on larger length scales. The latter are indeed those to which the well-known Donth's scheme refers. Local readjustments are described in the framework of a canonical formalism on a stationary ensemble of small-scale regions, distributed over all possible energy thresholds for rearrangement. Large-scale configurational changes, instead, are described as spontaneous processes. Two main regimes are envisaged, depending on whether the role played by the configurational degrees of freedom in the regression of the energy fluctuation is significant or not. It is argued that the latter case is related to the occurrence of an Arrhenian dependence of the central relaxation rate. Consistency with Donth's scheme is demonstrated, and data from the literature confirm the agreement of the two methods of analysis when configurational degrees of freedom are relevant for the fluctuation regression. Poly($n$-butyl methacrylate) is chosen in order to show how CRR size and temperature fluctuations at rearrangement can be estimated from stress relaxation experiments carried out by means of an atomic force microscopy setup. Cases in which the configurational pathway for regression is significantly hindered are considered. Relaxation in poly(dimethyl siloxane) confined in nanopores is taken as an example to suggest how a more complete view of the effects of configurational constraints would be possible if direct measurements of temperature fluctuations were combined with the proposed analysis.

Replica-exchange molecular dynamics (REMD) is a popular sampling method in the molecular simulation. By frequently exchanging the replicas at different temperatures, the molecule can jump out of the minima and sample efficiently in the conformational space. Although REMD has been shown to be practical in a lot of applications, it does have a critical limitation. All the replicas at all the temperatures must be simulated for a period between the replica-exchange steps. This may be problematic for the reaction with high free energy barriers. In that case, too many replicas are required in the simulation. To reduce the calculation quantity and improve its performance, in this paper we propose a modified REMD method. During the simulation, each replica at each temperature can stay in either the active or inactive state and only switch between the states at the exchange step. In the active state, the replica moves freely in the canonical ensemble by the normal molecular dynamics, and in the inactive state, the replica is frozen temporarily until the next exchange step. The number of the replicas in the active states (active replicas) depends on the number of CPUs in the computer. Using the additional inactive replicas, one can perform an REMD simulation in a wider temperature space. The practical applications show that the modified REMD method is reliable. With the same number of active replicas, this REMD method can produce a more reasonable free energy surface around the free energy minima than the standard REMD method.

Interdependent networks have been shown to be extremely vulnerable based on the percolation model. Parshani et al. [Europhys. Lett. 92, 68002 (2010)] further indicated that the more intersimilar networks are, the more robust they are to random failures. When traffic load is considered, how do the coupling patterns impact cascading failures in interdependent networks? This question has been largely unexplored until now. In this paper, we address this question by investigating the robustness of interdependent Erdös-Rényi random graphs and Barabási-Albert scale-free networks under either random failures or intentional attacks. It is found that interdependent Erdös-Rényi random graphs are robust yet fragile under either random failures or intentional attacks. Interdependent Barabási-Albert scale-free networks, however, are only robust yet fragile under random failures but fragile under intentional attacks. We further analyze the interdependent communication network and power grid and achieve similar results. These results advance our understanding of how interdependency shapes network robustness.

We develop accurate analytical relations for the droplet volume ratio, droplet length during breakup process, and pressure drop of asymmetric T junctions with a valve in each of the branches for producing unequal-sized droplets. An important advantage of this system is that after manufacturing the system, the size of the generated droplets can be changed simply by adjusting the valves. The results indicate that if the valve ratio is smaller than 0.65, the system enters a nonbreakup regime. Also the pressure drop does not depend on the time and decreases by increasing the valve ratio, namely, opening the degree of valve 1 to valve 2. In addition, the results reveal that by decreasing (increasing) the valve ratio, the droplet length of branch 1 decreases (increases) and the droplet length of branch 2 increases (decreases) linearly while the whole length of the droplet remains unchanged.

Basic flow generated in a viscous incompressible fluid by a “point” sink (source) of mass is revised. In practice, such flow can be modeled by sucking (pushing) fluid from a thin tube with a small porous sphere at one end. Intuitively, by sucking (pushing) fluid, one creates low (high) pressure near the origin and a positive (negative) radial pressure gradient drives the fluid to (from) the origin. A simple analysis, however, shows that the pressure distribution for both steady flows is the same. Then a question arises: How does the fluid “know” in what direction to flow? To explain this “paradox” an unsteady flow is considered and the pressure terms responsible for the flow direction are derived.

An analog to the equations of compressible flow that is based on the inviscid Burgers equation is utilized to investigate the effect of spatial discreteness of energy release on the propagation of a detonation wave. While the traditional Chapman-Jouguet (CJ) treatment of a detonation wave assumes that the energy release of the medium is homogeneous through space, the system examined here consists of sources represented by $\delta$ functions embedded in an otherwise inert medium. The sources are triggered by the passage of the leading shock wave following a delay that is either of fixed period or randomly generated. The solution for wave propagation through a large array $({10}^3–{10}^4)$ of sources in one dimension can be constructed without the use of a finite difference approximation by tracking the interaction of sawtooth-profiled waves for which an analytic solution is available. A detonation-like wave results from the interaction of the shock and rarefaction waves generated by the sources. The measurement of the average velocity of the leading shock front for systems of both regular, fixed-period and randomized sources is found to be in close agreement with the velocity of the equivalent CJ detonation in a uniform medium, wherein the sources have been spatially homogenized. This result may have implications for the applicability of the CJ criterion to detonations in highly heterogeneous media (e.g., polycrystalline, solid explosives) and unstable detonations with a transient and multidimensional structure (e.g., gaseous detonation waves).

In this paper we consider resonant ion acceleration by a plasma jet originating from the magnetic reconnection region. Such jets propagate in the background magnetic field with significantly curved magnetic-field lines. Decoupling of ion and electron motions at the leading edge of the jet results in generation of strong electrostatic fields. Ions can be trapped by this field and get accelerated along the jet front. This mechanism of resonant acceleration resembles surfing acceleration of charged particles at a shock wave. To describe resonant acceleration of ions, we use adiabatic theory of resonant phenomena. We show that particle motion along the curved field lines significantly influences the acceleration rate. The maximum gain of energy is determined by the particle's escape from the system due to this motion. Applications of the proposed mechanism to charged-particle acceleration in the planetary magnetospheres and the solar corona are discussed.

Quantum discord is a measure of quantum correlations beyond the entanglement-separability paradigm. It is conceptualized by using the von Neumann entropy as a measure of disorder. We introduce a class of quantum correlation measures as differences between total and classical correlations, in a shared quantum state, in terms of the sandwiched relative Rényi and Tsallis entropies. We compare our results with those obtained by using the traditional relative entropies. We find that the measures satisfy all the plausible axioms for quantum correlations. We evaluate the measures for shared pure as well as paradigmatic classes of mixed states. We show that the measures can faithfully detect the quantum critical point in the transverse quantum Ising model and find that they can be used to remove an unquieting feature of nearest-neighbor quantum discord in this respect. Furthermore, the measures provide better finite-size scaling exponents of the quantum critical point than the ones for other known order parameters, including entanglement and information-theoretic measures of quantum correlations.

When many independent users try to route traffic through a network, the flow can easily become suboptimal as a consequence of congestion of the most efficient paths. The degree of this suboptimality is quantified by the so-called price of anarchy (POA), but so far there are no general rules for when to expect a large POA in a random network. Here I address this question by introducing a simple model of flow through a network with randomly placed congestible and incongestible links. I show that the POA is maximized precisely when the fraction of congestible links matches the percolation threshold of the lattice. Both the POA and the total cost demonstrate critical scaling near the percolation threshold.

In this work we study the phase behavior of magnetic particles suspended in a simple nonmagnetic solvent. Magnetic particles are modelled as spherical particles carrying a three-dimensional, classical Heisenberg spin, whereas solvent molecules are treated as spherically symmetric Lennard-Jones particles. The binary mixture of magnetic particles and solvent is studied within the framework of classical density functional theory (DFT). Within DFT pair correlations are treated at the modified mean-field level at which they are approximated by orientation dependent Mayer $f$ functions. In the absence of an external magnetic field four generic types of phase diagrams are observed depending on the concentration of magnetic particles. In this case we observe liquid-liquid phase coexistence between an orientationally ordered (polarized) and a disordered phase characterized by slightly different concentrations of magnetic particles. Liquid-liquid phase coexistence is suppressed by an external field and vanishes completely if the strength of the field is sufficiently large.

We investigate nonequilibrium underdamped Langevin dynamics of Brownian particles that interact through a harmonic potential with coupling constant $K$ and are in thermal contact with two heat baths at different temperatures. The system is characterized by a net heat flow and an entropy production in the steady state. We compare the entropy production of the harmonic system with that of Brownian particles linked with a rigid rod. The harmonic system may be expected to reduce to the rigid rod system in the infinite $K$ limit. However, we find that the harmonic system in the $K\to \infty$ limit produces more entropy than the rigid rod system. The harmonic system has the center-of-mass coordinate as a slow variable and the relative coordinate as a fast variable. By identifying the contributions of the degrees of freedom to the total entropy production, we show that the hidden entropy production by the fast variable is responsible for the extra entropy production. We discuss the $K$ dependence of each contribution.

Mixed feedback loops combining transcriptional and posttranscriptional regulations are common in cellular regulatory networks. They consist of two genes, encoding a transcription factor and a small noncoding RNA (sRNA), which mutually regulate each other's expression. We present a theoretical and numerical study of coherent mixed feedback loops of this type, in which both regulations are negative. Under suitable conditions, these feedback loops are expected to exhibit bistability, namely, two stable states, one dominated by the transcriptional repressor and the other dominated by the sRNA. We use deterministic methods based on rate equation models, in order to identify the range of parameters in which bistability takes place. However, the deterministic models do not account for the finite lifetimes of the bistable states and the spontaneous, fluctuation-driven transitions between them. Therefore, we use stochastic methods to calculate the average lifetimes of the two states. It is found that these lifetimes strongly depend on rate coefficients such as the transcription rates of the transcriptional repressor and the sRNA. In particular, we show that the fraction of time the system spends in the sRNA-dominated state follows a monotonically decreasing sigmoid function of the transcriptional repressor transcription rate. The biological relevance of these results is discussed in the context of such mixed feedback loops in Escherichia coli. It is shown that the fluctuation-driven transitions and the dependence of some rate coefficients on the biological conditions enable the cells to switch to the state which is better suited for the existing conditions and to remain in that state as long as these conditions persist.