Chimera states in the systems of nonlocally coupled phase oscillators are considered stable in the continuous limit of spatially distributed oscillators. However, it is reported that in the numerical simulations without taking such limit, chimera states are chaotic transient and finally collapse into the completely synchronous solution. In this Rapid Communication, we numerically study chimera states by using the coupling function different from the previous studies and obtain the result that chimera states can be stable even without taking the continuous limit, which we call the persistent chimera state.

We report a charge-patterning phase transition on two-dimensional square lattices of titratable sites, here regarded as protonation sites, placed in a low-dielectric medium just below the planar interface between this medium and a salt solution. We calculate the work-of-charging matrix of the lattice with use of a linear Debye-Hückel model, as input to a grand-canonical partition function for the distribution of occupancy patterns. For a large range of parameter values, this model exhibits an approximate inverse cubic power-law decrease of the voltage produced by an individual charge, as a function of its in-lattice separation from neighboring titratable sites. Thus, the charge coupling voltage biases the local probabilities of proton binding as a function of the occupancy of sites for many neighbors beyond the nearest ones. We find that even in the presence of these longer-range interactions, the site couplings give rise to a phase transition in which the site occupancies exhibit an alternating, checkerboard pattern that is an analog of antiferromagnetic ordering. The overall strength $W$ of this canonical charge coupling voltage, per unit charge, is a function of the Debye length, the charge depth, the Bjerrum length, and the dielectric coefficients of the medium and the solvent. The alternating occupancy transition occurs above a curve of thermodynamic critical points in the $(p\text{H}\text{−}pK,W)$ plane, the curve representing a charge-regulation analog of variation of the Néel temperature of an Ising antiferromagnet as a function of an applied, uniform magnetic field. The analog of a uniform magnetic field in the antiferromagnet problem is a combination of $p\text{H}\text{−}pK$ and $W$, and $1/W$ is the analog of the temperature in the antiferromagnet problem. We use Monte Carlo simulations to study the occupancy patterns of the titratable sites, including interactions out to the 37th nearest-neighbor category (a distance of $\sqrt{74}$ lattice constants), first validating simulations through comparison with exact and approximate results for the nearest-neighbor case. We then use the simulations to map the charge-patterning phase boundary in the $(p\text{H}\text{−}pK,W)$ plane. The physical parameters that determine $W$ provide a framework for identifying and designing real surfaces that could exhibit charge-patterning phase transitions.

The likelihood that a mutant fixates in the wild population, i.e., fixation probability, has been intensively studied in evolutionary game theory, where individuals' fitness is frequency dependent. However, it is of limited interest when it takes long to take over. Thus the speed of evolution becomes an important issue. In general, it is still unclear how fixation times are affected by the population structure, although the fixation times have already been addressed in the well-mixed populations. Here we theoretically address this issue by pair approximation and diffusion approximation on regular graphs. It is shown (i) that under neutral selection, both unconditional and conditional fixation time are shortened by increasing the number of neighbors; (ii) that under weak selection, for the simplified prisoner's dilemma game, if benefit-to-cost ratio exceeds the degree of the graph, then the unconditional fixation time of a single cooperator is slower than that in the neutral case; and (iii) that under weak selection, for the conditional fixation time, limited neighbor size dilutes the counterintuitive stochastic slowdown which was found in well-mixed populations. Interestingly, we find that all of our results can be interpreted as that in the well-mixed population with a transformed payoff matrix. This interpretation is also valid for both death-birth and birth-death processes on graphs. This interpretation bridges the fixation time in the structured population and that in the well-mixed population. Thus it opens the avenue to investigate the challenging fixation time in structured populations by the known results in well-mixed populations.

Analysis of the temperature dependence of the structural relaxation time ${\tau }_{\alpha }(T)$ in supercooled liquids revealed a qualitatively distinct feature—a sharp, cusplike maximum in the second derivative of $\log {\tau }_{\alpha }\left(T\right)$at some $T_{\max }$. It suggests that the super-Arrhenius temperature dependence of ${\tau }_{\alpha }\left(T\right)$ in glass-forming liquids eventually crosses over to an Arrhenius behavior at $T<T_{\max }$, and there is no divergence of ${\tau }_{\alpha }\left(T\right)$ at nonzero $T$. $T_{\max }$ can be above or below $T_{\mathrm{g}}$, depending on the sensitivity of $\tau (T)$ to a change in the liquid's density quantified by the exponent $\gamma$ in the scaling ${\tau }_{\alpha }(T)\sim \exp (A/T{\rho }^{-\gamma })$. These results might turn the discussion of the glass transition in a different direction—toward the origin of the limiting activation energy for structural relaxation at low $T$.

