The transition to turbulence in many shear flows proceeds along two competing routes, one linked with finite-amplitude disturbances and the other one originating from a linear instability, as in, e.g., boundary layer flows. The dynamical systems concept of an edge manifold has been suggested in the subcritical case to explain the partition of the state space of the system. This investigation is devoted to the evolution of the edge manifold when linear stability is added in such subcritical systems, a situation poorly studied despite its prevalence in realistic fluid flows. In particular, the fate of the edge state as a mediator of transition is unclear. A deterministic three-dimensional model is suggested, parametrized by the linear instability growth rate. The edge manifold evolves topologically, via a global saddle-loop bifurcation of the underlying invariant sets, from the separatrix between two attraction basins to the mediator between two transition routes. For larger instability rates, the stable manifold of the saddle point increases in codimension from 1 to 2 after an additional local pitchfork node bifurcation, causing the collapse of the edge manifold. As the growth rate is increased, three different regimes of this model are identified, each one associated with a flow case from the recent hydrodynamic literature. A simple nonautonomous generalization of the model is also suggested in order to capture the complexity of spatially developing flows.

The dynamic structure factor and other dynamic characteristics of strongly coupled one-component plasmas have been studied [Yu. V. Arkhipov et al., Phys. Rev. Lett. 119, 045001 (2017)] using the self-consistent version of the method of moments. Within any version of the latter, the system dielectric function satisfies all involved sum rules and other exact relations automatically, and the advantage of this version is that, in addition, the dynamic characteristics (the dynamic structure factor, the dispersion, and decay parameters of the collective modes) are all expressed in terms of the static ones (the static structure factor) without any adjustment to the simulation data. The approach outlined in the aforementioned Letter is justified in detail and applied mainly to the classical Coulomb systems achieving satisfactory agreement with new numerical simulation data. It is shown how the realm of applicability of the method can be extended to partly degenerate and multicomponent systems, even to simple liquids. Some additional theoretical results are presented in the Supplemental Material.

Using an analytical model and computer simulation, we show that the wakefield driven by an ultrashort laser pulse in high-density plasma periodically reverses its polarity due to the carrier-envelope phase shift of the driver. The wakefield polarity reversal occurs on spatial scales shorter than the typical length considered for electron acceleration with the laser-wakefield mechanism. Consequently, the energies of accelerated electrons are significantly affected. The results obtained are important for the laser-wakefield acceleration under the conditions relevant to present-day high-repetition-rate laser systems.

An autonomous out-of-equilibrium Maxwell's demon is used to reverse the natural direction of the heat flux between two electric circuits kept at different temperatures and coupled by the electric thermal noise. The demon does not process any information, but it achieves its goal by using a frequency-dependent coupling with the two reservoirs of the system. There is no mean energy flux between the demon and the system, but the total entropy production $(\mathrm{system}+\mathrm{demon})$ is positive. The demon can be power supplied by thermocouples. The system and the demon are ruled by equations similar to those of two coupled Brownian particles and of the Brownian gyrator. Thus our results pave the way to the application of autonomous out-of-equilibrium Maxwell's demons to coupled nanosystems at different temperatures.

Budó’s generalization [A. Budó, J. Chem. Phys. 17, 686 (1949)] of the Debye rotational diffusion model of dielectric relaxation of polar molecules to an assembly with internal interacting polar groups is extended to inertial anomalous diffusion. Thus, the theory can be applied both in the GHz and the THz regions, accounting for anomalous behavior as well as the necessary return to optical transparency at very high frequencies. The linking of both dispersion regions in a single model including anomalous effects is accomplished via a fractional Fokker-Planck equation in phase space based on the continuous time random walk ansatz. The latter is written via the Langevin equations for the stochastic dynamics of pairs of interacting heavy polar groups embedded in the frame of reference of a particular molecule or molecular dimer rotating about a space-fixed axis. The fractional Fokker-Planck equation is then converted to a three-term matrix differential recurrence equation for the statistical moments. This is solved in the frequency domain for the linear dielectric response using matrix continued fractions. Thus, one has the complex susceptibility $\chi \left(\omega \right)$ for extensive ranges of damping, group dipole moment ratio, and friction. The susceptibility, as inferred from the small oscillation limit, inherently comprises a low frequency (GHz) band with width depending on the anomalous parameter and a far-infrared (THz) or Poley peak of resonant character with a comblike structure of harmonic peaks. This behavior is due to the double transcendental nature of the after-effect function.

