We investigate a phase model that includes both locally attractive and globally repulsive coupling in one dimension. This model exhibits nontrivial spatiotemporal patterns that have not been observed in systems that contain only local or global coupling. Depending on the relative strengths of the local and global coupling and on the form of global coupling, the system can show a spatially uniform state (in-phase synchronization), a monotonically increasing state (traveling wave), and three types of oscillations of relative phase difference. One of the oscillations of relative phase difference has the characteristic of being locally unstable but globally attractive. That is, any small perturbation to the periodic orbit in phase space destroys its periodic motion, but after a long time the system returns to the original periodic orbit. This behavior is closely related to the emergence of saddle two-cluster states for global coupling only, which are connected to each other by attractive heteroclinic orbits. The mechanism of occurrence of this type of oscillation is discussed.

The nonlinear optical Galton board (NLOGB), a quantum walk like (but nonlinear) discrete time quantum automaton, is shown to admit a complex evolution leading to long time thermalized states. The continuous limit of the Galton board is derived and shown to be a nonlinear Dirac equation (NLDE). The (Galerkin-truncated) NLDE evolution is shown to thermalize toward states qualitatively similar to those of the NLOGB. The NLDE conserved quantities are derived and used to construct a stochastic differential equation converging to grand canonical distributions that are shown to reproduce the (microcanonical) NLDE thermalized statistics. Both the NLOGB and the Galerkin-truncated NLDE are thus demonstrated to exhibit spontaneous thermalization.

Large amplitude traveling waves of the Korteweg–de-Vries (KdV) equation can be excited and controlled by a chirped frequency driving perturbation. The process involves capturing the wave into autoresonance (a continuous nonlinear synchronization) with the drive by passage through the linear resonance in the problem. The transition to autoresonance has a sharp threshold on the driving amplitude. In all previously studied autoresonant problems the threshold was found via a weakly nonlinear theory and scaled as ${\alpha }^{3/4},\alpha$ being the driving frequency chirp rate. It is shown that this scaling is violated in a long wavelength KdV limit because of the increased role of the nonlinearity in the problem. A fully nonlinear theory describing the phenomenon and applicable to all wavelengths is developed.

Following the exact decomposition in eigenstates of helicity for the Navier-Stokes equations in Fourier space [F. Waleffe, Phys. Fluids A 4, 350 (1992)], we introduce a modified version of helical shell models for turbulence with nonlocal triadic interactions. By using both an analytical argument and numerical simulation, we show that there exists a class of models, with a specific helical structure, that exhibits a statistically stable inverse energy cascade, in close analogy with that predicted for the Navier-Stokes equations restricted to the same helical interactions. We further support the idea that turbulent energy transfer is the result of a strong entanglement among triads possessing different transfer properties.

We investigate soliton mobility in the disordered Ablowitz-Ladik (AL) model and the standard nonlinear Schrödinger (NLS) lattice with the help of an effective potential generalizing the Peierls-Nabarro potential. This potential results from a deviation from integrability, which is due to randomness for the AL model, and both randomness and lattice discreteness for the NLS lattice. The statistical properties of such a potential are analyzed, and it is shown how the soliton mobility is affected by its size. The usefulness of this effective potential in studying soliton dynamics is demonstrated numerically. Furthermore, we propose two ways to enhance soliton transport in the presence of disorder: one is to use specific realizations of randomness, and the other is to consider a specific soliton pair.

We study the fluctuation-induced Casimir interactions in colloidal suspensions, especially between colloids immersed in a binary liquid close to its critical demixing point. To simulate these systems, we present a highly efficient cluster Monte Carlo algorithm based on geometric symmetries of the Hamiltonian. Utilizing the principle of universality, the medium is represented by an Ising system while the colloids are areas of spins with fixed orientation. Our results for the Casimir interaction potential between two particles at the critical point in two dimensions perfectly agree with the exact predictions. However, we find that in finite systems the behavior strongly depends on whether the $Z_2$ symmetry of the system is broken by the particles. We present Monte Carlo results for the three-body Casimir interaction potential and take a close look onto the case of one particle in the vicinity of two adjacent particles, which can be calculated from the two-particle interaction by a conformal mapping. These results emphasize the failure of the common decomposition approach for many-particle critical Casimir interactions.

