Granular fluidity has been central to the development of nonlocal constitutive equations, which are necessary for characterizing nonlocal effects observed in experimental granular flow data. However, validation of these equations has been largely computational due to challenges in laboratory experiments. Specifically, the origin of the fluidity on a microscopic, single-particle level is still unproven. In this work, we present an experimental validation of a microscopic definition of granular fluidity, and show the importance of basal boundary conditions to the validity of the theory.

The response of a two-dimensional plasma crystal to an externally imposed initial perturbation has been explored using molecular dynamics (MD) simulations. A two-dimensional (2D) monolayer of micron-sized charged particles (dust) is formed in the plasma environment under certain conditions. The particles interacting via Yukawa pair potential are confined in the vertical ($\widehat{z}$) direction by an external parabolic confinement potential, which mimics the combined effect of gravity and the sheath electric field typically present in laboratory dusty plasma experiments. An external perturbation is introduced in the medium by displacing a small central region of particles in the vertical direction. The displaced particles start to oscillate in the vertical direction, and their dynamics get modulated through a parametric decay process generating beats. It has also been shown that the same motion is excited in the dynamics of unperturbed particles. A simple theoretical model is provided to understand the origin of the beat motions of particles. Additionally, in our simulations, concentric circular wavefronts propagating radially outward are observed on the surface of the monolayer. The physical mechanism and parametric dependence of the observed phenomena are discussed in detail. This research sheds light on the medium's ability to exhibit macroscopic softness, a pivotal characteristic of soft matter, while sustaining surface wave modes. Our findings are also relevant to other strongly coupled systems, such as colloids and classical one-component plasmas.

We report on the experimental observation of the focusing effect of a $50\mathrm{MeV}$ accelerator electron beam in a gas-discharge plasma target. The plasma is generated by igniting an electric discharge in two collinear quartz tubes, with the currents up to $1.5\mathrm{kA}$ flowing in opposite directions in either of the two tubes. In such plasma current configuration, the electron beam is defocused in the first discharge tube and focused with a stronger force in the second one. With symmetric plasma currents, asymmetric effects are, however, induced on the beam transport process and the beam radius is reduced by a factor of 2.6 compared to the case of plasma discharge off. Experimental results are supported by two-dimensional particle-in-cell simulations.

A lifting relation connecting the distribution function explicitly with the hydrodynamic variables is necessary for the Boltzmann equation-based mesoscopic approaches in order to correctly initialize a nonuniform hydrodynamic flow. We derive two lifting relations for Guo et al.'s total-energy double-distribution-function (DDF) kinetic model [Z. L. Guo et al., Phys. Rev. E 75, 036704 (2007)], one from the Hermite expansion of the conserved and nonconserved moments, and the second from the $O\left(\tau \right)$ Chapman-Enskog (CE) approximation of the Maxwellian exponential equilibrium. While both forms are consistent to the compressible Navier-Stokes-Fourier system theoretically, we stress that the latter may introduce numerical oscillations under the recently optimized discrete velocity models [Y. M. Qi et al., Phys. Fluids 34, 116101 (2022)], namely a 27 discrete velocity model of the seventh-order Gauss-Hermite quadrature (GHQ) accuracy (D3V27A7) for the velocity field combined with a 13 discrete velocity model of the fifth-order GHQ accuracy (D3V13A5) for the total energy. It is shown that the Hermite-expansion–based lifting relation can be alternatively derived from the latter approach using the truncated Hermite-polynomial equilibrium. Additionally, a relationship between the order of CE expansions and the truncated order of Hermite equilibria is developed to determine the minimal order of a Hermite equilibria required to recover any multiple-timescale macroscopic system. Next, three-dimensional compressible Taylor-Green vortex flows with different initial conditions and Ma numbers are simulated to demonstrate the effectiveness and potential issues of these lifting relations. The Hermite-expansion–based lifting relation works well in all cases, while the Chapman-Enskog-expansion–based lifting relation may produce numerical oscillations and a theoretical model is developed to predict such oscillations. Furthermore, the corresponding lifting relations for Qi et al.'s total energy DDF model [Y. M. Qi et al., Phys. Fluids 34, 116101 (2022)] are derived, and additional simulations are performed to illustrate the generality of our approach.