We experimentally study the structural properties of a sediment of size distributed colloids. By determining each particle size using a size estimation algorithm, we are able to investigate the relationship between local environment and local ordering. Our results show that ordered environments of particles tend to generate where the local particle-size distribution is within 5%. In addition, we show that particles whose size is close to the average size have 12 coordinate neighbors, which matches the coordination number of the fcc and hcp crystals. On the other hand, bcc structures are observed around larger particles. Our results represent experiments to show a size dependence of the specific ordering in colloidal systems.

Hydrophobic poly(methyl methacrylate) (PMMA) colloidal particles, when dispersed in oil with a relatively high dielectric constant, can become highly charged. In the presence of an interface with a conducting aqueous phase, image-charge effects lead to strong binding of colloidal particles to the interface, even though the particles are wetted very little by the aqueous phase. We study both the behavior of individual colloidal particles as they approach the interface and the interactions between particles that are already interfacially bound. We demonstrate that using particles which are minimally wetted by the aqueous phase allows us to isolate and study those interactions which are due solely to charging of the particle surface in oil. Finally, we show that these interactions can be understood by a simple image-charge model in which the particle charge $q$ is the sole fitting parameter.

We study the effects of uniform time delays on the extreme fluctuations in stochastic synchronization and coordination problems with linear couplings in complex networks. We obtain the average size of the fluctuations at the nodes from the behavior of the underlying modes of the network. We then obtain the scaling behavior of the extreme fluctuations with system size, as well as the distribution of the extremes on complex networks, and compare them to those on regular one-dimensional lattices. For large complex networks, when the delay is not too close to the critical one, fluctuations at the nodes effectively decouple, and the limit distributions converge to the Fisher-Tippett-Gumbel density. In contrast, fluctuations in low-dimensional spatial graphs are strongly correlated, and the limit distribution of the extremes is the Airy density. Finally, we also explore the effects of nonlinear couplings on the stability and on the extremes of the synchronization landscapes.

We consider the coaction of two distinct noise sources on the activation process of a single excitable unit and two interacting excitable units, which are mathematically described by the Fitzhugh-Nagumo equations. We determine the most probable activation paths around which the corresponding stochastic trajectories are clustered. The key point lies in introducing appropriate boundary conditions that are relevant for a class II excitable unit, which can be immediately generalized also to scenarios involving two coupled units. We analyze the effects of the two noise sources on the statistical features of the activation process, in particular demonstrating how these are modified due to the linear or nonlinear form of interactions. Universal properties of the activation process are qualitatively discussed in the light of a stochastic bifurcation that underlies the transition from a stochastically stable fixed point to continuous oscillations.

We study the activation process in large assemblies of type II excitable units whose dynamics is influenced by two independent noise terms. The mean-field approach is applied to explicitly demonstrate that the assembly of excitable units can itself exhibit macroscopic excitable behavior. In order to facilitate the comparison between the excitable dynamics of a single unit and an assembly, we introduce three distinct formulations of the assembly activation event. Each formulation treats different aspects of the relevant phenomena, including the thresholdlike behavior and the role of coherence of individual spikes. Statistical properties of the assembly activation process, such as the mean time-to-first pulse and the associated coefficient of variation, are found to be qualitatively analogous for all three formulations, as well as to resemble the results for a single unit. These analogies are shown to derive from the fact that global variables undergo a stochastic bifurcation from the stochastically stable fixed point to continuous oscillations. Local activation processes are analyzed in the light of the competition between the noise-led and the relaxation-driven dynamics. We also briefly report on a system-size antiresonant effect displayed by the mean time-to-first pulse.

