In linear time-invariant dynamical and acoustical systems, reciprocity holds by the Onsager-Casimir principle of microscopic reversibility, and this can be broken only by odd external biases, nonlinearities, or time-dependent properties. A concept is proposed in this work for breaking dynamic reciprocity based on irreversible nonlinear energy transfers from large to small scales in a system with nonlinear hierarchical internal structure, asymmetry, and intentional strong stiffness nonlinearity. The resulting nonreciprocal large-to-small scale energy transfers mimic analogous nonlinear energy transfer cascades that occur in nature (e.g., in turbulent flows), and are caused by the strong frequency-energy dependence of the essentially nonlinear small-scale components of the system considered. The theoretical part of this work is mainly based on action-angle transformations, followed by direct numerical simulations of the resulting system of nonlinear coupled oscillators. The experimental part considers a system with two scales—a linear large-scale oscillator coupled to a small scale by a nonlinear spring—and validates the theoretical findings demonstrating nonreciprocal large-to-small scale energy transfer. The proposed study promotes a paradigm for designing nonreciprocal acoustic materials harnessing strong nonlinearity, which in a future application will be implemented in designing lattices incorporating nonlinear hierarchical internal structures, asymmetry, and scale mixing.

In this paper we investigate the complex dynamics originated by a cross-diffusion-induced subharmonic destabilization of the fundamental subcritical Turing mode in a predator-prey reaction-diffusion system. The model we consider consists of a two-species Lotka-Volterra system with linear diffusion and a nonlinear cross-diffusion term in the predator equation. The taxis term in the search strategy of the predator is responsible for the onset of complex dynamics. In fact, our model does not exhibit any Hopf or wave instability, and on the basis of the linear analysis one should only expect stationary patterns; nevertheless, the presence of the nonlinear cross-diffusion term is able to induce a secondary instability: due to a subharmonic spatial resonance, the stationary primary branch bifurcates to an out-of-phase oscillating solution. Noticeably, the strong resonance between the harmonic and the subharmonic is able to generate the oscillating pattern albeit the subharmonic is below criticality. We show that, as the control parameter is varied, the oscillating solution ($\mathrm{sub}T$ mode) can undergo a sequence of secondary instabilities, generating a transition toward chaotic dynamics. Finally, we investigate the emergence of $\mathrm{sub}T$-mode solutions on two-dimensional domains: when the fundamental mode describes a square pattern, subharmonic resonance originates oscillating square patterns. In the case of subcritical Turing hexagon solutions, the internal interactions with a subharmonic mode are able to generate the so-called “twinkling-eyes” pattern.

The nature of particle-wall interactions is shown to have a profound impact on the well-known “Knudsen paradox” [or the “Knudsen minimum” effect, which refers to the decrease of the mass-flow rate of a gas with increasing Knudsen number Kn, reaching a minimum at $\text{Kn}\sim O\left(1\right)$ and increasing logarithmically with Kn as $\text{Kn}\to \infty ]$ in the acceleration-driven Poiseuille flow of rarefied gases. The nonmonotonic variation of the flow rate with Kn occurs even in a granular or dissipative gas in contact with thermal walls. The latter result is in contradiction with recent work [Alam et al., J. Fluid Mech. 782, 99 (2015)] that revealed the absence of the Knudsen minimum in granular Poiseuille flow for which the flow rate was found to decrease at large values of Kn. The above conundrum is resolved by distinguishing between “thermal” and “athermal” walls, and it is shown that, for both molecular and granular gases, the momentum transfer to athermal walls is much different than that to thermal walls which is directly responsible for the anomalous flow-rate variation with Kn in the rarefied regime. In the continuum limit of $\text{Kn}\to 0$, the athermal walls are shown to be closely related to “no-flux” (“adiabatic”) walls for which the Knudsen minimum does not exist either. A possible characterization of athermal walls in terms of (1) an effective specularity coefficient for the slip velocity and (2) a flux-type boundary condition for granular temperature is suggested based on simulation results.

We study a friction-controlled slide of a body excited by random motions of the foundation it is placed on. Specifically, we are interested in such quantities as displacement, traveled distance, and energy loss due to friction. We assume that the random excitation is switched off at some time (possibly infinite) and show that the problem can be treated in an analytic, explicit, manner. Particularly, we derive formulas for the moments of the displacement and distance, and also for the average energy loss. To accomplish that we use the Pugachev-Sveshnikov equation for the characteristic function of a continuous random process given by a system of SDEs. This equation is solved by reduction to a parametric Riemann boundary value problem of complex analysis.

