It is known that diffusion provokes substantial changes in continuous absorbing phase transitions. Conversely, its effect on discontinuous transitions is much less understood. In order to shed light in this direction, we study the inclusion of diffusion in the simplest one-dimensional model with a discontinuous absorbing phase transition, namely, the long-range contact process ($\sigma$-CP). Particles interact as in the usual CP, but the transition rate depends on the length $\ell$ of inactive sites according to $1+a{\ell }^{-\sigma }$, where $a$ and $\sigma$ are control parameters. The inclusion of diffusion in this model has been investigated by numerical simulations and mean-field calculations. Results show that there exists three distinct regimes. For sufficiently low and large $\sigma$'s the transition is, respectively, always discontinuous or continuous, independently of the strength of the diffusion. On the other hand, in an intermediate range of $\sigma$'s, the diffusion causes a suppression of the phase coexistence leading to a continuous transition belonging to the directed percolation universality class.

We have shown that the critical temperature of a Bose-Einstein condensate to a normal phase transition of noninteracting bosons in cubic optical lattices has a linear dependence on the filling factor, especially at large densities. The condensed fraction exhibits a linear power law dependence on temperature in contrast to the case of ideal homogeneous Bose gases.

We simulate with hydrodynamics a suspension of active disks squirming through a Newtonian fluid. We explore numerically the full range of squirmer area fractions from dilute to close packed and show that “motility induced phase separation,” which was recently proposed to arise generically in active matter, and which has been seen in simulations of active Brownian disks, is strongly suppressed by hydrodynamic interactions. We give an argument for why this should be the case and support it with counterpart simulations of active Brownian disks in a parameter regime that provides a closer counterpart to hydrodynamic suspensions than in previous studies.

In this paper we study the process of transition from passive to excitable behavior due to interaction between nonlinear dynamical systems. We show that under certain conditions a passive unit may demonstrate qualitatively new excitable dynamics. We study the properties of an excitable medium constructed on the basis of the proposed transition. The effects are demonstrated with the realistic Luo-Rudy model. Application to the cardiac dynamics and functioning is discussed. The qualitative analytic and numerical description is also given for the phenomenological FitzHugh-Nagumo system.

Red blood cells (RBCs) are the major component of blood, and the flow of blood is dictated by that of RBCs. We employ vesicles, which consist of closed bilayer membranes enclosing a fluid, as a model system to study the behavior of RBCs under a confined Poiseuille flow. We extensively explore two main parameters: (i) the degree of confinement of vesicles within the channel and (ii) the flow strength. Rich and complex dynamics for vesicles are revealed, ranging from steady-state shapes (in the form of parachute and slipper shapes) to chaotic dynamics of shape. Chaos occurs through a cascade of multiple periodic oscillations of the vesicle shape. We summarize our results in a phase diagram in the parameter plane (degree of confinement and flow strength). This finding highlights the level of complexity of a flowing vesicle in the small Reynolds number where the flow is laminar in the absence of vesicles and can be rendered turbulent due to elasticity of vesicles.

We present an experimental investigation of the statistical properties of the position fluctuations of low-loss oscillators in nonequilibrium steady states. The oscillators are coupled to a heat bath, and a nonequilibrium steady state is produced by flowing a constant heat flux, setting a temperature difference across the oscillators. We investigated the distribution of the measurements of the square of the oscillator position and searched for signs of changes with respect to the equilibrium case. We found that, after normalization by the mean value, the second, third, and fourth standardized statistical moments are not modified by the underlying thermodynamic state. This differs from the behavior of the absolute, i.e., not normalized, second moment, which is strongly affected by temperature gradients and heat fluxes. We illustrate this with a numerical experiment in which we study via molecular dynamics the fluctuations of the length of a one-dimensional chain of identical particles interacting via anharmonic interparticle potentials, with the extremes thermostated at different temperatures: we use the variance of the length in correspondence to its first elastic mode of resonance to define an effective temperature which we observe to depart from the thermodynamic one in the nonequilibrium states. We investigate the effect of changing the interparticle potential and show that the qualitative behavior of the nonequilibrium excess is unchanged. Our numerical results are consistent with the chain length being Gaussian distributed in the nonequilibrium states. Our experimental investigation reveals that the position variance is the only, and crucially easily accessible, observable for distinguishing equilibrium from nonequilibrium steady states. The consequences of this fact for the design of interferometric gravitational wave detectors are discussed.