We find a dark component guiding the practically interesting bright-bright vector one-soliton to two different parametric domains giving rise to different physical situations by constructing a more general form of three-component dark-bright-bright mixed vector one-soliton solution of the generalized Manakov model with nine free real parameters. Moreover our main investigation of the collision dynamics of such mixed vector solitons by constructing the multisoliton solution of the generalized Manakov model with the help of Hirota technique reveals that the dark-bright-bright vector two-soliton supports energy-exchange collision dynamics. In particular the dark component preserves its initial form and the energy-exchange collision property of the bright-bright vector two-soliton solution of the Manakov model during collision. In addition the interactions between bound state dark-bright-bright vector solitons reveal oscillations in their amplitudes. A similar kind of breathing effect was also experimentally observed in the Bose-Einstein condensates. Some possible ways are theoretically suggested not only to control this breathing effect but also to manage the beating, bouncing, jumping, and attraction effects in the collision dynamics of dark-bright-bright vector solitons. The role of multiple free parameters in our solution is examined to define polarization vector, envelope speed, envelope width, envelope amplitude, grayness, and complex modulation of our solution. It is interesting to note that the polarization vector of our mixed vector one-soliton evolves in sphere or hyperboloid depending upon the initial parametric choices.

We present hydrodynamic and magnetohydrodynamic (MHD) simulations of liquid sodium flow with the PLUTO compressible MHD code to investigate influence of magnetic boundary conditions on the collimation of helicoidal motions. We use a simplified cartesian geometry to represent the flow dynamics in the vicinity of one cavity of a multiblades impeller inspired by those used in the Von-Kármán-sodium (VKS) experiment. We show that the impinging of the large-scale flow upon the impeller generates a coherent helicoidal vortex inside the blades, located at a distance from the upstream blade piloted by the incident angle of the flow. This vortex collimates any existing magnetic field lines leading to an enhancement of the radial magnetic field that is stronger for ferromagnetic than for conducting blades. The induced magnetic field modifies locally the velocity fluctuations, resulting in an enhanced helicity. This process possibly explains why dynamo action is more easily triggered in the VKS experiment when using soft iron impellers.

Plastic materials (CH) doped with mid-$Z$ elements are used as ablators in inertial confinement fusion (ICF) capsules and in their surrogates. Hugoniot equation of state (EOS) and electronic properties of CH doped with germanium (at 2.5% and 13% dopant fractions) are investigated experimentally up to 7 Mbar using velocity and reflectivity measurements of shock fronts on the GEKKO laser at Osaka University. Reflectivity and temperature measurements were updated using a quartz standard. Shocked quartz reflectivity was measured at 532 and 1064 nm. Theoretical investigation of shock pressure and reflectivity was then carried out by ab initio simulations using the quantum molecular dynamics (QMD) code abinit and compared with tabulated average atom EOS models. We find that shock states calculated by QMD are in better agreement with experimental data than EOS models because of a more accurate description of ionic structure. We finally discuss electronic properties by comparing reflectivity data to a semiconductor gap closure model and to QMD simulations.

We apply the event-chain Monte Carlo algorithm to the three-dimensional ferromagnetic Heisenberg model. The algorithm is rejection-free and also realizes an irreversible Markov chain that satisfies global balance. The autocorrelation functions of the magnetic susceptibility and the energy indicate a dynamical critical exponent $z\approx 1$ at the critical temperature, while that of the magnetization does not measure the performance of the algorithm. We show that the event-chain Monte Carlo algorithm substantially reduces the dynamical critical exponent from the conventional value of $z\simeq 2$.

Labyrinthine patterns arise in two-dimensional physical systems submitted to competing interactions, in fields ranging from solid-state physics to hydrodynamics. For systems of interacting particles, labyrinthine and stripe phases were studied in the context of colloidal particles confined into a monolayer, both numerically by means of Monte Carlo simulations and experimentally using superparamagnetic particles. Here we report an experimental observation of a labyrinthine phase in an out-of-equilibrium system constituted of macroscopic particles. Once sufficiently magnetized, they organize into short chains of particles in contact and randomly orientated. We characterize the transition from a granular gas state towards a solid labyrinthine phase, as a function of the ratio of the interaction strength to the kinetic agitation. The spatial local structure is analyzed by means of accurate particle tracking. Moreover, we explain the formation of these chains using a simple model.

