In this paper, a lattice Boltzmann model is developed to solve a class of the nonlinear wave equations. Through selecting equilibrium distribution function and an amending function properly, the governing evolution equation can be recovered correctly according to our proposed scheme, in which the Chapman-Enskog expansion is employed. We validate the algorithm on some problems where analytic solutions are available, including the second-order telegraph equation, the nonlinear Klein-Gordon equation, and the damped, driven sine-Gordon equation. It is found that the numerical results agree well with the analytic solutions, which indicates that the present algorithm is very effective and can be used to solve more general nonlinear problems.

We report a hybrid numerical method for the solution of the Model H fluctuating hydrodynamic equations for binary mixtures. The momentum conservation equations with Landau-Lifshitz stresses are solved using the fluctuating lattice Boltzmann equation while the order parameter conservation equation with Langevin fluxes is solved using stochastic method of lines. Two methods, based on finite difference and finite volume, are proposed for spatial discretization of the order parameter equation. Special care is taken to ensure that the fluctuation-dissipation theorem is maintained at the lattice level in both cases. The methods are benchmarked by comparing static and dynamic correlations and excellent agreement is found between analytical and numerical results. The Galilean invariance of the model is tested and found to be satisfactory. Thermally induced capillary fluctuations of the interface are captured accurately, indicating that the model can be used to study nonlinear fluctuations.

Chaotically spiking attractors in semiconductor lasers with optoelectronic feedback have been recently observed to be the result of canard phenomena in three-dimensional phase space (incomplete homoclinic scenarios). Since light-emitting diodes display the same dynamics and are much more easily controllable, we use one of these systems to complete the attractor analysis demonstrating experimentally and theoretically the occurrence of complex sequences of periodic mixed-mode oscillations. In particular, we investigate the transition between periodic and chaotic mixed-mode states and analyze the effects of the unavoidable experimental noise on these transitions.

The distribution of the maximal relative height (MRH) of self-affine one-dimensional elastic interfaces in a random potential is studied. We analyze the ground-state configuration at zero driving force, and the critical configuration exactly at the depinning threshold, both for the random-manifold and random-periodic universality classes. These configurations are sampled by exact numerical methods, and their MRH distributions are compared with those with the same roughness exponent and boundary conditions, but produced by independent Fourier modes with normally distributed amplitudes. Using Pickands’ theorem we derive an exact analytical description for the right tail of the latter. After properly rescaling the MRH distributions we find that corrections from the Gaussian independent modes approximation are, in general, small, as previously found for the average width distribution of depinning configurations. In the large size limit all corrections are finite except for the ground state in the random-periodic class whose MRH distribution becomes, for periodic boundary conditions, indistinguishable from the Airy distribution. We find that the MRH distributions are, in general, sensitive to changes of boundary conditions.

We extend recent analyses of stochastic effects in game dynamical learning to cases of multiplayer games and to games defined on networked structures. By means of an expansion in the noise strength we consider the weak-noise limit and present an analytical computation of spectral properties of fluctuations in multiplayer public goods games. This extends existing work on two-player games. In particular we show that coherent cycles may emerge driven by noise in the adaptation dynamics. These phenomena are not too dissimilar from cyclic strategy switching observed in experiments of behavioral game theory.

The coefficient of restitution of a spherical particle in contact with a flat plate is investigated as a function of the impact velocity. As an experimental observation we notice nontrivial (non-Gaussian) fluctuations of the measured values. For a fixed impact velocity, the probability density of the coefficient of restitution, $p\left(\varepsilon \right)$, is formed by two exponential functions (one increasing, one decreasing) of different slope. This behavior may be explained by a certain roughness of the particle which leads to energy transfer between the linear and rotational degrees of freedom.

We present an investigation of the ordering in a nematic liquid-crystal film confined between two cylindrical surfaces with antagonistic (radial and planar) anchoring alignments. A Monte Carlo study of a Lebwohl-Lasher model with suitable boundary conditions has been performed to calculate the ordering and the molecular organization for different film thicknesses. The simulation results are compared with some theoretical predictions obtained with the elastic continuum approach. The agreement between theory and simulation is improved as the thickness decreases.

This paper studies the formation of the large dark patterns, known as parr marks, that form on the Amago trout as it grows from the early larval stages to adulthood. The Amago trout, known as Oncorhynchus masou ishikawa, exhibits stripes during the early stages of development that in turn evolve (through reorientation and peak insertion) to form zigzag spot patterns as the fish grows to adulthood. By considering a standard representation of the Turing model for biological self-organization via interacting and diffusing morphogens, we illustrate that a diffusively driven instability can generate transient patterns consistent with those experimentally observed during the process of parr-mark formation in the early development of the Amago trout. Surface evolution is modeled through an experimentally driven growth function. Our studies conclude that the surface evolution profile, the surface geometry, and the curvature are key factors that play a pivotal role in reaction-diffusion systems in a study motivated by observations of Amago trout parr-mark pattern formation.