We study the two-species symbiotic contact process, recently proposed by de Oliveira, Santos, and Dickman [Phys. Rev. E 86, 011121 (2012)]. In this model, each site of a lattice may be vacant or host single individuals of species $A$ and/or $B$. Individuals at sites with both species present interact in a symbiotic manner, having a reduced death rate $\mu <1$. Otherwise, the dynamics follows the rules of the basic contact process, with individuals reproducing to vacant neighbor sites at rate $\lambda$ and dying at a rate of unity. We determine the full phase diagram in the $\lambda \text{−}\mu$ plane in one and two dimensions by means of exact numerical quasistationary distributions, cluster approximations, and Monte Carlo simulations. We also study the effects of asymmetric creation rates and diffusion of individuals. In two dimensions, for sufficiently strong symbiosis (i.e., small $\mu$), the absorbing-state phase transition becomes discontinuous for diffusion rates $D$ within a certain range. We report preliminary results on the critical surface and tricritical line in the $\lambda \text{−}\mu \text{−}D$ space. Our results raise the possibility that strongly symbiotic associations of mobile species may be vulnerable to sudden extinction under increasingly adverse conditions.

We investigate the emergence of self-sustained oscillations in excitable Erdös-Rényi random networks (EERRNs). Interestingly, periodical self-sustained oscillations have been found at a moderate connection probability $P$. For smaller or larger $P$, the system evolves into a homogeneous rest state with distinct mechanisms. One-dimensional Winfree loops are discovered as the sources to maintain the oscillations. Moreover, by analyzing these oscillation sources, we propose two criteria to explain the spatiotemporal dynamics obtained in EERRNs. Finally, the two critical connection probabilities for which self-sustained oscillations can emerge are approximately predicted based on these two criteria.

In this paper we study the susceptible-infected-susceptible epidemic dynamics, considering a specialized setting where popular places (termed passive entities) are visited by agents (termed active entities). We consider two types of spreading dynamics: direct spreading, where the active entities infect each other while visiting the passive entities, and indirect spreading, where the passive entities act as carriers and the infection is spread via them. We investigate in particular the effect of selection strategy, i.e., the way passive entities are chosen, in the spread of epidemics. We introduce a mathematical framework to study the effect of an arbitrary selection strategy and derive formulas for prevalence, extinction probabilities, and epidemic thresholds for both indirect and direct spreading. We also obtain a very simple relationship between the extinction probability and the prevalence. We pay special attention to preferential selection and derive exact formulas. The analysis reveals that an increase in the diversity in the selection process lowers the epidemic thresholds. Comparing the direct and indirect spreading, we identify regions in the parameter space where the prevalence of the indirect spreading is higher than the direct one.

We introduce a formalism of fractional diffusion on networks based on a fractional Laplacian matrix that can be constructed directly from the eigenvalues and eigenvectors of the Laplacian matrix. This fractional approach allows random walks with long-range dynamics providing a general framework for anomalous diffusion and navigation, and inducing dynamically the small-world property on any network. We obtained exact results for the stationary probability distribution, the average fractional return probability, and a global time, showing that the efficiency to navigate the network is greater if we use a fractional random walk in comparison to a normal random walk. For the case of a ring, we obtain exact analytical results showing that the fractional transition and return probabilities follow a long-range power-law decay, leading to the emergence of Lévy flights on networks. Our general fractional diffusion formalism applies to regular, random, and complex networks and can be implemented from the spectral properties of the Laplacian matrix, providing an important tool to analyze anomalous diffusion on networks.