We study impulse response in one-dimensional homogeneous microgranular chains on a linear elastic substrate. Microgranular interactions are analytically described by the Schwarz contact model which includes nonlinear compressive as well as snap-to and from-contact adhesive effects forming a hysteretic loop in the force deformation relationship. We observe complex transient dynamics, including disintegration of solitary pulses, local clustering, and low-to-high-frequency energy transfers resulting in enhanced energy dissipation. We study in detail the underlying dynamics of cluster formation in the impulsively loaded medium and relate enhanced energy dissipation to the rate of cluster formation. These unusual and interesting dynamical phenomena are shown to be robust over a range of physically feasible conditions and are solely scale effects since they are attributed to surface forces, which have no effect at the macroscale. We establish a universal relation between the reclustering rate and the effective damping in these systems. Our findings demonstrate that scale effects generating new nonlinear features can drastically affect the dynamics and acoustics of microgranular materials.

Epidemics in networks can be affected by cooperation in transmission of infection and also connectivity between nodes. An interplay between these two properties and their influence on epidemic spread are addressed in the paper. A particular type of cooperative effects (called synergy effects) is considered, where the transmission rate between a pair of nodes depends on the number of infected neighbors. The connectivity effects are studied by constructing networks of different topology, starting with lattices with only local connectivity and then with networks that have both local and global connectivity obtained by random bond-rewiring to nodes within a certain distance. The susceptible-infected-removed epidemics were found to exhibit several interesting effects: (i) for epidemics with strong constructive synergy spreading in networks with high local connectivity, the bond rewiring has a negative role in epidemic spread, i.e., it reduces invasion probability; (ii) in contrast, for epidemics with destructive or weak constructive synergy spreading on networks of arbitrary local connectivity, rewiring helps epidemics to spread; (iii) and, finally, rewiring always enhances the spread of epidemics, independent of synergy, if the local connectivity is low.

Social activities display bursty behavior characterized by heavy-tailed interevent time distributions. We examine the bursty behavior of airplanes' arrivals in hub airports. The analysis indicates that the air transportation system universally follows a power-law interarrival time distribution with an exponent $\alpha =2.5$ and an exponential cutoff. Moreover, we investigate the mechanism of this bursty behavior by introducing a simple model to describe it. In addition, we compare the extent of the hub-and-spoke structure and the burstiness of various airline networks in the system. Remarkably, the results suggest that the hub-and-spoke network of the system and the carriers' strategy to facilitate transit are the origins of this universality.

Systems in which the heat flux depends on the direction of the flow are said to present thermal rectification. This effect has attracted much theoretical and experimental interest in recent years. However, in most theoretical models the effect is found to vanish in the thermodynamic limit, in disagreement with experiment. The purpose of this paper is to show that the rectification may be restored by including an energy-conserving noise which randomly flips the velocity of the particles with a certain rate $\lambda$. It is shown that as long as $\lambda$ is nonzero, the rectification remains finite in the thermodynamic limit. This is illustrated in a classical harmonic chain subject to a quartic pinning potential (the ${\mathrm{\Phi }}^4$ model) and coupled to heat baths by Langevin equations.

Recently we showed that the critical nonequilibrium relaxation in the Swendsen-Wang algorithm is widely described by the stretched-exponential relaxation of physical quantities in the Ising or Heisenberg models. Here we make a similar analysis in the Berezinsky-Kosterlitz-Thouless phase transition in the two-dimensional (2D) $XY$ model and in the first-order phase transition in the 2D $q=5$ Potts model and find that these phase transitions are described by the simple exponential relaxation and power-law relaxation of physical quantities, respectively. We compare the relaxation behaviors of these phase transitions with those of the second-order phase transition in the three- and four-dimensional $XY$ models and in the 2D $q$-state Potts models for $2\le q\le 4$ and show that the species of phase transitions can be clearly characterized by the present analysis. We also compare the size dependence of relaxation behaviors of the first-order phase transition in the 2D $q=5$ and 6 Potts models and propose a quantitative criterion on “weakness” of the first-order phase transition.

In this work we study opinion formation in a population participating in a public debate with two distinct choices. We consider three distinct mechanisms of social interactions and individuals' behavior: conformity, nonconformity, and inflexibility. The conformity is ruled by the majority-rule dynamics, whereas the nonconformity is introduced in the population as an independent behavior, implying the failure of attempted group influence. Finally, the inflexible agents are introduced in the population with a given density. These individuals present a singular behavior, in a way that their stubbornness makes them reluctant to change their opinions. We consider these effects separately and all together, with the aim to analyze the critical behavior of the system. We perform numerical simulations in some lattice structures and for distinct population sizes. Our results suggest that the different formulations of the model undergo order-disorder phase transitions in the same universality class as the Ising model. Some of our results are complemented by analytical calculations.