A trajectory-based formulation for fractional dynamics is presented and the trajectories are generated deterministically. In this theoretical framework, we derive a new class of estimators in terms of confluent hypergeometric function (${}_{1}F_1$) to represent the Riesz fractional derivative. Using this method, the simulation of free and confined Lévy flight are in excellent agreement with the exact numerical and analytical results. In addition, the barrier crossing in a bistable potential driven by Lévy noise of index $\alpha$ is investigated. In phase space, the behavior of trajectories reveal the feature of Lévy flight in a better perspective.

An unusual correlation function was conjectured by Campostrini et al. [Phys. Rev. E 91, 042123 (2015)] for the ground state of a transverse Ising chain with geometrical frustration. Later, we provided a rigorous proof for it and demonstrated its nonlocal nature based on an evaluation of a Toeplitz determinant in the thermodynamic limit [J. Stat. Mech. (2016) 113102]. In this paper, we further prove that all the low excited energy states forming the gapless kink phase share the same asymptotic correlation function with the ground state. As a consequence, the thermal correlation function almost remains constant at low temperatures if one assumes a canonical ensemble.

The turnstile provides us a useful tool to describe the flux in twist Hamiltonian systems. Thus, its determination allows us to find the areas where the trajectories flux through barriers. We show that the mechanism of the turnstile can increase the flux in nontwist Hamiltonian systems. A model which captures the essence of these systems is the standard nontwist map, introduced by del Castillo-Negrete and Morrison. For selected parameters of this map, we show that chaotic trajectories entering in resonances zones can be explained by turnstiles formed by a set of homoclinic points. We argue that for nontwist systems, if the heteroclinic points are sufficiently close, they can connect twin-islands chains. This provides us a scenario where the trajectories can cross the resonance zones and increase the flux. For these cases the escape basin boundaries are nontrivial, which demands the use of an appropriate characterization. We applied the uncertainty exponent and the entropies of the escape basin boundary in order to quantify the degree of unpredictability of the asymptotic trajectories.

We analyze the stability of the network's giant connected component under impact of adverse events, which we model through the link percolation. Specifically, we quantify the extent to which the largest connected component of a network consists of the same nodes, regardless of the specific set of deactivated links. Our results are intuitive in the case of single-layered systems: the presence of large degree nodes in a single-layered network ensures both its robustness and stability. In contrast, we find that interdependent networks that are robust to adverse events have unstable connected components. Our results bring novel insights to the design of resilient network topologies and the reinforcement of existing networked systems.

The influence of zealots on the noisy voter model is studied theoretically and numerically at the mean-field level. The noisy voter model is a modification of the voter model that includes a second mechanism for transitions between states: Apart from the original herding processes, voters may change their states because of an intrinsic noisy-in-origin source. By increasing the importance of the noise with respect to the herding, the system exhibits a finite-size phase transition from a quasiconsensus state, where most of the voters share the same opinion, to one with coexistence. Upon introducing some zealots, or voters with fixed opinion, the latter scenario may change significantly. We unveil new situations by carrying out a systematic numerical and analytical study of a fully connected network for voters, but allowing different voters to be directly influenced by different zealots. We show that this general system is equivalent to a system of voters without zealots, but with heterogeneous values of their parameters characterizing herding and noisy dynamics. We find excellent agreement between our analytical and numerical results. Noise and herding or zealotry acting together in the voter model yields a nontrivial mixture of the scenarios with the two mechanisms acting alone: It represents a situation where the global-local (noise-herding) competition is coupled to a symmetry breaking (zealots). In general, the zealotry enhances the effective noise of the system, which may destroy the original quasiconsensus state, and can introduce a bias towards the opinion of the majority of zealots, hence breaking the symmetry of the system and giving rise to new phases. In the most general case we find two different transitions: a discontinuous transition from an asymmetric bimodal phase to an extreme asymmetric phase and a second continuous transition from the extreme asymmetric phase to an asymmetric unimodal phase.

The formation of a functional vascular network depends on biological, chemical, and physical processes being extremely well coordinated. Among them, the mechanical properties of the extracellular matrix and cell adhesion are fundamental to achieve a functional network of endothelial cells, able to fully cover a required domain. By the use of a Cellular Potts Model and Finite Element Method it is shown that there exists a range of values of endothelial traction forces, cell-cell adhesion, and matrix rigidities where the network can spontaneously be formed, and its properties are characterized. We obtain the analytical relation that the minimum traction force required for cell network formation must obey. This minimum value for the traction force is approximately independent on the considered cell number and cell-cell adhesion. We quantify how these two parameters influence the morphology of the resulting networks (size and number of meshes).