We study an Otto heat machine whose working substance is a single two-level system interacting with a cold thermal reservoir and with a squeezed hot thermal reservoir. By adjusting the squeezing or the adiabaticity parameter (the probability of transition) we show that our two-level system can function as a universal heat machine, either producing net work by consuming heat or consuming work that is used to cool or heat environments. Using our model we study the performance of these machine in the finite-time regime of the isentropic strokes, which is a regime that contributes to make them useful from a practical point of view.

We study a one-dimensional classical stochastic kinetically constrained model (KCM) inspired by Rydberg atoms in their “facilitated” regime, where sites can flip only if a single of their nearest neighbors is excited. We call this model “XOR-FA” to distinguish it from the standard Fredrickson-Andersen (FA) model. We describe the dynamics of the XOR-FA model, including its relation to simple exclusion processes in its domain wall representation. The interesting relaxation dynamics of the XOR-FA is related to the prominence of large dynamical fluctuations that lead to phase transitions between active and inactive dynamical phases as in other KCMs. By means of numerical tensor network methods we study in detail such transitions in the dynamical large deviation regime.

We compute rarefactive solitary wave solutions in a nonlinear lattice with nearest-neighbor interaction forces that are sublinear near the undeformed state. This setting includes bistable bonds governed by a double-well potential. In contrast to the prototypical Korteweg–de Vries–type delocalization, the obtained solutions feature a nontrivial sonic limit (Chapman-Jouguet regime) with nonzero energy and algebraic decay at infinity. In the bistable case the waves are strongly localized and have high energy over the entire velocity range. Direct numerical simulations suggest stability of the computed solitary waves. We consider several quasicontinuum models that mimic some features of the obtained solutions, including the nontrivial nature of the sonic limit, but fail to accurately approximate their core structure for all velocities in the bistable regime.

Bioelectrical patterns are established by spatiotemporal correlations of cell membrane potentials at the multicellular level, being crucial to development, regeneration, and tumorigenesis. We have conducted multicellular simulations on bioelectrical community effects and intercellular coupling in multicellular aggregates. The simulations aim at establishing under which conditions a local heterogeneity consisting of a small patch of cells can be stabilized against a large aggregate of surrounding identical cells which are in a different bioelectrical state. In this way, instructive bioelectrical information can be persistently encoded in spatiotemporal patterns of separated domains with different cell polarization states. The multicellular community effects obtained are regulated both at the single-cell and intercellular levels, and emerge from a delicate balance between the degrees of intercellular coupling in: (i) the small patch, (ii) the surrounding bulk, and (iii) the interface that separates these two regions. The model is experimentally motivated and consists of two generic voltage-gated ion channels that attempt to establish the depolarized and polarized cell states together with coupling conductances whose individual and intercellular different states permit a dynamic multicellular connectivity. The simulations suggest that community effects may allow the reprogramming of single-cell bioelectrical states, in agreement with recent experimental data. A better understanding of the resulting electrical regionalization can assist the electroceutical correction of abnormally depolarized regions initiated in the bulk of normal tissues as well as suggest new biophysical mechanisms for the establishment of target patterns in multicellular engineering.