We consider the nonlinear Dirac equation in 1 + 1 dimension with scalar-scalar self interaction $\frac{g^2}{\kappa +1}{\left(\overline{\Psi }\Psi \right)}^{\kappa +1}$ and with mass $m$. Using the exact analytic form for rest frame solitary waves of the form $\Psi (x,t)=\psi \left(x\right)e^{-i\omega t}$ for arbitrary $\kappa$, we discuss the validity of various approaches to understanding stability that were successful for the nonlinear Schrödinger equation. In particular we study the validity of a version of Derrick's theorem and the criterion of Bogolubsky as well as the Vakhitov-Kolokolov criterion, and find that these criteria yield inconsistent results. Therefore, we study the stability by numerical simulations using a recently developed fourth-order operator splitting integration method. For different ranges of $\kappa$ we map out the stability regimes in $\omega$. We find that all stable nonlinear Dirac solitary waves have a one-hump profile, but not all one-hump waves are stable, while all waves with two humps are unstable. We also find that the time $t_c$, it takes for the instability to set in, is an exponentially increasing function of $\omega$ and $t_c$ decreases monotonically with increasing $\kappa$.

We explore a wide variety of patterns of closed surfaces that minimize the elastic bending energy with fixed surface area and volume. To avoid complicated discretization and numerical instabilities for sharp surfaces, we reformulate the underlying constrained minimization problem by constructing phase-field functionals of bending energy with penalty terms for the constraints and develop stable numerical methods to relax these functionals. We report our extensive computational results with different initial surfaces. These results are discussed in terms of the reduced volume and are compared with the known results obtained using the sharp-interface approach. Finally, we discuss the implications of our numerical findings.

We investigate the dynamics of single microparticles immersed in water that are driven out of equilibrium in the presence of an additional external colored noise. As a case study, we trap a single polystyrene particle in water with optical tweezers and apply an external electric field with flat spectrum but a finite bandwidth of the order of kHz. The intensity of the external noise controls the amplitude of the fluctuations of the position of the particle and therefore of its effective temperature. Here we show, in two different nonequilibrium experiments, that the fluctuations of the work done on the particle obey the Crooks fluctuation theorem at the equilibrium effective temperature, given that the sampling frequency and the noise cutoff frequency are properly chosen.

Time-integrated quantities such as work and heat increase incessantly in time during nonequilibrium processes near steady states. In the long-time limit, the average values of work and heat become asymptotically equivalent to each other, since they only differ by a finite energy change in average. However, the fluctuation theorem (FT) for the heat is found not to hold with the equilibrium initial ensemble, while the FT for the work holds. This reveals an intriguing effect of everlasting initial memory stored in rare events. We revisit the problem of a Brownian particle in a harmonic potential dragged with a constant velocity, which is in contact with a thermal reservoir. The heat and work fluctuations are investigated with initial Boltzmann ensembles at temperatures generally different from the reservoir temperature. We find that, in the infinite-time limit, the FT for the work is fully recovered for arbitrary initial temperatures, while the heat fluctuations significantly deviate from the FT characteristics except for the infinite initial-temperature limit (a uniform initial ensemble). Furthermore, we succeed in calculating finite-time corrections to the heat and work distributions analytically, using the modified saddle point integral method recently developed by us. Interestingly, we find noncommutativity between the infinite-time limit and the infinite-initial-temperature limit for the probability distribution function (PDF) of the heat.

We deduce a nonlinear and inhomogeneous Fokker-Planck equation within a generalized Stratonovich, or stochastic $\alpha$, prescription ($\alpha =0$, $1/2$, and 1, respectively, correspond to the Itô, Stratonovich and anti-Itô prescriptions). We obtain its stationary state $p_{st}\left(x\right)$ for a class of constitutive relations between drift and diffusion and show that it has a $q$-exponential form, $p_{st}\left(x\right)=N_q{[1-(1-q)\beta V\left(x\right)]}^{1/(1-q)}$, with an index $q$ which does not depend on $\alpha$ in the presence of any nonvanishing nonlinearity. This is in contrast with the linear case, for which the index $q$ is $\alpha$ dependent.