We consider many-site mutation-recombination models of evolution with selection. We are looking for situations where the recombination increases the mean fitness of the population, and there is an optimal recombination rate. We found two fitness landscapes supporting such nonmonotonic behavior of the mean fitness versus the recombination rate. The first case is related to the evolution near the error threshold on a neutral-network-like fitness landscape, for moderate genome lengths and large population. The more realistic case is the second one, in which we consider the evolutionary dynamics of a finite population on a rugged fitness landscape (the smooth fitness landscape plus some random contributions to the fitness). We also give the solution to the horizontal gene transfer model in the case of asymmetric mutations. To obtain nonmonotonic behavior for both mutation and recombination, we need a specially designed (ideal) fitness landscape.

The motility of a cell can be triggered or inhibited not only by an applied force but also by a mechanically neutral force couple. This type of loading, represented by an applied stress and commonly interpreted as either squeezing or stretching, can originate from extrinsic interaction of a cell with its neighbors. To quantify the effect of applied stresses on cell motility we use an analytically transparent one-dimensional model accounting for active myosin contraction and induced actin turnover. We show that stretching can polarize static cells and initiate cell motility while squeezing can symmetrize and arrest moving cells. We show further that sufficiently strong squeezing can lead to the loss of cell integrity. The overall behavior of the system depends on the two dimensionless parameters characterizing internal driving (chemical activity) and external loading (applied stress). We construct a phase diagram in this parameter space distinguishing between static, motile, and collapsed states. The obtained results are relevant for the mechanical understanding of contact inhibition and the epithelial-to-mesenchymal transition.

We perform statistical mechanics calculations to analyze the global phase diagram of a fully connected version of a Maier-Saupe-Zwanzig lattice model with the inclusion of couplings to an elastic strain field. We point out the presence of uniaxial and biaxial nematic structures, depending on temperature $T$ and on the applied stress $\sigma$. Under uniaxial extensive tension, applied stress favors uniaxial orientation, and we obtain a first-order boundary along which there is a coexistence of two uniaxial paranematic phases, and which ends at a simple critical point. Under uniaxial compressive tension, stress favors biaxial orientation; for small values of the coupling parameters, the first-order boundary ends at a tricritical point, beyond which there is a continuous transition between a paranematic and a biaxially ordered structure. For some representative choices of the model parameters, we obtain a number of analytic results, including the location of critical and tricritical points and the line of stability of the biaxial phase.

We consider the motion of a gravity-driven flow down a vertical fiber subjected to a radial electric field. This flow exhibits rich dynamics including the formation of droplets, or beads, driven by a Rayleigh-Plateau mechanism modified by the presence of gravity as well as the Maxwell stress at the interface. A spatiotemporal stability analysis is performed to investigate the effect of electric field on the absolute-convective instability (AI-CI) characteristics. We performed a numerical simulation on the nonlinear evolution of the film to examine the transition from CI to AI regime. The numerical results are in excellent agreement with the spatiotemporal stability analysis. The blowup behavior of nonlinear simulation predicts the formation of touchdown singularity of the interface due to the effect of electric field. We try to connect the blowup behavior with the AI-CI characteristics. It is found that the singularities mainly occur in the AI regime. The results indicate that the film may have a tendency to form very sharp tips due to the enhancement of the absolute instability induced by the electric field. We perform a theoretical analysis to study the behaviors of the singularities. The results show that there exists a self-similarity between the temporal and spatial distances from the singularities.