We develop a minimal model for the stochastic dynamics of microorganisms where individuals communicate via autochemotaxis. This means that microorganisms, such as bacteria, amoebae, or cells, follow the gradient of a chemical that they produce themselves to attract or repel each other. A microorganism is represented as a self-propelled particle or walker with constant speed while its velocity direction diffuses on the unit circle. We study the autochemotactic response of a single self-propelled walker whose dynamics is non-Markovian. We show that its long-time dynamics is always diffusive by deriving analytic expressions for its diffusion coefficient in the weak- and strong-coupling case. We confirm our findings by numerical simulations.

We study time series and the spontaneous miniature end-plate potentials (MEPPs) of mammals recorded at neuromuscular junctions using two different approaches: generalized thermostatistics and detrended fluctuation analysis (DFA). Classical concepts establish that the magnitude of these potentials is characterized by Gaussian statistics and that their intervals are randomly displayed. First we show that MEPP distributions adequately satisfy the $q$-Gaussian distributions that maximize the Tsallis entropy, indicating their nonextensive and nonequilibrium behavior. We then examine the intervals between the miniature potentials via DFA, where the profile of the intervals between events configures a deviation from the expected random behavior. Some possible physiological substrates for these findings are discussed.

We present numerically exact results for the quantum coherent energy transfer in the Fenna-Matthews-Olson molecular aggregate under realistic physiological conditions, including vibrational fluctuations of the protein and the pigments for an experimentally determined fluctuation spectrum. We find coherence times shorter than observed experimentally. Furthermore, we determine the energy transfer current and quantify its “quantumness” as the distance of the density matrix to the classical pointer states for the energy current operator. Most importantly, we find that the energy transfer happens through a “Schrödinger-cat-like” superposition of energy current pointer states.

This paper describes a systematic approach to suppressing cardiac alternans in simulated Purkinje fibers using localized current injections. We investigate the controllability and observability of the periodically paced Noble model for different locations of the recording and control electrodes. In particular, we show that the loss of controllability causes the failure of the control approach introduced by Echebarria and Karma [Chaos 12, 923 (2002)] for longer fiber lengths. Furthermore, we explain how the optimal locations for the recording and control electrodes and the timing of the feedback current can be selected, accounting for both linear and nonlinear effects, effectively doubling the length of fibers that can be controlled with previous methods.

We compute the distribution of relative velocities for a one-dimensional model of heavy particles suspended in a turbulent flow, quantifying the caustic contribution to the moments of relative velocities. The same principles determine the corresponding caustic contribution in $d$ spatial dimensions. The distribution of relative velocities $\Delta \mathbf{v}$ at small separations $R$ acquires the universal form $\rho (\Delta \mathbf{v},R)\sim R^{d-1}{\left|\Delta \mathbf{v}\right|}^{D_2-2d}$ for large (but not too large) values of $\left|\Delta \mathbf{v}\right|$. Here $D_2$ is the phase-space correlation dimension. Our conclusions are in excellent agreement with numerical simulations of particles suspended in a randomly mixing flow in two dimensions, and in quantitative agreement with published data on direct numerical simulations of particles in turbulent flows.

We present a simple cellular automaton (CA) model for two-lane roads explaining the physics of traffic breakdown, highway capacity, and synchronized flow. The model consists of the rules “acceleration,” “deceleration,” “randomization,” and “motion” of the Nagel-Schreckenberg CA model as well as “overacceleration through lane changing to the faster lane,” “comparison of vehicle gap with the synchronization gap,” and “speed adaptation within the synchronization gap” of Kerner's three-phase traffic theory. We show that these few rules of the CA model can appropriately simulate fundamental empirical features of traffic breakdown and highway capacity found in traffic data measured over years in different countries, like characteristics of synchronized flow, the existence of the spontaneous and induced breakdowns at the same bottleneck, and associated probabilistic features of traffic breakdown and highway capacity. Single-vehicle data derived in model simulations show that synchronized flow first occurs and then self-maintains due to a spatiotemporal competition between speed adaptation to a slower speed of the preceding vehicle and passing of this slower vehicle. We find that the application of simple dependences of randomization probability and synchronization gap on driving situation allows us to explain the physics of moving synchronized flow patterns and the pinch effect in synchronized flow as observed in real traffic data.

Evolutionary game theory on dynamical networks has received much attention. Most of the work has been focused on $2\times 2$ games such as prisoner's dilemma and snowdrift, with general $n\times n$ games seldom addressed. In particular, analytical methods are still lacking. Here we generalize the stochastic linking dynamics proposed by Wu, Zhou, Fu, Luo, Wang, and Traulsen [PLoS ONE 5, e11187 (2010)] to $n\times n$ games. We analytically obtain that the fast linking dynamics results in the replicator dynamics with a rescaled payoff matrix. In the rescaled matrix, intuitively, each entry is the product of the original entry and the average duration time of the corresponding link. This result is shown to be robust to a wide class of imitation processes. As applications, we show both analytically and numerically that the biodiversity, modeled as the stability of a zero-sum rock-paper-scissors game, cannot be altered by the fast linking dynamics. In addition, we show that the fast linking dynamics can stabilize tit-for-tat as an evolutionary stable strategy in the repeated prisoner's dilemma game provided the interaction between the identical strategies happens sufficiently often. Our method paves the way for an analytical study of the multiple-strategy coevolutionary dynamics.