We demonstrate that an insulating diblock copolymer system (DBC) filled with conductive fillers can be used as an electrically responsive soft composite that changes its conductivity in response to temperature-induced changes in its morphology. By combining the phase field model describing the morphology of the DBC system with the Monte Carlo simulations and the resistor network model describing electrical properties of the filler network, we calculate the conductivity of this composite. The calculated conductivity is found to essentially depend, in particular, on the temperature of the composite. Changing the temperature is shown to result in morphological changes in the DBC system causing the structural changes in the filler network. In particular, the order-disorder transition in the host DBC system is found to be accompanied by the conductor-insulator transition in the filler network. The effect of the difference between the affinities of the fillers for dissimilar copolymer blocks on the composite conductivity, as well as the effect of the repulsive and attractive interaction between fillers, is considered in detail.

To what extent the orientational order of self-propelling particles affects the dynamics of active bath-immersed ratchet devices remains unclear. We report experimental results of an inverse $\lambda$-like transition of the angular velocity of two ratchet gears in an active bath of self-propelling granular rods (SPRs) at different gear distances. The transition is caused by the phase transition in the orientational order of those SPRs located in the space between the gears. Brownian dynamics simulation of confined SPRs supports these observations. Thus, conditions for the upper bound efficiency of systems of active ratchet gears were obtained.

We introduce a two-dimensional active nematic with quenched disorder. We write the coarse-grained hydrodynamic equations of motion for slow variables, viz. density and orientation. Disorder strength is tuned from zero to large values. Results from the numerical solution of equations of motion as well as the calculation of two-point orientation correlation function using linear approximation shows that the ordered steady state follows a disorder dependent crossover from quasi-long-range order to short-range order. Such crossover is due to the pinning of $\pm 1/2$ topological defects in the presence of finite disorder, which breaks the system in uncorrelated domains. Finite disorder slows the dynamics of $+1/2$ defect, and it leads to slower growth dynamics. The two-point correlation functions for the density and orientation fields show good dynamic scaling but no static scaling for the different disorder strengths. Our findings can motivate experimentalists to verify the results and find applications in living and artificial apolar systems in the presence of a quenched disorder.

A collisional model of a confined quasi-two-dimensional granular mixture is considered to analyze homogeneous steady states. The model includes an effective mechanism to transfer the kinetic energy injected by vibration in the vertical direction to the horizontal degrees of freedom of grains. The set of Enskog kinetic equations for the velocity distribution functions of each component is derived first to analyze the homogeneous state. As in the one-component case, an exact scaling solution is found where the time dependence of the distribution functions occurs entirely through the granular temperature $T$. As expected, the kinetic partial temperatures $T_i$ of each component are different and, hence, energy equipartition is broken down. In the steady state, explicit expressions for the temperature $T$ and the ratio of partial kinetic temperatures $T_i/T_j$ are obtained by considering Maxwellian distributions defined at the partial temperatures $T_i$. The (scaled) granular temperature and the temperature ratios are given in terms of the coefficients of restitution, the solid volume fraction, the (scaled) parameters of the collisional model, and the ratios of mass, concentration, and diameters. In the case of a binary mixture, the theoretical predictions are exhaustively compared with both direct simulation Monte Carlo and molecular dynamics simulations with a good agreement. The deviations are identified to be originated in the non-Gaussianity of the velocity distributions and on microsegregation patterns, which induce spatial correlations not captured in the Enskog theory.

We investigate entanglement growth for a pair of coupled kicked rotors. For weak coupling, the growth of the entanglement entropy is found to be initially linear followed by a logarithmic growth. We calculate analytically the time after which the entanglement entropy changes its profile, and a good agreement with the numerical result is found. We further show that the different regimes of entanglement growth are associated with different rates of energy growth displayed by a rotor. At a large time, energy grows diffusively, which is preceded by an intermediate dynamical localization. The time span of intermediate dynamical localization decreases with increasing coupling strength. We argue that the observed diffusive energy growth is the result of one rotor acting as an environment to the other, which destroys the coherence. We show that the decay of the coherence is initially exponential followed by a power law.