This paper investigates the three-dimensional instabilities of the flow past a periodically heaving airfoil. By comparison with a pitching foil [Deng et al., Phys. Rev. E 92, 063013 (2015)], here we present distinctive characteristics for the heaving foil, particulary regarding its Floquet modes. By increasing the frequency (Sr), or equivalently decreasing the amplitude ($A_D$) along the marginal stability curve in the ($\mathrm{Sr},A_D$) phase space, the critical Floquet mode emerges sequentially as A, quasiperiodic (QP), and B. It is interesting to note that both modes A and B are synchronous with the base flow, in contrast to the quasiperiodic mode QP. To further investigate the instability across the marginal curve, we fix the frequency at $\mathrm{Sr}=0.187$, of which the critical Floquet mode is located in the synchronous regime, while varying $A_D$ around the critical point. We find that the dominant mode switches from mode A to mode B, while mode QP never becomes critical as we increase $A_D$. We note that mode S, a subharmonic mode, can also be unstable, which, however, is not physically realizable, because the magnitude of its Floquet multiplier is always smaller than that of mode B. We have also studied the influence of various Reynolds numbers at the same critical point on the marginal stability curve, with the results resembling that by varying the amplitude $A_D$.

A model to study the dynamics of colloidal particles in nonuniform electric fields is proposed. For an isolated sphere, the conditions and threshold for sustained (Quincke) rotation in a linear direct current (dc) field are determined. Particle dynamics becomes more complex with increasing electric field strength, changing from steady spinning around the particle center to time-dependent orbiting motion around the minimum field location. Pairs of particles exhibit intricate trajectories, which are a combination of translation, due to dielectrophoresis, and rotation, due to the Quincke effect. Our model provides a basis to study the collective dynamics of many particles in a general electric field.

We present a quantum theory for the dynamic structure factors in nonequilibrium, correlated, two-component systems such as plasmas or warm dense matter. The polarization function, which is needed as the input for the calculation of the structure factors, is calculated in nonequilibrium based on a perturbation expansion in the interaction strength. To make our theory applicable for x-ray scattering, a generalized Chihara decomposition for the total electron structure factor in nonequilibrium is derived. Examples are given and the influence of correlations and exchange on the structure and the x-ray-scattering spectrum are discussed for a model nonequilibrium distribution, as often encountered during laser heating of materials, as well as for two-temperature systems.

We derive an exact law for three-dimensional (3D) homogeneous compressible isothermal Hall magnetohydrodynamic turbulence, without the assumption of isotropy. The Hall current is shown to introduce new flux and source terms that act at the small scales (comparable or smaller than the ion skin depth) to significantly impact the turbulence dynamics. The law provides an accurate means to estimate the energy cascade rate over a broad range of scales covering the magnetohydrodynamic inertial range and the sub-ion dispersive range in 3D numerical simulations and in in situ spacecraft observations of compressible turbulence. This work is particularly relevant to astrophysical flows in which small-scale density fluctuations cannot be ignored such as the solar wind, planetary magnetospheres, and the interstellar medium.

For a one-dimensional (1D) momentum conserving system, intensive studies have shown that generally its heat current autocorrelation function (HCAF) tends to decay in a power-law manner and results in the breakdown of the Fourier heat conduction law in the thermodynamic limit. This has been recognized to be a dominant hydrodynamic effect. Here we show that, instead, the kinetic effect can be dominant in some cases and leads to the Fourier law for finite-size systems. Usually the HCAF undergoes a fast decaying kinetic stage followed by a long slowly decaying hydrodynamic tail. In a finite range of the system size, we find that whether the system follows the Fourier law depends on whether the kinetic stage dominates. Our Rapid Communication is illustrated by the 1D hard-core gas models with which the HCAF is derived analytically and verified numerically by molecular dynamics simulations.

Proper location of the mitotic spindle is critical for chromosome segregation and the selection of the cell division plane. However, how mitotic spindles sense cell size and shape to regulate their own position and orientation is still largely unclear. To investigate this question systematically, we used a general model by considering chromosomes, microtubule dynamics, and forces of various molecular motors. Our results show that in cells of various sizes and shapes, spindles can always be centered and oriented along the long axis robustly in the absence of other specified mechanisms. We found that the characteristic time of positioning and orientation processes increases with cell size. Spindles sense the cell size mainly by the cortical force in small cells and by the cytoplasmic force in large cells. In addition to the cell size, the cell shape mainly influences the orientation process. We found that more slender cells have a faster orientation process, and the final orientation is not necessarily along the longest axis but is determined by the radial profile and the symmetry of the cell shape. Finally, our model also reproduces the separation and repositioning of the spindle poles during the anaphase. Therefore, our work provides a general tool for studying the mitotic spindle across the whole mitotic phase.