The problem of random motion of harmonically trapped charged particles in a constant external magnetic field is studied. A generalized three-dimensional Langevin equation with a power-law memory kernel is used to model the interaction of Brownian particles with the complex structure of viscoelastic media (e.g., dusty plasmas). The influence of a fluctuating environment is modeled by an additive fractional Gaussian noise. In the long-time limit the exact expressions of the first-order and second-order moments of the fluctuating position for the Brownian particle subjected to an external periodic force in the plane perpendicular to the magnetic field have been calculated. Also, the particle's angular momentum is found. It is shown that an interplay of external periodic forcing, memory, and colored noise can generate a variety of cooperation effects, such as memory-induced sign reversals of the angular momentum, multiresonance versus Larmor frequency, and memory-induced particle confinement in the absence of an external trapping field. Particularly in the case without external trapping, if the memory exponent is lower than a critical value, we find a resonancelike behavior of the anisotropy in the particle position distribution versus the driving frequency, implying that it can be efficiently excited by an oscillating electric field. Similarities and differences between the behaviors of the models with internal and external noises are also discussed.

The application of noise to a nonlinear system can have the effect of increasing the signal transmission of the system through the phenomenon of stochastic resonance (SR). This paper presents an analytical characterization of the dependence of the signal transmission performance of an organic field-effect transistor (OFET) on external noise. Similarly to the threshold of a nonlinear system, the additive internal noise of the system can be used to control the emergence of SR. Internal noise or the addition of random numbers to the system enables one to observe the SR phenomenon in an OFET under an intrinsically nonresonant condition. Internal noise plays a thresholdlike role, but it functions in a different manner. The fluctuations in performance due to external noise become smaller when the effect of internal noise becomes dominant compared with that of the threshold. In conclusion, it is found that internal noise plays a robust thresholdlike role with respect to variations in external noise intensity.

This paper presents numerical studies of cavity polariton solitons (CPSs) in a resonantly pumped semiconductor microcavity with an imbedded spatial defect. In the bistable regime of the well-known homogeneous polariton condensate, with proper incident wave vector and pump strength, bright and/or dark cavity solitons can be found in the presence of a spatially confined potential. The minimum pump strength required to observe the CPSs or nonlinear localized states in this parametric pump scheme is therefore reported.

Polyhedral nematics are examples of exotic orientational phases that possess a complex internal symmetry, representing highly nontrivial ways of rotational symmetry breaking, and are subject to current experimental pursuits in colloidal and molecular systems. The classification of these phases has been known for a long time; however, their transitions to the disordered isotropic liquid phase remain largely unexplored, except for a few symmetries. In this work, we utilize a recently introduced non-Abelian gauge theory to explore the nature of the underlying nematic-isotropic transition for all three-dimensional polyhedral nematics. The gauge theory can readily be applied to nematic phases with an arbitrary point-group symmetry, including those where traditional Landau methods and the associated lattice models may become too involved to implement owing to a prohibitive order-parameter tensor of high rank or (the absence of) mirror symmetries. By means of exhaustive Monte Carlo simulations, we find that the nematic-isotropic transition is generically first-order for all polyhedral symmetries. Moreover, we show that this universal result is fully consistent with our expectation from a renormalization group approach, as well as with other lattice models for symmetries already studied in the literature. We argue that extreme fine tuning is required to promote those transitions to second-order ones. We also comment on the nature of phase transitions breaking the $\text{O}\left(3\right)$ symmetry in general cases.

We present a continuum model describing dissolution and growth of a crystal contact confined against a substrate. Diffusion and hydrodynamics in the liquid film separating the crystal and the substrate are modeled within the lubrication approximation. The model also accounts for the disjoining pressure and surface tension. Within this framework, we obtain evolution equations which govern the nonequilibrium dynamics of the crystal interface. Based on this model, we explore the problem of dissolution under an external load, known as pressure solution. We find that in steady state, diverging (power-law) crystal-surface repulsions lead to flat contacts with a monotonic increase of the dissolution rate as a function of the load. Forces induced by viscous dissipation then surpass those due to disjoining pressure at large enough loads. In contrast, finite repulsions (exponential) lead to sharp pointy contacts with a dissolution rate independent of the load and the liquid viscosity. Ultimately, in steady state, the crystal never touches the substrate when pressed against it. This result is independent from the nature of the crystal-surface interaction due to the combined effects of viscosity and surface tension.

