We investigate the defocusing coupled nonlinear Schrödinger equations from a $3\times 3$ Lax pair. The Darboux transformations with the nonzero plane-wave solutions are presented to derive the newly localized wave solutions including dark-dark and bright-dark solitons, breather-breather solutions, and different types of new vector rogue wave solutions, as well as interactions between distinct types of localized wave solutions. Moreover, we analyze these solutions by means of parameters modulation. Finally, the perturbed wave propagations of some obtained solutions are explored by means of systematic simulations, which demonstrates that nearly stable and strongly unstable solutions. Our research results could constitute a significant contribution to explore the distinct nonlinear waves (e.g., dark solitons, breather solutions, and rogue wave solutions) dynamics of the coupled system in related fields such as nonlinear optics, plasma physics, oceanography, and Bose-Einstein condensates.

We perform experiments on a nonautonomous Mercury beating heart system, which is forced to pulsate using an external square wave potential. At suitable frequencies and volumes, the drop exhibits pulsation with polygonal shapes having $n$ corners. We find the scaling dependence of the forcing frequency ${\nu }_n$ on the volume $V$ of the drop and establish the relationship ${\nu }_n\propto \frac{n}{\sqrt{V}}$. It is shown that the geometrical shape of substrate is important for obtaining closer match to these scaling relationships. Furthermore, we study synchronization of two nonidentical drops driven by the same frequency and establish that synchrony happens when the relationship $n_2/n_1=\sqrt{V_2/V_1}$ is satisfied.

We present theoretical and experimental studies on pattern formation with bistable dynamical units coupled in a star network configuration. By applying a localized perturbation to the central or the peripheral elements, we demonstrate the subsequent spreading, pinning, or retraction of the activations; such analysis enables the characterization of the formation of stationary patterns of localized activity. The results are interpreted with a theoretical analysis of a simplified bistable reaction-diffusion model. Weak coupling results in trivial pinned states where the activation cannot propagate. At strong coupling, a uniform state is expected with active or inactive elements at small or large degree networks, respectively. A nontrivial stationary spatial pattern, corresponding to an activation pinning, is predicted to occur at an intermediate number of peripheral elements and at intermediate coupling strengths, where the central activation of the network is pinned, but the peripheral activation propagates toward the center. The results are confirmed in experiments with star networks of bistable electrochemical reactions. The experiments confirm the existence of the stationary spatial patterns and the dependence of coupling strength on the number of peripheral elements for transitions between pinned and retreating or spreading fronts in forced network configurations (where the central or periphery elements are forced to maintain their states).

In a recent work, mean-field analysis and computer simulations were employed to analyze critical self-organization in networks of excitable cellular automata where randomly chosen synapses in the network were depressed after each spike (the so-called annealed dynamics). Calculations agree with simulations of the annealed version, showing that the nominal branching ratio $\sigma$ converges to unity in the thermodynamic limit, as expected of a self-organized critical system. However, the question remains whether the same results apply to the biological case where only the synapses of firing neurons are depressed (the so-called quenched dynamics). We show that simulations of the quenched model yield significant deviations from $\sigma =1$ due to spatial correlations. However, the model is shown to be critical, as the largest eigenvalue of the synaptic matrix approaches unity in the thermodynamic limit, that is, ${\lambda }_c=1$. We also study the finite size effects near the critical state as a function of the parameters of the synaptic dynamics.

We generalize the original majority-vote model by incorporating inertia into the microscopic dynamics of the spin flipping, where the spin-flip probability of any individual depends not only on the states of its neighbors, but also on its own state. Surprisingly, the order-disorder phase transition is changed from a usual continuous or second-order type to a discontinuous or first-order one when the inertia is above an appropriate level. A central feature of such an explosive transition is a strong hysteresis behavior as noise intensity goes forward and backward. Within the hysteresis region, a disordered phase and two symmetric ordered phases are coexisting and transition rates between these phases are numerically calculated by a rare-event sampling method. A mean-field theory is developed to analytically reveal the property of this phase transition.

The basic principles of self-organization of one-component charged particles, confined in disk and circular parabolic potentials, are proposed. A system of equations is derived, which allows us to determine equilibrium configurations for an arbitrary, but finite, number of charged particles that are distributed over several rings. Our approach reduces significantly the computational effort in minimizing the energy of equilibrium configurations and demonstrates a remarkable agreement with the values provided by molecular dynamics calculations. With the increase of particle number $n>180$ we find a steady formation of a centered hexagonal lattice that smoothly transforms to valence circular rings in the ground-state configurations for both potentials.