We perform computational studies of static packings of a variety of nonspherical particles including circulo-lines, circulo-polygons, ellipses, asymmetric dimers, dumbbells, and others to determine which shapes form packings with fewer contacts than degrees of freedom (hypostatic packings) and which have equal numbers of contacts and degrees of freedom (isostatic packings), and to understand why hypostatic packings of nonspherical particles can be mechanically stable despite having fewer contacts than that predicted from naive constraint counting. To generate highly accurate force- and torque-balanced packings of circulo-lines and cir-polygons, we developed an interparticle potential that gives continuous forces and torques as a function of the particle coordinates. We show that the packing fraction and coordination number at jamming onset obey a masterlike form for all of the nonspherical particle packings we studied when plotted versus the particle asphericity $\mathcal{A}$, which is proportional to the ratio of the squared perimeter to the area of the particle. Further, the eigenvalue spectra of the dynamical matrix for packings of different particle shapes collapse when plotted at the same $\mathcal{A}$. For hypostatic packings of nonspherical particles, we verify that the number of “quartic” modes along which the potential energy increases as the fourth power of the perturbation amplitude matches the number of missing contacts relative to the isostatic value. We show that the fourth derivatives of the total potential energy in the directions of the quartic modes remain nonzero as the pressure of the packings is decreased to zero. In addition, we calculate the principal curvatures of the inequality constraints for each contact in circulo-line packings and identify specific types of contacts with inequality constraints that possess convex curvature. These contacts can constrain multiple degrees of freedom and allow hypostatic packings of nonspherical particles to be mechanically stable.

The lattice Boltzmann method has become one of the standard techniques for simulating a wide range of fluid flows. However, the intrinsic coupling of momentum and space discretization restricts the traditional lattice Boltzmann method to regular lattices. Alternative off-lattice Boltzmann schemes exist for both single- and multiphase flows that decouple the velocity discretization from the underlying spatial grid. The current study extends the applicability of these off-lattice methods by introducing a finite element formulation that enables simulating contact line dynamics for partially wetting fluids. This work exemplifies the implementation of the scheme and furthermore presents benchmark experiments that show the scheme reduces spurious currents at the liquid-vapor interface by at least two orders of magnitude compared to a nodal implementation and allows for predicting the equilibrium states accurately in the range of moderate contact angles.

According to Landauer's principle, the minimum heat emission required for computing is linked to logical entropy, or logical reversibility. The validity of Landauer's principle has been investigated for several decades and was finally demonstrated in recent experiments by showing that the minimum heat emission is associated with the reduction in logical entropy during a logically irreversible operation. Although the relationship between minimum heat emission and logical reversibility is being revealed, it is not clear how much free energy is required to be dissipated for a logically irreversible operation. In the present study, in order to reveal the connection between logical reversibility and free energy dissipation, we numerically demonstrated logically irreversible protocols using adiabatic superconductor logic. The calculation results of work during the protocol showed that, while the minimum heat emission conforms to Landauer's principle, the free energy dissipation can be arbitrarily reduced by performing the protocol quasistatically. The above results show that logical reversibility is not associated with thermodynamic reversibility, and that heat is not only emitted from logic devices but also absorbed by logic devices. We also formulated the heat emission from adiabatic superconductor logic during a logically irreversible operation at a finite operation speed.

The behavior of a stochastic resonate-and-fire neuron model based on a reduction of a fractional noise-driven generalized Langevin equation (GLE) with a power-law memory kernel is considered. The effect of temporally correlated random activity of synaptic inputs, which arise from other neurons forming local and distant networks, is modeled as an additive fractional Gaussian noise in the GLE. Using a first-passage-time formulation, in certain system parameter domains exact expressions for the output interspike interval (ISI) density and for the survival probability (the probability that a spike is not generated) are derived and their dependence on input parameters, especially on the memory exponent, is analyzed. In the case of external white noise, it is shown that at intermediate values of the memory exponent the survival probability is significantly enhanced in comparison with the cases of strong and weak memory, which causes a resonancelike suppression of the probability of spike generation as a function of the memory exponent. Moreover, an examination of the dependence of multimodality in the ISI distribution on input parameters shows that there exists a critical memory exponent ${\alpha }_c\approx 0.402$, which marks a dynamical transition in the behavior of the system. That phenomenon is illustrated by a phase diagram describing the emergence of three qualitatively different structures of the ISI distribution. Similarities and differences between the behavior of the model at internal and external noises are also discussed.