We report patterns consisting of coexistence of synchronous and asynchronous states [for example, a granular gas co-existing with (i) bouncing bed, (ii) undulatory subharmonic waves, and (iii) Leidenfrost-like states] in experiments on vertically vibrated binary granular mixtures in a Hele-Shaw cell. Most experiments have been carried out with equimolar binary mixtures of glass and steel balls of same diameter by varying the total layer height ($F$) for a range of shaking acceleration ($\mathrm{\Gamma }$). All patterns as well as the related phase diagram in the ($\mathrm{\Gamma },F$) plane have been reproduced via molecular dynamics simulations of the same system. The segregation of heavier and lighter particles along the horizontal direction is shown to be the progenitor of such phase-coexisting patterns as confirmed in both experiment and simulation. At strong shaking we uncover a partial convection state in which a pair of convection rolls is found to coexist with a Leidenfrost-like state. The crucial role of the relative number density of two species on controlling the buoyancy-driven granular convection is demonstrated. The onset of horizontal segregation can be explained in terms of an anisotropic diffusion tensor.

We introduce a multispecies exclusion model where length-conserving probabilistic fusion and fission of the hard rods are allowed. Although all rods enter the system with the same initial length $\ell =1$, their length can keep changing, because of fusion and fission, as they move in a step-by-step manner towards the exit. Two neighboring hard rods of lengths ${\ell }_1$ and ${\ell }_2$ can fuse into a single rod of longer length $\ell ={\ell }_1+{\ell }_2$ provided $\ell \le N$. Similarly, length-conserving fission of a rod of length ${\ell }^{\prime }\le N$ results in two shorter daughter rods. Based on the extremum current hypothesis, we plot the phase diagram of the model under open boundary conditions utilizing the results derived for the same model under periodic boundary condition using mean-field approximation. The density profile and the flux profile of rods are in excellent agreement with computer simulations. Although the fusion and fission of the rods are motivated by similar phenomena observed in intraflagellar transport (IFT) in eukaryotic flagella, this exclusion model is too simple to account for the quantitative experimental data for any specific organism. Nevertheless, the concepts of “flux profile” and “transition zone” that emerge from the interplay of fusion and fission in this model are likely to have important implications for IFT and for other similar transport phenomena in long cell protrusions.

The field of disordered systems in statistical physics provides many simple models in which the competing influences of thermal and nonthermal disorder lead to new phases and nontrivial thermal behavior of order parameters. In this paper, we add a model to the subject by considering a disordered system where the state space consists of various orderings of a list. As in spin glasses, the disorder of such “permutation glasses” arises from a parameter in the Hamiltonian being drawn from a distribution of possible values, thus allowing nominally “incorrect orderings” to have lower energies than “correct orderings” in the space of permutations. We analyze a Gaussian, uniform, and symmetric Bernoulli distribution of energy costs, and, by employing Jensen's inequality, derive a simple condition requiring the permutation glass to always transition to the correctly ordered state at a temperature lower than that of the nondisordered system, provided that this correctly ordered state is accessible. We in turn find that in order for the correctly ordered state to be accessible, the probability that an incorrectly ordered component is energetically favored must be less than the inverse of the number of components in the system. We show that all of these results are consistent with a replica symmetric ansatz of the system. We conclude by arguing that there is no distinct permutation glass phase for the simplest model considered here and by discussing how to extend the analysis to more complex Hamiltonians capable of novel phase behavior and replica symmetry breaking. Finally, we outline an apparent correspondence between the presented system and a discrete-energy-level fermion gas. In all, the investigation introduces a class of exactly soluble models into statistical mechanics and provides a fertile ground to investigate statistical models of disorder.

Motivated by the qualitative picture of canonical typicality, we propose a refined formulation of the eigenstate thermalization hypothesis (ETH) for chaotic quantum systems. This formulation, which we refer to as subsystem ETH, is in terms of the reduced density matrix of subsystems. This strong form of ETH outlines the set of observables defined within the subsystem for which it guarantees eigenstate thermalization. We discuss the limits when the size of the subsystem is small or comparable to its complement. In the latter case we outline the way to calculate the leading volume-proportional contribution to the von Neumann and Renyi entanglment entropies. Finally, we provide numerical evidence for the proposal in the case of a one-dimensional Ising spin chain.