In this work we study the energy transport in a one-dimensional system composed of two dissimilar Frenkel-Kontorova lattices connected by a time-modulated coupling and in contact with two heat reservoirs operating at different temperature by means of molecular dynamics simulations. There is a value of the driving frequency at which the heat flux takes its maximum value, a phenomenon termed thermal resonance. Structural modifications in the lattice strongly alter the way in which the external driving interacts with the phonon bands. The overlap of the latter in the harmonic regime of the model determines the frequency range wherein resonance emerges. Parameter dependencies by which the incoming heat flux can be directed to either of the heat reservoirs are examined as well. Our results may be conductive to further developments in designing thermal devices.

We study the position distribution of an active Brownian particle (ABP) in the presence of stochastic resetting in two spatial dimensions. We consider three different resetting protocols: (1) where both position and orientation of the particle are reset, (2) where only the position is reset, and (3) where only the orientation is reset with a certain rate $r$. We show that in the first two cases, the ABP reaches a stationary state. Using a renewal approach, we calculate exactly the stationary marginal position distributions in the limiting cases when the resetting rate $r$ is much larger or much smaller than the rotational diffusion constant $D_R$ of the ABP. We find that, in some cases, for a large resetting rate, the position distribution diverges near the resetting point; the nature of the divergence depends on the specific protocol. For the orientation resetting, there is no stationary state, but the motion changes from a ballistic one at short times to a diffusive one at late times. We characterize the short-time non-Gaussian marginal position distributions using a perturbative approach.

Networks of coupled nonlinear oscillators allow for the formation of nontrivial partially synchronized spatiotemporal patterns, such as chimera states, in which there are coexisting coherent (synchronized) and incoherent (desynchronized) domains. These complementary domains form spontaneously, and it is impossible to predict where the synchronized group will be positioned within the network. Therefore, possible ways to control the spatial position of the coherent and incoherent groups forming the chimera states are of high current interest. In this work we investigate how to control chimera patterns in multiplex networks of FitzHugh-Nagumo neurons, and in particular we want to prove that it is possible to remotely control chimera states exploiting the multiplex structure. We introduce a pacemaker oscillator within the network: this is an oscillator that does not receive input from the rest of the network but is sending out information to its neighbors. The pacemakers can be positioned in one or both layers. Their presence breaks the spatial symmetry of the layer in which they are introduced and allows us to control the position of the incoherent domain. We demonstrate how the remote control is possible for both uni- and bidirectional coupling between the layers. Furthermore we show which are the limitations of our control mechanisms when it is generalized from single-layer to multilayer networks.

Dynamical systems with long-delay feedback can exhibit complicated temporal phenomena, which once reorganized in a two-dimensional space are reminiscent of spatiotemporal behavior. In this framework, a normal forms description has been developed to reproduce the dynamics, and the opportunity to treat the corresponding variables as true space and time has since been established. However, recently, an alternative approach has been proposed [F. Marino and G. Giacomelli, Phys. Rev. E 98, 060201(R) (2018)] with a different interpretation of the variables involved, which better takes into account their physical character and allows for an easier determination of the normal forms. In this paper, we extend such idea and apply it to a number of paradigmatic examples, paving the way to a rethinking of the concept of spatiotemporal representation of long-delayed systems.

Highway freight transportation (HFT) plays an important role in the economic activities. Predicting HFT networks is not only scientifically significant in the understanding of the mechanism governing the formation and dynamics of these networks, but also of practical significance in highway planning and design for policymakers and truck allocation and route planning for logistic companies. In this work we apply parameter-free radiation models to predict the HFT network in mainland China and assess their predictive performance using metrics based on links and fluxes, which can be done in reference to the real directed and weighted HFT network between 338 Chinese cities constructed from about 15.06 million truck transportation records in five months. It is found that the radiation models exhibit relatively high accuracy in predicting links but low accuracy in predicting fluxes on links. Similar to gravity models, radiation models also suffer difficulty in predicting long-distance links and the fluxes on them. Nevertheless, the radiation models perform well in reproducing several scaling laws of the HFT network. The adoption of population or gross domestic product in the model has a minor impact on the results, and replacing the geographic distance by the path length taken by most truck drivers does not improve the results.