Foams made of complex fluids such as particle suspensions have a great potential for the development of advanced aerated materials. In this paper, we study the rheological behavior of liquid foams loaded with granular suspensions. We focus on the effect of small particles, i.e., particle-to-bubble size ratio smaller than 0.1, and we measure the complex modulus as a function of particle size and particle volume fraction. With respect to previous work, the results highlight a new elastic regime characterized by unequaled modulus values as well as independence of size ratio. A careful investigation of the material microstructure reveals that particles organize through the network between the gas bubbles and form a granular skeleton structure with tightly packed particles. The latter is proven to be responsible for the reported new elastic regime. Rheological probing performed by strain sweep reveals a two-step yielding of the material: The first one occurs at small strain and is clearly attributed to yielding of the granular skeleton; the second one corresponds to the yielding of the bubble assembly, as observed for particle-free foams. Moreover, the elastic modulus measured at small strain is quantitatively described by models for solid foams in assuming that the granular skeleton possesses a bulk elastic modulus of order 100 kPa. Additional rheology experiments performed on the bulk granular material indicate that this surprisingly high value can be understood as soon as the magnitude of the confinement pressure exerted by foam bubbles on packed grains is considered.

We consider the electrophoresis of a charged colloid for a generalized situation in which the particle is considered to be polarizable and the surface exhibits hydrophobicity. The dielectric polarization of the particle creates a nonlinear dependence of the electrophoretic velocity on the applied electric field, and the core hydrophobicity amplifies the fluid convection in the Debye layer. Thus, a linear analysis is no longer applicable for this situation. The present analysis is based on the numerical solution of the nonlinear electrokinetic equations based on the Navier-Stokes-Nernst-Planck-Poisson equations coupled with the Laplace equation for the electric field within the dielectric particle. The hydrophobicity of the particle may influence its electric polarization by enhancing the convective transport of ions. The nonlinear effects, such as double-layer polarization and relaxation, are also influenced by the hydrophobicity of the particle surface. The present results compare well for a lower range of the applied electric field and surface charge density with the existing results for a perfectly dielectric particle with a hydrophobic surface based on the first-order perturbation analysis due to Khair and Squires [Phys. Fluids 21, 042001 (2009)]. Dielectric polarization creates a reduction in particle electrophoretic velocity, and its impact is strong for a moderate range of Debye length. A quantitative measure of the nonlinear effects is demonstrated by comparing the electrophoretic velocity with an existing linear model.

We numerically examine the dynamic phases and pattern formation of two-dimensional monodisperse repulsive disks driven over random quenched disorder. We show that there is a series of distinct dynamic regimes as a function of increasing drive, including a clogged or pile-up phase near depinning, a homogeneous disordered flow state, and a dynamically phase separated regime consisting of high-density crystalline regions surrounded by a low density of disordered disks. At the highest drives the disks arrange into one-dimensional moving chains. The phase separated regime has parallels with the phase separation observed in active matter systems, but arises from a distinct mechanism consisting of the combination of nonequilibrium fluctuations with density-dependent mobility. We discuss the pronounced differences between this system and previous studies of driven particles with longer-range repulsive interactions moving over random substrates, such as superconducting vortices or electron crystals, where dynamical phase separation and distinct one-dimensional moving chains are not observed. Our results should be generic to a broad class of systems in which the particle-particle interactions are short ranged, such as sterically interacting colloids or Yukawa particles with strong screening driven over random pinning arrays, superconducting vortices in the limit of small penetration depths, or quasi-two-dimensional granular matter flowing over rough landscapes.

The macroscopic laws controlling the advection and diffusion of solute at the scale of the porous continuum are derived in a general manner that does not place limitations on the geometry and time evolution of the pore space. Special focus is given to the definition and symmetry of the dispersion tensor that is controlling how a solute plume spreads out. We show that the dispersion tensor is not symmetric and that the asymmetry derives from the advective derivative in the pore-scale advection-diffusion equation. When flow is spatially variable across a voxel, such as in the presence of a permeability gradient, the amount of asymmetry can be large. As first shown by Auriault [J.-L. Auriault et al. Transp. Porous Med. 85, 771 (2010)] in the limit of low Péclet number, we show that at any Péclet number, the dispersion tensor $D_{ij}$ satisfies the flow-reversal symmetry $D_{ij}(+\mathbf{q})=D_{ji}(-\mathbf{q})$ where $\mathbf{q}$ is the mean flow in the voxel under analysis; however, Reynold's number must be sufficiently small that the flow is reversible when the force driving the flow changes sign. We also demonstrate these symmetries using lattice-Boltzmann simulations and discuss some subtle aspects of how to measure the dispersion tensor numerically. In particular, the numerical experiments demonstrate that the off-diagonal components of the dispersion tensor are antisymmetric which is consistent with the analytical dependence on the average flow gradients that we propose for these off-diagonal components.

