We show that the one-dimensional discrete nonlinear Schrödinger chain (DNLS) at finite temperature has three different dynamical regimes (ultralow-, low-, and high-temperature regimes). This has been established via (i) one-point macroscopic thermodynamic observables (temperature $T$, energy density $\epsilon$, and the relationship between them), (ii) emergence and disappearance of an additional almost conserved quantity (total phase difference), and (iii) classical out-of-time-ordered correlators and related quantities (butterfly speed and Lyapunov exponents). The crossover temperatures $T_{l-\mathrm{ul}}$ (between low- and ultra-low-temperature regimes) and $T_{h-l}$ (between the high- and low-temperature regimes) extracted from these three different approaches are consistent with each other. The analysis presented here is an important step forward toward the understanding of DNLS which is ubiquitous in many fields and has a nonseparable Hamiltonian form. Our work also shows that the different methods used here can serve as important tools to identify dynamical regimes in other interacting many-body systems.

We consider a population of two-dimensional oscillators with random couplings and explore the collective states. The coupling strength between oscillators is randomly quenched with two values, one of which is positive while the other is negative, and the oscillators can spatially move depending on the state variables for phase and position. We find that the system shows the phase transition from the incoherent state to the fully synchronized one at a proper ratio of the number of positive couplings to the total. The threshold is numerically measured and analytically predicted by the linear stability analysis of the fully synchronized state. It is found that the random couplings induce the long-term state patterns appearing for constant strength. The oscillators move to the places where the randomly quenched couplings work as if annealed. We further observe that the system with mixed randomnesses for quenched couplings shows the combination of the deformed patterns understandable with each annealed average.

This study evaluates data-driven models from a dynamical system perspective, such as unstable fixed points, periodic orbits, chaotic saddle, Lyapunov exponents, manifold structures, and statistical values. We find that these dynamical characteristics can be reconstructed much more precisely by a data-driven model than by computing directly from training data. With this idea, we predict the laminar lasting time distribution of a particular macroscopic variable of chaotic fluid flow, which cannot be calculated from a direct numerical simulation of the Navier-Stokes equation because of its high computational cost.

In this work we study shear reversals of dense non-Brownian suspensions composed of cohesionless elliptical particles. By numerical simulations, we show that a new fragility appears for frictionless ellipses in the flowing states, where particles can flow indefinitely in one direction at applied shear stresses but shear jam in the other direction upon shear stress reversal. This new fragility, absent in the isotropic particle case, is linked to the directional order of the elongated particles at steady shear and its reorientation at shear stress reversal, which forces the suspensions to pass through a more disordered state with an increased number of contacts in which it might get arrested.

The electrostatic screening of charge in one-dimensional confinement leads to long-range breakdown in electroneutrality within a nanopore. Through a series of continuum simulations, we demonstrate the principles of electroneutrality breakdown for electrolytes in one-dimensional confinement. We show how interacting pores in a membrane can counteract the phenomenon of electroneutrality breakdown, eventually returning to electroneutrality. Emphasis is placed on applying simplifying formulas to reduce the multidimensional partial differential equations into a single ordinary differential equation for the electrostatic potential. Dielectric mismatch between the electrolyte and membrane, pore aspect ratio, and confinement dimensionality are studied independently, outlining the relevance of electroneutrality breakdown in nanoporous membranes for selective ion transport and separations.

We study the realizability of simplicial complexes with a given pair of integer sequences, representing the node degree distribution and the facet size distribution, respectively. While the $s$-uniform variant of the problem is NP-complete when $s\ge 3$, we identify two populations of input sequences, most of which can be solved in polynomial time using a recursive algorithm that we contribute. Combining with a sampler for the simplicial configuration model [J.-G. Young et al., Phys. Rev. E 96, 032312 (2017)], we facilitate the efficient sampling of simplicial ensembles from arbitrary degree and size distributions. We find that, contrary to expectations based on dyadic networks, increasing the nodes' degrees reduces the number of loops in simplicial complexes. Our work unveils a fundamental constraint on the degree-size sequences and sheds light on further analyses of higher-order phenomena based on local structures.