We investigate the force of flowing granular material on an obstacle. A sphere suspended in a discharging silo experiences both the weight of the overlaying layers and drag of the surrounding moving grains. In experiments with frictional hard glass beads, the force on the obstacle was practically flow-rate independent. In contrast, flow of nearly frictionless soft hydrogel spheres added drag to the gravitational force. The dependence of the total force on the obstacle diameter is qualitatively different for the two types of material: It grows quadratically with the obstacle diameter in the soft, low-friction material, while it grows much weaker, nearly linearly with the obstacle diameter, in the bed of glass spheres. In addition to the drag, the obstacle embedded in flowing low-friction soft particles experiences a total force from the top as if immersed in a hydrostatic pressure profile, but a much lower counterforce acting from below. In contrast, when embedded in frictional, hard particles, a strong pressure gradient forms near the upper obstacle surface.

A semi-Markov process method for obtaining general counting statistics for open quantum systems is extended to the scenario of resetting. The simultaneous presence of random resets and wave function collapses means that the quantum jump trajectories are no longer semi-Markov. However, focusing on trajectories and using simple probability formulas, general counting statistics can still be constructed from reset-free statistics. An exact tilted matrix equation is also obtained. The inputs of these methods are the survival distributions and waiting-time density distributions instead of quantum operators. In addition, a continuous-time cloning algorithm is introduced to simulate the large-deviation properties of open quantum systems. Several quantum optics systems are used to demonstrate these results.

The delayed Duffing equation, $x^{″}+\varepsilon x^{\prime }+x+x^3+cx(t-\tau )=0,$ admits a Hopf bifurcation which becomes singular in the limit $\varepsilon \to 0$ and $\tau =O\left(\varepsilon \right)\to 0$. To resolve this singularity, we develop an asymptotic theory where $x(t-\tau )$ is Taylor expanded in powers of $\tau .$ We derive a minimal system of ordinary differential equations that captures the Hopf bifurcation branch of the original delay differential equation. An unexpected result of our analysis is the necessity of expanding $x(t-\tau )$ up to third order rather than first order. Our work is motivated by laser stability problems exhibiting the same bifurcation problem as the delayed Duffing oscillator [Kovalev et al., Phys. Rev. E 103, 042206 (2021)]. Here we substantiate our theory based on the short delay limit by showing the overlap (matching) between our solution and two different asymptotic solutions derived for arbitrary fixed delays.

Using computer simulation and analytical theory, we study an active analog of the well-known Tonks gas, where active Brownian particles are confined to a periodic one-dimensional (1D) channel. By introducing the notion of a kinetic temperature, we derive an accurate analytical expression for the pressure and clarify the paradoxical behavior where active Brownian particles confined to 1D exhibit anomalous clustering but no motility-induced phase transition. More generally, this work provides a deeper understanding of pressure in active systems as we uncover a unique link between the kinetic temperature and swim pressure valid for active Brownian particles in higher dimensions.

While it is often assumed that liquid-vapor interfaces in nonequilibrium processes are in states of local thermodynamic equilibrium, this might not be the case for strong deviations from equilibrium. Clausius-Clapeyron equations for bulk properties yield a consistently defined temperature of the interface that is close to the liquid bulk temperature. The alternative interface temperature defined through the surface tension will be different for stronger nonequilibrium processes. Structural variables are introduced to extend the thermodynamic description of interfaces to a wider range of processes. Interfacial resistivities will depend on interface temperature as well as mass and heat flux through the interface.

Jamming is an emergent phenomenon wherein the local stability of individual particles percolates to form a globally rigid structure. However, the onset of rigidity does not imply that every particle becomes rigid, and indeed some remain locally unstable. These particles, if they become unmoored from their neighbors, are called rattlers, and their identification is critical to understanding the rigid backbone of a packing, as these particles cannot bear stress. The accurate identification of rattlers, however, can be a time-consuming process, and the currently accepted method lacks a simple geometric interpretation. In this manuscript, we propose two simpler classifications of rattlers in hard sphere systems based on the convex hull of contacting neighbors and the maximum inscribed sphere of the radical Voronoi cell, each of which provides geometric insight into the source of their instability. Furthermore, the convex hull formulation can be generalized to explore stability in hyperstatic soft sphere packings, spring networks, nonspherical packings, and mean-field non-central-force potentials.