Leaving the joint enterprise when defection is unveiled is always a viable option to avoid being exploited. Although loner strategy helps the population not to be trapped into the tragedy of the commons state, it could offer only a modest income for nonparticipants. In this paper we demonstrate that showing some tolerance toward defectors could not only save cooperation in harsh environments but in fact results in a surprisingly high average payoff for group members in public goods games. Phase diagrams and the underlying spatial patterns reveal the high complexity of evolving states where cyclic dominant strategies or two-strategy alliances can characterize the final state of evolution. We identify microscopic mechanisms which are responsible for the superiority of global solutions containing tolerant players. This phenomenon is robust and can be observed both in well-mixed and in structured populations highlighting the importance of tolerance in our everyday life.

The nonequilibrium dynamics of an isolated quantum system after a sudden quench to a dynamical critical point is expected to be characterized by scaling and universal exponents due to the absence of time scales. We explore these features for a quench of the parameters of a Hamiltonian with $O\left(N\right)$ symmetry, starting from a ground state in the disordered phase. In the limit of infinite $N$, the exponents and scaling forms of the relevant two-time correlation functions can be calculated exactly. Our analytical predictions are confirmed by the numerical solution of the corresponding equations. Moreover, we find that the same scaling functions, yet with different exponents, also describe the coarsening dynamics for quenches below the dynamical critical point.

We present a two-dimensional lattice model of water to study the effects of ion hydration on the properties of water. We map the water molecules as lattice particles consisting of a single oxygen atom at the center of a site and two hydrogen atoms on each side. The internal state of the system, such as the dipole moment at a site, is defined with respect to the location of the hydrogen atoms at the site depending on their role in hydrogen bonds (H bonds) being a donor or an acceptor. We study the influence of the charge and the radius of the ion on the insertion energy and on the H bonds in the first and second hydration layers around the ion and in the bulk. In particular we analyze how the competing interactions of the short-ranged H bonds and the long-ranged electrostatics influence the hydration properties. The role of the ion both as a source of the electrostatic interactions as well as a defect is also discussed. Our model also shows the well-known fact that the polarizability of the water molecules destroys the hydrogen bond network and increases the dipole moment of the molecules near the ion.

Using a thermodynamic integration scheme, we compute the free energy cost per unit area, $\gamma$, of forming an interface between a crystal and a frozen structured wall, formed by particles frozen into the same equilibrium structure as the crystal. Even though the structure and potential energy of the crystalline phase in the vicinity of the wall is the same as in the bulk, $\gamma$ has a nonzero value and increases with increasing density of the crystal and the wall. Investigating the effect of several interaction potentials between the particles, we observe a positive $\gamma$ at all crystalline densities if the potential is purely repulsive. For models with attractive interactions, such as the Lennard-Jones potential, a negative value for $\gamma$ is obtained at low densities. A qualitative explanation for the change of sign of $\gamma$ when going from repulsive to attractive interactions, at low crystal densities, is suggested.

We propose a theoretical description of the interaction mediated by a nematic-liquid-crystal host with different Frank elastic constants. A general expression for the energy of such an interaction between colloidal particles of arbitrary size and shape suspended in a homeotropic cell is obtained. In the cells of large thickness, the presented potential converges to that found previously for small particles in the nematic bulk. In general, our results confirm the validity of the one-constant approximation for weakly elastically anisotropic nematic liquid crystals. For nematics with a high splay-to-bend ratio we predict a larger range of the interaction. Using the dependence of this range on the elastic constants, we show that there exists a qualitative similarity between the interactions in a nematic and in a smectic-$A$ phase. It manifests itself, in particular, in a decrease of the angle between a chain of quadrupole particles and the uniform far-field director across a nematic–smectic-$A$ phase transition. We also demonstrate that the anisotropy of the elastic constants can lead to the formation of thermodynamically stable linear superstructures of asymmetric particles (elastic monopoles) with large, compared to usual dipole chains, interparticle distances.

We study a binary dynamical process that is a representation of the voter model with two candidates and opinion makers. The voters are represented by nodes of a network of social contacts with internal states labeled 0 or 1 and nodes that are connected can influence each other. The network is also perturbed by opinion makers, a set of external nodes whose states are frozen in 0 or 1 and that can influence all nodes of the network. The quantity of interest is the probability of finding $m$ nodes in state 1 at time $t$. Here we study this process on star networks, which are simple representations of hubs found in complex systems, and compare the results with those obtained for networks that are fully connected. In both cases a transition from disordered to ordered equilibrium states is observed as the number of external nodes becomes small. For fully connected networks the probability distribution becomes uniform at the critical point. For star networks, on the other hand, we show that the equilibrium distribution splits in two peaks, reflecting the two possible states of the central node. We obtain approximate analytical solutions for the equilibrium distribution that clarify the role of the central node in the process. We show that the network topology also affects the time scale of oscillations in single realizations of the dynamics, which are much faster for the star network. Finally, extending the analysis to two stars we compare our results with simulations in simple scale-free networks.