The study of dynamical large deviations allows for a characterization of stationary states of lattice gas models out of equilibrium conditioned on averages of dynamical observables. The application of this framework to the two-dimensional random walk conditioned on partial currents reveals the existence of a dynamical phase transition between delocalized band dynamics and localized vortex dynamics. We present a numerical microscopic characterization of the phases involved and provide analytical insight based on the macroscopic fluctuation theory. A spectral analysis of the microscopic generator shows that the continuous phase transition is accompanied by spontaneous ${\mathbb{Z}}_2$-symmetry breaking whereby the stationary solution loses the reflection symmetry of the generator. Dynamical phase transitions similar to this one, which do not rely on exclusion effects or interactions, are likely to be observed in more complex nonequilibrium physics models.

In this work we study the thermal rectification efficiency, i.e., asymmetric heat flow, of a three-dimensional mass-graded anharmonic lattice of length $N$ and width $W$ by means of nonequilibrium molecular dynamics simulations. The obtained rectification, which is of the same order of magnitude as that of the corresponding one-dimensional lattice, saturates at low values of the aspect ratio $W/N$, consistent with the already known behavior of the corresponding heat fluxes of the homogeneous system under analogous conditions. The maximum rectification is obtained in the temperature range wherein no rectification could be obtained in other one-dimensional systems, as well as in the corresponding one-dimensional instance of the model studied herein.

We consider a spatially extended box-shaped wave field that consists of a plane wave (the condensate) in the middle and equals zero at the edges, in the framework of the focusing one-dimensional nonlinear Schrodinger equation. Within the inverse scattering transform theory, the scattering data for this wave field is presented by the continuous spectrum of the nonlinear radiation and the soliton eigenvalues together with their norming constants; the number of solitons $N$ is proportional to the box width. We remove the continuous spectrum from the scattering data and find analytically the specific corrections to the soliton norming constants that arise due to the removal procedure. The corrected soliton parameters correspond to symmetric in space $N$-soliton solution, as we demonstrate analytically in the paper. Generating this solution numerically for $N$ up to 1024, we observe that, at large $N$, it converges asymptotically to the condensate, representing its solitonic model. Our methods can be generalized for other strongly nonlinear wave fields, as we demonstrate for the hyperbolic secant potential, building its solitonic model as well.

We study the outflow of soft particles through quasi-two-dimensional hoppers with both experiments and simulations. The experiments utilize spheres made with hydrogel, silicone rubber, and glass. The hopper chamber has an adjustable exit width and tilt angle (the latter to control the magnitude of gravitational forcing). Our simulation mimics the experiments using purely two-dimensional soft particles with viscous interactions but no friction. Results from both simulations and experiments demonstrate that clogging is easier for reduced gravitational force or stiffer particles. For particles with low or no friction, the average number of particles in a clogging arch depends only on the ratio between hopper exit width and the mean particle diameter. In contrast, for the silicone rubber particles with larger frictional interactions, arches have more particles than the low friction cases. Additionally, an analysis of the number of particles left in the hopper when clogging occurs provides evidence for a hydrostatic pressure effect that is relevant for the clogging of soft particles, but less so for the harder (glass) or frictional (silicone rubber) particles.

Granular flows during a shear-induced mixing process are studied using discrete element methods. The aim is to understand the underlying elementary mechanisms of transition from unmixed to mixed phases for a granular material featuring a broad distribution of particles, which we investigate systematically by varying the strain rate and system size. Here the strain rate varies over four orders of magnitude and the system size varies from ten thousand to more than a million granules. A strain rate-dependent transition from quasistatic to purely inertial flow is observed. At the macroscopic scale, the contact stresses drop due to the formation of shear-induced instabilities that serves as an onset of granular flows and initiates mixing between the granules. The stress-drop displays a profound system size dependence. At the granular scale, mixing dynamics are correlated with the formation of shear bands, which result in significantly different timescales of mixing, especially for those regions that are close to the system walls and the bulk. Overall, our results reveal that although the transient dynamics display a generic behavior, these have a significant finite-size effect. In contrast, macroscopic behaviors at steady states have negligible system size dependence.