We use classical molecular dynamics (MD) to study electron-ion temperature equilibration in two-component plasmas in regimes for which the presence of coupled collective modes has been predicted to substantively reduce the equilibration rate. Guided by previous kinetic theory work, we examine hydrogen plasmas at a density of $n={10}^{26}{\mathrm{cm}}^{-3}, T_i={10}^5\mathrm{K}$, and ${10}^7\mathrm{K}<T_e<{10}^9\mathrm{K}$. The nonequilibrium classical MD simulations are performed with interparticle interactions modeled by quantum statistical potentials (QSPs). Our MD results indicate (i) a large effect from time-varying potential energy, which we quantify by appealing to an adiabatic two-temperature equation of state, and (ii) a notable deviation in the energy equilibration rate when compared to calculations from classical Lenard-Balescu theory including the QSPs. In particular, it is shown that the energy equilibration rates from MD are more similar to those of the theory when coupled modes are neglected. We suggest possible reasons for this surprising result and propose directions of further research along these lines.

The percolation threshold of fracture networks is investigated by extensive direct numerical simulations. The fractures are randomly located and oriented in three-dimensional space. A very wide range of regular, irregular, and random fracture shapes is considered, in monodisperse or polydisperse networks containing fractures with different shapes and/or sizes. The results are rationalized in terms of a dimensionless density. A simple model involving a new shape factor is proposed, which accounts very efficiently for the influence of the fracture shape. It applies with very good accuracy in monodisperse or moderately polydisperse networks, and provides a good first estimation in other situations. A polydispersity index is shown to control the need for a correction, and the corrective term is modelled for the investigated size distributions.

I study the heat transport properties of a chain of coupled quantum harmonic oscillators in contact at its ends with two heat reservoirs at distinct temperatures. My approach is based on the use of an evolution equation for the density operator which is a canonical quantization of the classical Fokker-Planck-Kramers equation. I set up the evolution equation for the covariances and obtain the stationary covariances at the stationary states from which I determine the thermal conductance in closed form when the interparticle interaction is small. The conductance is finite in the thermodynamic limit implying an infinite thermal conductivity.

In this paper, we investigate and develop an alternative approach to the numerical analysis and characterization of random fluctuations with the heavy-tailed probability distribution function (PDF), such as turbulent heat flow and solar flare fluctuations. We identify the heavy-tailed random fluctuations based on the scaling properties of the tail exponent of the PDF, power-law growth of $q\mathrm{th}$ order correlation function, and the self-similar properties of the contour lines in two-dimensional random fields. Moreover, this work leads to a substitution for the fractional Edwards-Wilkinson (EW) equation that works in the presence of $\mu$-stable Lévy noise. Our proposed model explains the configuration dynamics of the systems with heavy-tailed correlated random fluctuations. We also present an alternative solution to the fractional EW equation in the presence of $\mu$-stable Lévy noise in the steady state, which is implemented numerically, using the $\mu$-stable fractional Lévy motion. Based on the analysis of the self-similar properties of contour loops, we numerically show that the scaling properties of contour loop ensembles can qualitatively and quantitatively distinguish non-Gaussian random fields from Gaussian random fluctuations.

Two-dimensional electrophoresis was used to analyze the mobility of DNA fragments in micellar gels of pluronic F127 $\left(\mathrm{E}{\mathrm{O}}_{100}\mathrm{P}{\mathrm{O}}_{70}\mathrm{E}{\mathrm{O}}_{100}\right)$ and pluronic P123 $\left(\mathrm{E}{\mathrm{O}}_{20}\mathrm{P}{\mathrm{O}}_{70}\mathrm{E}{\mathrm{O}}_{20}\right)$. The 20–3500 base pair DNA fragments were separated by size first in agarose gels, and then in pluronic gels at room temperature. In agarose gels, the DNA mobility decreases monotonically with increasing DNA length. In pluronic gels, however, the mobility varies nonmonotonically according to fragment lengths that are strongly correlated with the diameter of the spherical micelles. Brownian dynamics (BD) simulations with short-ranged intra-DNA hydrodynamic interactions were performed to numerically calculate the length-dependent mobility in pluronic lattices. The rising and falling trends, as well as the oscillations of mobility, were captured by the coarse-grained BD simulations. Molecular dynamics simulations in pluronic F127, with explicitly modeled micelle coronas, justified that the hydrodynamic interactions mediated by the complex fluid of hydrated poly(ethylene oxide) are a possible reason for the initial rise of mobility with DNA length.