Chimera states consisting of domains of coherently and incoherently oscillating oscillators in a two-dimensional periodic array of nonlocally coupled phase oscillators are studied. In addition to the one-dimensional chimera states familiar from one spatial dimension, two-dimensional structures termed twisted chimera states and spiral wave chimera states are identified in simulations. The properties of many of these states, including stability, are determined using an evolution equation for a complex order parameter and are found to be in agreement with the simulations.

A central topic in Hele-Shaw flow research is the inclusion of physical effects on the interface between fluids. In this context, the addition of surface tension restrains the emergence of high interfacial curvatures, while consideration of kinetic undercooling effects inhibits the occurrence of high interfacial velocities. By connecting kinetic undercooling to the action of the dynamic contact angle, we show in a quantitative manner that the kinetic undercooling contribution varies as a linear function of the normal velocity at the interface. A perturbative weakly nonlinear analysis is employed to extract valuable information about the influence of kinetic undercooling on the shape of the emerging fingered structures. Under radial Hele-Shaw flow, it is found that kinetic undercooling delays, but does not suppress, the development of finger tip-broadening and finger tip-splitting phenomena. In addition, our results indicate that kinetic undercooling plays a key role in determining the appearance of tip splitting in rectangular Hele-Shaw geometry.

The complete cascade of bifurcations from steady to chaotic convection, as the Rayleigh number is varied, is considered numerically inside an air-filled differentially heated cavity. The system is assumed to be two-dimensional and is invariant under a generalized reflection about the center of the cavity. In the neighborhood of several codimension-two points, two main routes emerge, characterized by different symmetries of the first oscillatory eigenstate. Along these two competing routes, different sequences of bifurcations and symmetry breakings lead from the steady base flow to the hyperchaotic regime. Several families of two- and three-frequency tori have been identified via the computation of the leading Lyapunov exponents. Modal structures extracted from time series reveal the occurrence of slow internal oscillations in the center of the cavity and faster wall modes confined to vertical boundary layers. Further quasiperiodicity windows have been detected on each route. The different regimes eventually disappear in a boundary crisis in favor of a single, globally symmetric, hyperchaotic regime.

We present three-dimensional realizations of a class of lattice Boltzmann models introduced recently by the authors [I. V. Karlin, F. Bösch, and S. S. Chikatamarla, Phys. Rev. E 90, 031302(R) (2014)] and review the role of the entropic stabilizer. Both coarse- and fine-grid simulations are addressed for the Kida vortex flow benchmark. We show that the outstanding numerical stability and performance is independent of a particular choice of the moment representation for high-Reynolds-number flows. We report accurate results for low-order moments for homogeneous isotropic decaying turbulence and second-order grid convergence for most assessed statistical quantities. It is demonstrated that all the three-dimensional lattice Boltzmann realizations considered herein converge to the familiar lattice Bhatnagar-Gross-Krook model when the resolution is increased. Moreover, thanks to the dynamic nature of the entropic stabilizer, the present model features less compressibility effects and maintains correct energy and enstrophy dissipation. The explicit and efficient nature of the present lattice Boltzmann method renders it a promising candidate for both engineering and scientific purposes for highly turbulent flows.

A shortcut to adiabaticity is a finite-time process that produces the same final state as would result from infinitely slow driving. We show that such shortcuts can be found for weak perturbations from linear response theory. With the help of phenomenological response functions, a simple expression for the excess work is found—quantifying the nonequilibrium excitations. For two specific examples, i.e., the quantum parametric oscillator and the spin 1/2 in a time-dependent magnetic field, we show that finite-time zeros of the excess work indicate the existence of shortcuts. Finally, we propose a degenerate family of protocols, which facilitates shortcuts to adiabaticity for specific and very short driving times.

Maximum-entropy (ME) inference of state probabilities using state-dependent constraints is popular in the study of complex systems. In stochastic systems, how state space topology and path-dependent constraints affect ME-inferred state probabilities remains unknown. To that end, we derive the transition probabilities and the stationary distribution of a maximum path entropy Markov process subject to state- and path-dependent constraints. A main finding is that the stationary distribution over states differs significantly from the Boltzmann distribution and reflects a competition between path multiplicity and imposed constraints. We illustrate our results with particle diffusion on a two-dimensional landscape. Connections with the path integral approach to diffusion are discussed.