While density functional theory with integral equations techniques are very efficient tools in the numerical analysis of complex fluids, analytical insight into the phenomenon of effective interactions is still limited. In this paper, we propose a theory of binary systems that results in a relatively simple analytical expression combining arbitrary microscopic potentials into effective interaction. The derivation is based on translating a many-particle Hamiltonian including particle-depletant and depletant-depletant interactions into the occupation field language, which turns the partition function into multiple Gaussian integrals, regardless of what microscopic potentials are chosen. As a result, we calculate the effective Hamiltonian and discuss when our formula is a dominant contribution to the effective interactions. Our theory allows us to analytically reproduce several important characteristics of systems under scrutiny. In particular, we analyze the following: the effective attraction as a demixing factor in the binary systems of Gaussian particles, the screening of charged spheres by ions, which proves equivalent to Derjaguin-Landau-Verwey-Overbeek (DLVO) theory, effective interactions in the binary mixtures of Yukawa particles, and the system of particles consisting of both a repulsive core and an attractive/repulsive Yukawa interaction tail. For this last case, we reproduce the “attraction-through-repulsion” and “repulsion-through-attraction” effects previously observed in simulations.

Numerous real-world networks have been observed to interact with each other, resulting in interconnected networks that exhibit diverse, nontrivial behavior with dynamical processes. Here we investigate epidemic spreading on interconnected networks at the level of metapopulation. Through a mean-field approximation for a metapopulation model, we find that both the interaction network topology and the mobility probabilities between subnetworks jointly influence the epidemic spread. Depending on the interaction between subnetworks, proper controls of mobility can efficiently mitigate epidemics, whereas an extremely biased mobility to one subnetwork will typically cause a severe outbreak and promote the epidemic spreading. Our analysis provides a basic framework for better understanding of epidemic behavior in related transportation systems as well as for better control of epidemics by guiding human mobility patterns.

Thermally induced mechanical loading has been shown to have significant effects on micro- and nano-objects immersed in a gas with a nonuniform temperature field. While the majority of existing studies and related applications focus on forces, we investigate the torque, and thus the rotational motion, produced by such a mechanism. Our study has found that a torque can be induced if the configuration of the system is asymmetric. In addition, both the magnitude and the direction of the torque depend highly on the system configuration, indicating the possibility of manipulating the rotational motion via geometrical design. Based on this feature, two types of rotational micromotor that are of practical importance, namely pendulum motor and unidirectional motor, are designed. The magnitude of the torque at $\mathrm{Kn}=0.5$ can reach to around $2\mathrm{nN}\times \mu \mathrm{m}$ for a rectangular microbeam with a length of $100\mu \mathrm{m}$.

We investigate hydrodynamic interaction between two neutrally buoyant circular drops in a confined shear flow based on a computational fluid dynamics simulation using the volume-of-fluid method. The rheological behaviors of interactive drops and the flow regimes are explored with a focus on elucidation of underlying physical mechanisms. We find that two types of drop behaviors during interaction occur, including passing-over motion and reversing motion, which are governed by the competition between the drag of passing flow and the entrainment of reversing flow in matrix fluid. With the increasing confinement, the drop behavior transits from the passing-over motion to reversing motion, because the entrainment of the reversing-flow matrix fluid turns to play the dominant role. The drag of the ambient passing flow is increased by enlarging the initial lateral separation due to the departure of the drop from the reversing flow in matrix fluid, resulting in the emergence of passing-over motion. In particular, a corresponding phase diagram is plotted to quantitatively illustrate the dependence of drop morphologies during interaction on confinement and initial lateral separation.