The high temperature expansion (HTE) of the specific heat of a spin system fails at low temperatures, even if it is combined with a Padé approximation. On the other hand, we often have information about the low-temperature asymptotics (LTA) of the system. Interpolation methods combine both kind of information, HTE and LTA, in order to obtain an approximation of the specific heat that holds for the whole temperature range. Here we revisit the entropy method that has been previously published and propose two variants that better cope with problems of the entropy method for gapped systems. We compare all three methods applied to the antiferromagnetic Haldane spin-one chain and especially apply the second variant, called log Z method, to the cuboctahedron for different spin quantum numbers. In particular, we demonstrate that the interpolation method is able to detect an extra low-temperature maximum in the specific heat that may appear if a separation of two energy scales is present in the considered system. Finally, we illustrate how interpolation also works for classical spin systems.

In this paper, we consider entanglement in a system out of equilibrium, adopting the viewpoint given by the formalism of superstatistics. Such an approach yields a good effective description for a system in a slowly fluctuating environment within a weak interaction between the system and the environment. For this purpose, we introduce an alternative version of the formalism within a quantum mechanical picture and use it to study entanglement in the Heisenberg $XY$ model, subject to temperature fluctuations. We consider both isotropic and anisotropic cases and explore the effect of different temperature fluctuations (${\chi }^2$, log-normal, and $F$ distributions). Our results suggest that particular fluctuations may enhance entanglement and prevent it from vanishing at higher temperatures than those predicted for the same system at thermal equilibrium.

We investigate the short-time universal behavior of the two-dimensional Ashkin-Teller model at the Baxter line by performing time-dependent Monte Carlo simulations. First, as preparatory results, we obtain the critical parameters by searching the optimal power-law decay of the magnetization. Thus, the dynamic critical exponents ${\theta }_m$ and ${\theta }_p$, related to the magnetic and electric order parameters, as well as the persistence exponent ${\theta }_g$, are estimated using heat-bath Monte Carlo simulations. In addition, we estimate the dynamic exponent $z$ and the static critical exponents $\beta$ and $\nu$ for both order parameters. We propose a refined method to estimate the static exponents that considers two different averages: one that combines an internal average using several seeds with another, which is taken over temporal variations in the power laws. Moreover, we also performed the bootstrapping method for a complementary analysis. Our results show that the ratio $\beta /\nu$ exhibits universal behavior along the critical line corroborating the conjecture for both magnetization and polarization.

We investigate the phase transitions in a coupled system of Ising spins and a fluctuating network. Each spin interacts with $q$ neighbors through links of the rewiring network. The Ising spins and the network are in thermal contact with the heat baths at temperatures $T_S$ and $T_L$, respectively, so the whole system is driven out of equilibrium for $T_S\ne T_L$. The model is a generalization of the $q$-neighbor Ising model [A. Jędrzejewski et al., Phys. Rev. E 92, 052105 (2015)], which corresponds to the limiting case of $T_L=\infty$. Despite the mean-field nature of the interaction, the $q$-neighbor Ising model was shown to display a discontinuous phase transition for $q\ge 4$. Setting up the rate equations for the magnetization and the energy density, we obtain the phase diagram in the $T_S-T_L$ parameter space. The phase diagram consists of a ferromagnetic phase and a paramagnetic phase. The two phases are separated by a continuous phase transition belonging to the mean-field universality class or by a discontinuous phase transition with an intervening coexistence phase. The equilibrium system with $T_S=T_L$ falls into the former case while the $q$-neighbor Ising model falls into the latter case. At the tricritical point, the system exhibits the mean-field tricritical behavior. Our model demonstrates a possibility that a continuous phase transition turns into a discontinuous transition by a nonequilibrium driving. Heat flow induced by the temperature difference between two heat baths is also studied.