We study the mixing and degeneracy of two unbound energy eigenstates (resonances) in a two coupled channel model of scattering and reactions. We derive the necessary and sufficient conditions for existence of an exceptional point in the extended spectrum of bound and resonance energy eigenvalues in this model and show that these are not the same as in the single channel case. When these conditions are satisfied, in the complex energy plane, the two simple resonance poles of the scattering matrix merge into one double pole at the exceptional point. In parameter space, the surface of the eigenenergies has a branch point of square root type and branch cuts in its real and imaginary parts that start at the exceptional point and extend in opposite directions. The rich phenomenology of crossings and anticrossings of energies and widths of the doublet of unbound states, as well as the changes of identity of the poles of the scattering matrix observed when one control parameter is varied while the other is kept constant, is fully explained in terms of sections of the eigenenergy surfaces.

Classical turbulent mixing paradigm—ingrained in scaling laws and closure models—is revisited in a variable-conductivity inhomogeneous medium. We perform direct numerical simulations to study the evolution of a passive scalar (temperature) field in a fluid with large conductivity gradients and investigate the behavior of scalar dissipation, conditional scalar dissipation, and velocity-to-scalar time scale ratio. Subject to the conditions of the investigation, it is found that these mixing characteristics become reasonably insensitive to conductivity after about one-third eddy turnover time. While the results support the classical paradigm, important preconditions and limitations are clearly identified.

The influence of an induced axial magnetic field on plasma dynamics and radiative characteristics of $Z$ pinches is investigated. An axial magnetic field was induced in a novel $Z$-pinch load: a double planar wire array with skewed wires (DPWAsk), which represents a planar wire array in an open magnetic configuration. The induced axial magnetic field suppressed magneto-Rayleigh-Taylor (MRT) instabilities (with $m$ $=$ 0 and $m$ $=$ 1 instability modes) in the $Z$-pinch plasma. The influence of the initial axial magnetic field on the structure of the plasma column at stagnation was manifested through the formation of a more uniform plasma column compared to a standard double planar wire array (DPWA) load [V. L. Kantsyrev et al., Phys. Plasmas 15, 030704 (2008)]. The DPWAsk load is characterized by suppression of MRT instabilities and by the formation of the sub-keV radiation pulse that occurs before the main x-ray peak. Gradients in plasma parameters along the cathode-anode gap were observed and analyzed for DPWAsk loads made from low atomic number $Z$ (Al) and mid-$Z$ (brass) wires.

We investigate buckling of soft elastic capsules under negative pressure or for reduced capsule volume. Based on nonlinear shell theory and the assumption of a hyperelastic capsule membrane, shape equations for axisymmetric and initially spherical capsules are derived and solved numerically. A rich bifurcation behavior is found, which is presented in terms of bifurcation diagrams. The energetically preferred stable configuration is deduced from a least-energy principle both for prescribed volume and prescribed pressure. We find that buckled shapes are energetically favorable already at smaller negative pressures and larger critical volumes than predicted by the classical buckling instability. By preventing self-intersection for strongly reduced volume, we obtain a complete picture of the buckling process and can follow the shape from the initial undeformed state through the buckling instability into the fully collapsed state. Interestingly, the sequences of bifurcations and stable capsule shapes differ for prescribed volume and prescribed pressure. In the buckled state, we find a relation between curvatures at the indentation rim and the bending modulus, which can be used to determine elastic moduli from experimental shape analysis.

Dynamics of a chain of interacting parity-time-invariant nonlinear dimers is investigated. A dimer is built as a pair of coupled elements with equal gain and loss. A relation between stationary soliton solutions of the model and solitons of the discrete nonlinear Schrödinger (DNLS) equation is demonstrated. Approximate solutions for solitons whose width is large in comparison to the lattice spacing are derived, using a continuum counterpart of the discrete equations. These solitons are mobile, featuring nearly elastic collisions. Stationary solutions for narrow solitons, which are immobile due to the pinning by the effective Peierls-Nabarro potential, are constructed numerically, starting from the anticontinuum limit. The solitons with the amplitude exceeding a certain critical value suffer an instability leading to blowup, which is a specific feature of the nonlinear parity-time-symmetric chain, making it dynamically different from DNLS lattices. A qualitative explanation of this feature is proposed. The instability threshold drops with the increase of the gain-loss coefficient, but it does not depend on the lattice coupling constant, nor on the soliton's velocity.