We derive the general relations between the maximum power, maximum efficiency, and minimum dissipation for the irreversible heat engine in a nonlinear response regime. In this context, we use the minimally nonlinear irreversible model and obtain the lower and upper bounds of the above relations for the asymmetric dissipation limits. These relations can be simplified further when the system possesses the time-reversal symmetry or antisymmetry. We find that our results are the generalization of various such relations obtained earlier for different heat engines.

We investigate equilibration and generalized thermalization of the quantum Harmonic chain under local quantum quench. The quench action we consider is connecting two disjoint harmonic chains of different sizes and the system jumps between two integrable settings. We verify the validity of the generalized Gibbs ensemble description for this infinite-dimensional Hilbert space system and also identify equilibration between the subsystems as in classical systems. Using Bogoliubov transformations, we show that the eigenstates of the system prior to the quench evolve toward the Gibbs Generalized Ensemble description. Eigenstates that are more delocalized (in the sense of inverse participation ratio) prior to the quench, tend to equilibrate more rapidly. Further, through the phase space properties of a generalized Gibbs ensemble and the strength of stimulated emission, we identify the necessary criterion on the initial states for such relaxation at late times and also find out the states that would potentially not be described by the generalized Gibbs ensemble description.

Nonlinear interaction between flexural modes is critical to heat conductivity and mechanical vibration of two-dimensional materials such as graphene. Much effort has been devoted to understand the underlying mechanism. In this paper, we examine solely the out-of-plane flexural modes and identify their energy flow pathway during thermalization process. The key is the development of a universal scheme that numerically characterizes the strength of nonlinear interactions between normal modes. In particular, for our square graphene system, the modes are grouped into four classes by their distinct symmetries. The couplings are significantly larger within a class than between classes. As a result, the equations for the normal modes in the same class as the initially excited one can be approximated by driven harmonic oscillators, therefore, they get energy almost instantaneously. Because of the hierarchical organization of the mode coupling, the energy distribution among the modes will arrive at a stable profile, where most of the energy is localized on a few modes, leading to the formation of “natural package” and metastable states. The dynamics for modes in other symmetry classes follows a Mathieu type of equation, thus, interclass energy flow, when the initial excitation energy is small, starts typically when there is a mode that lies in the unstable region in the parameter space of Mathieu equation. Due to strong coupling of the modes inside the class, the whole class will get energy and be lifted up by the unstable mode. This characterizes the energy flow pathway of the system. These results bring fundamental understandings to the Fermi-Pasta-Ulam problem in two-dimensional systems with complex potentials, and reveal clearly the physical picture of dynamical interactions between the flexural modes, which will be crucial to the understanding of their abnormal contribution to heat conduction and nonlinear mechanical vibrations.

Chimera states are complex spatio-temporal patterns which consist of coexisting domains of spatially coherent and incoherent dynamics in systems of coupled oscillators. In small networks, chimera states usually exhibit short lifetimes and erratic drifting of the spatial position of the incoherent domain. A tweezer feedback control scheme can stabilize and fix the position of chimera states. We analyze the action of the tweezer control in small nonlocally coupled networks of Van der Pol and FitzHugh-Nagumo oscillators, and determine the ranges of optimal control parameters. We demonstrate that the tweezer control scheme allows for stabilization of chimera states with different shapes, and can be used as an instrument for controlling the coherent domains size, as well as the maximum average frequency difference of the oscillators.

We have introduced a multiplex recurrence network approach by combining recurrence networks with the multiplex network approach in order to investigate multivariate time series. The potential use of this approach is demonstrated on coupled map lattices and a typical example from palaeobotany research. In both examples, topological changes in the multiplex recurrence networks allow for the detection of regime changes in their dynamics. The method goes beyond classical interpretation of pollen records by considering the vegetation as a whole and using the intrinsic similarity in the dynamics of the different regional vegetation elements. We find that the different vegetation types behave more similarly when one environmental factor acts as the dominant driving force.

Studying gene regulation networks in a single cell is an important, interesting, and hot research topic of molecular biology. Such process can be described by chemical master equations (CMEs). We propose a Hamilton-Jacobi equation method with finite-size corrections to solve such CMEs accurately at the intermediate region of switching, where switching rate is comparable to fast protein production rate. We applied this approach to a model of self-regulating proteins [H. Ge et al., Phys. Rev. Lett. 114, 078101 (2015)] and found that as a parameter related to inducer concentration increases the probability of protein production changes from unimodal to bimodal, then to unimodal, consistent with phenotype switching observed in a single cell.