We experimentally study the late-time, highly nonlinear regime of the Rayleigh-Taylor instability in a decelerating phase. A series of laser-driven experiments is performed on the LULI2000 laser, in which the initial Atwood number is varied by adjusting the decelerating medium density. The high-power laser is used in a direct drive configuration to put into motion a solid target. Its rear side, which initially possesses a two-dimensional machined sinusoidal perturbations, expands and decelerates into a foam leading to a Rayleigh-Taylor unstable situation. The interface position and morphology are measured by time-resolved x-ray radiography. We develop a simple Atwood-dependent model describing the motion of the decelerating interface, from which its acceleration history is obtained. The measured amplitude of the instability, or mixing zone width, is then compared with late-time acceleration-dependent Rayleigh-Taylor instability models. The shortcomings of this classical model, when applied to high-energy-density conditions, are shown. This calls into question their uses for systems, where a shock wave is present, such as those found in laboratory astrophysics or in inertial confinement fusion.

We present a mechanism for efficiency increase in quantum heat engines containing internal energy levels that do not couple to the external work sink. The gain is achieved by using these levels to channel heat in a direction opposite to the one dictated by the second law. No quantum coherence, quantum correlations or ergotropy are required. A similar mechanism allows the engine to run “in reverse” and still produce useful work. We illustrate these ideas using a simple quantum Otto cycle in a coupled-spin system. We find this engine also exhibits other counterintuitive phenomenology. For example, its efficiency may increase as the temperature difference between the heat baths decreases. Conversely, it may cease to operate if the hotter bath becomes too hot or the colder bath too cold.

Grid cells in the dorsolateral band of the medial entorhinal cortex (dMEC) display strikingly regular periodic firing patterns on a lattice of positions in two-dimensional (2D) space. This helps animals to encode relative spatial location without reference to external cues. The dMEC is damaged in the early stages of Alzheimer's disease, which affects navigation ability of a disease victim, reducing the synaptic density of neurons in the network. Within an established two-dimensional continuous attractor neural network model of grid cell activity, we introduce neural sheet damage parametrized by radius and by the strength of the synaptic output for neurons in the damaged region. The mean proportionality of the grid field flow rate in the dMEC to the velocity of the model animal is maintained, but there is a broadened distribution of flow rates in the damaged case. This flow-rate-to-velocity proportionality is essential to establish coherent grid firing fields for individual grid cells for a roaming animal. When we examine the coherence of the grid cell firing field by studying Bragg peaks of the Fourier transformed lattice firing field intensity in both damaged and undamaged regions, we find that, for a wide range of damage radius and reduced synaptic strength, for undamaged model grid cells there is an incoherent firing field structure with only a single central peak. In the radius-damage plane this is adjacent to narrow bands of striped lattices (two additional Bragg peaks), which abut an orthorhombic pattern (four additional Bragg peaks), that abuts the undamaged hexagonal region (six additional Bragg peaks). Within the damaged region, grid cells show no Bragg peaks outside the central one, which shows reduced intensity with increasing damage, and outside the damaged region the central Bragg peak strength is largely unaffected. There is a reentrant region of normal grid firing fields for very large damage area. We anticipate that the modified grid cell behavior can be observed in noninvasive functional magnetic resonance imaging (fMRI) of the dMEC.

By introducing the notion of a dynamic overlap concentration scale, we identify additional universal features of the mechanical properties of active colloids. We codify these features by recognizing that the characteristic length scale of an active particle's trajectory, the run length, introduces a concentration scale ${\varphi }^{*}$. Large-scale simulations of repulsive active Brownian particles (ABPs) confirm that this run-length dependent concentration, the trajectory-space analog of the overlap concentration in polymer solutions, delineates distinct concentration regimes in which interparticle collisions alter particle trajectories. Using ${\varphi }^{*}$ and concentration scales associated with colloidal jamming, the mechanical equation of state for ABPs collapses onto a set of principal curves that contain several overlooked features. The inclusion of these features qualitatively alters previous predictions of the behavior for active colloids, as we demonstrate by computing the spinodal for a suspension of purely repulsive ABPs. Our findings suggest that dynamic overlap concentration scales should help unravel the behavior of active and driven systems.

We consider a mixture of active and passive run-and-tumble disks in an inhomogeneous environment where only half of the sample contains quenched disorder or pinning. The disks are initialized in a fully mixed state of uniform density. We identify several distinct dynamical phases as a function of motor force and pinning density. At high pinning densities and high motor forces, there is a two-step process initiated by a rapid accumulation of both active and passive disks in the pinned region, which produces a large density gradient in the system. This is followed by a slower species phase separation process where the inactive disks are shepherded by the active disks into the pin-free region, forming a nonclustered fluid and producing a more uniform density with species phase separation. For higher pinning densities and low motor forces, the dynamics becomes very slow and the system maintains a strong density gradient. For weaker pinning and large motor forces, a floating clustered state appears, and the time-averaged density of the system is uniform. We illustrate the appearance of these phases in a dynamic phase diagram.