To fill a gap in the literature about the specific dynamics of thermovibrational flow in a square cavity filled with a viscoelastic fluid when vibrations and the imposed temperature gradient are concurrent, a parametric investigation has been conducted to investigate the response of this system over a relatively wide subregion of the space of parameters (${\text{Pr}}_g=10$; viscosity ratio $\xi =0.5$; nondimensional frequency $\mathrm{\Omega }=25$, 50, 75, and 100; and ${\text{Ra}}_{\omega }\in [{\text{Ra}}_{\omega ,\text{cr}},3.3\times {10}^4]$, where ${\text{Ra}}_{\omega ,\text{cr}}$ is the critical vibrational Rayleigh number). Through solution of the governing nonlinear equations formulated in the framework of the finitely extensible nonlinear elastic Chilcott-Rallison paradigm, it is shown that the flow is prone to develop a unique hierarchy of bifurcations where initially subharmonic spatiotemporal regimes can be taken over by more complex states driven by the competition of disturbances with different symmetries if certain conditions are considered. What drives a wedge between the cases with parallel and perpendicular vibrations is essentially the existence of a threshold to be exceeded to produce convection in the former case. Nevertheless, these two configurations share some interesting properties, which are reminiscent of the resonances and antiresonances typical of multicomponent mechanical structures. Additional insights into these behaviors are gained through consideration of quantities representative of the kinetic and elastic energy globally possessed by the system and its sensitivity to the initial conditions.

A remarkable feature of two-dimensional turbulence is the transfer of energy from small to large scales. This process can result in the self-organization of the flow into large, coherent structures due to energy condensation at the largest scales. We investigate the formation of this condensate in a quasigeostropic flow in the limit of small Rossby deformation radius, namely the large-scale quasigeostrophic model. In this model potential energy is transferred up-scale while kinetic energy is transferred down-scale in a direct cascade. We focus on a jet mean flow and carry out a thorough investigation of the second-order statistics for this flow, combining a quasilinear analytical approach with direct numerical simulations. We show that the quasilinear approach applies in regions where jets are strong and is able to capture all second-order correlators in that region, including those related to the kinetic energy. This is a consequence of the blocking of the direct cascade by the mean flow in jet regions, suppressing fluctuation-fluctuation interactions. The suppression of the direct cascade is demonstrated using a local coarse-graining approach allowing us to measure space dependent interscale kinetic energy fluxes, which we show are concentrated in between jets in our simulations. We comment on the possibility of a similar direct cascade arrest in other two-dimensional flows, arguing that it is a special feature of flows in which the fluid element interactions are local in space.

Magnetorotational instability–driven (MRI-driven) turbulence and dynamo phenomena are analyzed using direct statistical simulations. Our approach begins by developing a unified mean-field model that combines the traditionally decoupled problems of the large-scale dynamo and angular momentum transport in accretion disks. The model consists of a hierarchical set of equations, capturing up to the second-order correlators, while a statistical closure approximation is employed to model the three-point correlators. We highlight the web of interactions that connect different components of stress tensors—Maxwell, Reynolds, and Faraday—through shear, rotation, correlators associated with mean fields, and nonlinear terms. We determine the dominant interactions crucial for the development and sustenance of MRI turbulence. Our general mean-field model for the MRI-driven system allows for a self-consistent construction of the electromotive force, inclusive of inhomogeneities and anisotropies. Within the realm of large-scale magnetic field dynamo, we identify two key mechanisms—the rotation-shear-current effect and the rotation-shear-vorticity effect—that are responsible for generating the radial and vertical magnetic fields, respectively. We provide the explicit (nonperturbative) form of the transport coefficients associated with each of these dynamo effects. Notably, both of these mechanisms rely on the intrinsic presence of large-scale vorticity dynamo within MRI turbulence.