We investigate the macroscopic dynamics of a two-fluid system with tetrahedral order. As all normal-fluid two-fluid systems one has—compared to a simple fluid—the velocity difference between the two subsystems and the concentration of one component as additional macroscopic variables. Depending on the type of system, the concentration can either be a conserved quantity or relax on a long, but finite timescale. Due to the existence of the tetrahedral order such a system breaks parity symmetry. Here we discuss physical systems without preferred direction in real space, meaning that our description applies to optically isotropic materials. We find a number of reversible as well as dissipative dynamic cross-coupling terms due to the additional octupolar order, when compared to a fluid mixture. As the most interesting cross-coupling term from an experimental point of view, we identify a dissipative cross-coupling between the relative velocity and the usual velocity gradients. Applying a shear flow in a plane, this dissipative cross-coupling leads to a velocity difference perpendicular to the shear plane. As a result one can obtain a spatially homogeneous oscillation of the relative velocity. In addition, this induced relative velocity can couple as a function of time and space to the concentration, which gives rise to an overdamped propagating soundlike mode, where the overdamping arises from the fact that velocity difference is a macroscopic variable and not strictly conserved. We also show that electric field gradients are connected with an analogous reversible cross-coupling and can lead in a planar shear geometry to an overdamped propagating mode as well.

Ultrathin membranes with nanoporous conduits show promise for ionic separations and desalination applications, but the mechanisms underlying the nonlinear ionic transport observed in these systems are not well understood. Here, we demonstrate how induced charge at membrane interfaces can lead to nonlinear ionic transport and voltage-dependent conductance through such channels. The application of an electric field on a polarizable membrane leads to induced charges at the membrane interfaces. The induced charges in turn are screened by diffuse charges in the electrolyte, which are acted upon by the electric field. For extremely thin membranes, the induced charge effect can be significant even for moderate applied voltages commonly used in experiments. We apply a continuum Poisson-Nernst-Planck model to characterize the current-voltage behavior of ultrathin membranes over a wide parameter space. The predictions of the model are compared to recent experiments on graphene and ${\mathrm{MoS}}_2$ membranes in an electric field. We expect the role of induced charge to be especially pronounced in the limit of atomically thin membranes.

We describe a dynamical state observed shortly above onset of the frozen wave instability. The transition to drifting waves, which are repeatedly created and destroyed, is a marked departure from the usual behavior of frozen waves, which are generally understood to remain motionless (on average) in the reference frame of the vibrating container. The spatial inhomogeneity of the underlying base flow, due both to the presence of the lateral walls and to the associated vibroequilibria effect, provides the driving mechanism. Energy arguments are used to understand the initial outward drift and the existence of a critical threshold which is estimated from the dependence of the drift velocity on the applied forcing. The dependence on container aspect ratio $\mathrm{\Gamma }$ is investigated, and drifting is seen to occur only when $1.5\lesssim \mathrm{\Gamma }\lesssim 3.5$.

We investigated the self-excited as well as optogalvanically stimulated heartbeat instability in RF discharge complex plasma. Three video cameras measured the motion of the microparticles, the plasma emission, and the laser-induced fluorescence simultaneously. Comprehensive studies of the optogalvanic control of the heartbeat instability revealed that the microparticle suspension can be stabilized by a continuous laser, whereas a modulated laser beam induces the void contraction either transiently or resonantly. The resonance occurred when the laser modulation frequency coincided with the frequency of small breathing oscillations of the microparticle suspension, which are known to be a prerequisite to the heartbeat instability. Based on the experimental results we suggest that the void contraction during the instability is caused by an abrupt void transition from the dim to the bright regime [Pikalev et al., Plasma Sources Sci. Technol. 30, 035014 (2021)]. In the bright regime, a time-averaged electric field at the void boundary heats the electrons causing bright plasma emission inside the void. The dim void has much lower electric field at the boundary and exhibits therefore no emission feature associated with it.