We consider a dimer lattice of the Fermi-Pasta-Ulam-Tsingou (FPUT) type, where alternating linear couplings have a controllably small difference and the cubic nonlinearity ($\beta$-FPUT) is the same for all interaction pairs. We use a weakly nonlinear formal reduction within the lattice band gap to obtain a continuum, nonlinear Dirac-type system. We derive the Dirac soliton profiles and the model's conservation laws analytically. We then examine the cases of the semi-infinite and the finite domains and illustrate how the soliton solutions of the bulk problem can be glued to the boundaries for different types of boundary conditions. We thus explain the existence of various kinds of nonlinear edge states in the system, of which only one leads to the standard topological edge states observed in the linear limit. We finally examine the stability of bulk and edge states and verify them through direct numerical simulations, in which we observe a solitonlike wave setting into motion due to the instability.

We present a Hopfield-like autoassociative network for memories representing examples of concepts. Each memory is encoded by two activity patterns with complementary properties. The first is dense and correlated across examples within concepts, and the second is sparse and exhibits no correlation among examples. The network stores each memory as a linear combination of its encodings. During retrieval, the network recovers sparse or dense patterns with a high or low activity threshold, respectively. As more memories are stored, the dense representation at low threshold shifts from examples to concepts, which are learned from accumulating common example features. Meanwhile, the sparse representation at high threshold maintains distinctions between examples due to the high capacity of sparse, decorrelated patterns. Thus, a single network can retrieve memories at both example and concept scales and perform heteroassociation between them. We obtain our results by deriving macroscopic mean-field equations that yield capacity formulas for sparse examples, dense examples, and dense concepts. We also perform simulations that verify our theoretical results and explicitly demonstrate the capabilities of the network.

The pull-off adhesion force was measured by atomic force microscopy in sphere-plate geometry in water where a capillarylike behavior develops due to nanobubbles and was compared to the corresponding capillary adhesion in air. The sphere and the plate were coated with gold, and the pull-off adhesion force was measured as a function of the evolving surface roughness of the plate, and the retraction velocity of the interacting surfaces. In absolute magnitude, the pull-off force in air is larger than that in liquid by an order of magnitude or more, but in both cases, the pull-off force follows a monotonic decrease with increasing roughness. However, the relative decrement of the adhesion force in water was approximately 300%, and significantly higher than that in air for the same change of the rms roughness in the range $\sim 7\text{−}14$ nm. Finally, the adhesion force in water shows a relatively complex dependence on the retraction velocity of the interacting surfaces as the roughness increases due to possible deformation of the nanobubbles and the bridges they form between the surfaces.

It is well known that energy dissipation and finite size can deeply affect the dynamics of granular matter, often making usual hydrodynamic approaches problematic. Here we report on the experimental investigation of a small model system, made of ten beads constrained into a 1D geometry by a narrow vertical pipe and shaken at the base by a piston excited by a periodic wave. Recording the beads motion with a high frame rate camera allows to investigate in detail the microscopic dynamics and test hydrodynamic and kinetic models. Varying the energy, we explore different regimes from fully fluidized to the edge of condensation, observing good hydrodynamic behavior down to the edge of fluidization, despite the small system size. Density and temperature fields for different system energies can be collapsed by suitable space and time rescaling, and the expected constitutive equation holds very well when the particle diameter is considered. At the same time, the balance between dissipated and fed energy is not well described by commonly adopted dependence due to the up-down symmetry breaking. Our observations, supported by the measured particle velocity distributions, show a different phenomenological temperature dependence, which yields equation solutions in agreement with experimental results.

We investigate a symmetric logarithmic derivative (SLD) Fisher information for kinetic uncertainty relations (KURs) of open quantum systems described by the GKSL quantum master equation with and without the detailed balance condition. In a quantum kinetic uncertainty relation derived by Vu and Saito [Phys. Rev. Lett. 128, 140602 (2022)], the Fisher information of probability of quantum trajectory with a time-rescaling parameter plays an essential role. This Fisher information is upper bounded by the SLD Fisher information. For a finite time and arbitrary initial state, we derive a concise expression of the SLD Fisher information, which is a double time integral and can be calculated by solving coupled first-order differential equations. We also derive a simple lower bound of the Fisher information of quantum trajectory. We point out that the SLD Fisher information also appears in the speed limit based on the Mandelstam-Tamm relation by Hasegawa [Nat. Commun. 14, 2828 (2023)]. When the jump operators connect eigenstates of the system Hamiltonian, we show that the Bures angle in the interaction picture is upper bounded by the square root of the dynamical activity at short times, which contrasts with the classical counterpart.