Stochastic phenomena are often described by Langevin equations, which serve as a mesoscopic model for microscopic dynamics. It has been known since the work of Parisi and Sourlas that reversible (or equilibrium) dynamics present supersymmetries (SUSYs). These are revealed when the path-integral action is written as a function not only of the physical fields, but also of Grassmann fields representing a Jacobian arising from the noise distribution. SUSYs leave the action invariant upon a transformation of the fields that mixes the physical and the Grassmann ones. We show that contrary to common belief, it is possible to extend the known reversible construction to the case of arbitrary irreversible dynamics, for overdamped Langevin equations with additive white noise—provided their steady state is known. The construction is based on the fact that the Grassmann representation of the functional determinant is not unique, and can be chosen so as to present a generalization of the Parisi-Sourlas SUSY. We show how such SUSYs are related to time-reversal symmetries and allow one to derive modified fluctuation-dissipation relations valid in nonequilibrium. We give as a concrete example the results for the Kardar-Parisi-Zhang equation.

Frustration effects caused by the presence of the six-site interaction in the ferromagnetic spin-$1/2$ Ising system on the kagome lattice are investigated in detail using the star kagomelike recursive lattice approximation. It is shown that although the model always exhibits the existence of only two standard phases (the ferromagnetic phase and the paramagnetic one) separated by the curve of the second-order phase transitions, depending on the value of the multisite interaction, the ferromagnetic phase splits into three different ground states in the zero-temperature limit with different magnetization and thermodynamic properties. The free energy of the model is derived, the residual entropies of all ground states are determined, and it is shown that the presence of the six-site multisite interaction leads to the formation of two highly macroscopically degenerated ground states in the studied ferromagnetic system, one of which is realized only for unique ratio of the six-site interaction to the ferromagnetic interaction. It is demonstrated that the existence of this highly macroscopically degenerated single-point-like ground state leads to appearance of the Schottky anomaly in the low-temperature behavior of the specific heat capacity in the vicinity of this ground state. It is also shown that the simultaneous presence of the frustration and of the second-order phase transitions in the studied model is responsible for the existence of even three local maxima in the temperature behavior of the specific heat capacity.

We investigate the modulational instability (MI) of a continuous wave (cw) under the combined effects of higher-order dispersions, self steepening and self-frequency shift, cubic, quintic, and septic nonlinearities. Using Maxwell's theory, an extended nonlinear Schrödinger equation is derived. The linear stability analysis of the cw solution is employed to extract an expression for the MI gain, and we point out its sensitivity to both higher-order dispersions and nonlinear terms. In particular, we insist on the balance between the sixth-order dispersion and nonlinearity, septic self-steepening, and the septic self-frequency shift terms. Additionally, the linear stability analysis of cw is confronted with the stability conditions for solitons. Different combinations of the dispersion parameters are proposed that support the stability of solitons and the occurrence of MI. This is confronted with full numerical simulations where the input cw gives rise to a broad range of behaviors, mainly related to nonlinear patterns formation. Interestingly, under the activation of MI, a suitable balance between the sixth-order dispersion and the septic self-frequency shift term is found to highly influence the propagation direction of the optical wave patterns.

Particle-type solutions are observed in out-of-equilibrium systems. These states can be motionless, oscillatory, or propagative depending on the injection and dissipation of energy. We investigate a family of localized standing waves based on a liquid-crystal light valve with spatiotemporal modulated optical feedback. These states are nonlinear waves in which energy concentrates in a localized and oscillatory manner. The organization of the family of solutions is characterized as a function of the applied voltage. Close to the reorientation transition, an amplitude equation allows us to elucidate the origin of these localized states and establish their bifurcation diagram. Theoretical findings are in qualitative agreement with experimental observations. Our results open the possibility of manipulating localized states induced by light, which can be used to expand and improve the storage and manipulation of information.

A model of opinion dynamics is introduced in which each individual's opinion is measured on a bounded continuous spectrum. Each opinion is influenced heterogeneously by every other opinion in the population. It is demonstrated that consensus, polarization and a spread of moderate opinions are all possible within this model. Using dynamic mean-field theory, we are able to identify the statistical features of the interactions between individuals that give rise to each of the aforementioned emergent phenomena. The nature of the transitions between each of the observed macroscopic states is also studied. It is demonstrated that heterogeneity of interactions between individuals can lead to polarization, that mostly antagonistic or contrarian interactions can promote consensus at a moderate opinion, and that mostly reinforcing interactions encourage the majority to take an extreme opinion.

Low-dimensional descriptions of spiking neural network dynamics are an effective tool for bridging different scales of organization of brain structure and function. Recent advances in deriving mean-field descriptions for networks of coupled oscillators have sparked the development of a new generation of neural mass models. Of notable interest are mean-field descriptions of all-to-all coupled quadratic integrate-and-fire (QIF) neurons, which have already seen numerous extensions and applications. These extensions include different forms of short-term adaptation considered to play an important role in generating and sustaining dynamic regimes of interest in the brain. It is an open question, however, whether the incorporation of presynaptic forms of synaptic plasticity driven by single neuron activity would still permit the derivation of mean-field equations using the same method. Here we discuss this problem using an established model of short-term synaptic plasticity at the single neuron level, for which we present two different approaches for the derivation of the mean-field equations. We compare these models with a recently proposed mean-field approximation that assumes stochastic spike timings. In general, the latter fails to accurately reproduce the macroscopic activity in networks of deterministic QIF neurons with distributed parameters. We show that the mean-field models we propose provide a more accurate description of the network dynamics, although they are mathematically more involved. Using bifurcation analysis, we find that QIF networks with presynaptic short-term plasticity can express regimes of periodic bursting activity as well as bistable regimes. Together, we provide novel insight into the macroscopic effects of short-term synaptic plasticity in spiking neural networks, as well as two different mean-field descriptions for future investigations of such networks.

We consider a dynamic network of individuals that may hold one of two different opinions in a two-party society. As a dynamical model, agents can endlessly create and delete links to satisfy a preferred degree, and the network is shaped by homophily, a form of social interaction. Characterized by the parameter $J\in [-1,1]$, the latter plays a role similar to Ising spins: agents create links to others of the same opinion with probability $(1+J)/2$ and delete them with probability $(1-J)/2$. Using Monte Carlo simulations and mean-field theory, we focus on the network structure in the steady state. We study the effects of $J$ on degree distributions and the fraction of cross-party links. While the extreme cases of homophily or heterophily ($J=\pm 1$) are easily understood to result in complete polarization or anti-polarization, intermediate values of $J$ lead to interesting features of the network. Our model exhibits the intriguing feature of an “overwhelming transition” occurring when communities of different sizes are subject to sufficient heterophily: agents of the minority group are oversubscribed and their average degree greatly exceeds that of the majority group. In addition, we introduce an original measure of polarization which displays distinct advantages over the commonly used average edge homogeneity.

We derive a mode-coupling theory (MCT) to describe the dynamics of a tracer particle that is embedded in a dense system of active Brownian particles (ABPs) in two spatial dimensions. The ABP undergo translational and rotational Brownian motion and are equipped with a fixed self-propulsion speed along their orientational vector that describes their active motility. The resulting equations of motion for the tagged-particle density-correlation functions describe the various cases of tracer dynamics close to the glass transition: that of a single active particle in a glass-forming passive host suspensions, that of a passive colloidal particle in a suspension of ABP, and that of active tracers in a bath of active particles. Numerical results are presented for these cases assuming hard-sphere interactions among the particles. The qualitative and quantitative accuracy of the theory is tested against event-driven Brownian dynamics (ED-BD) simulations of active and passive hard disks. Simulation and theory are found in quantitative agreement, provided one adjusts the overall density (as known from the passive description of glassy dynamics), and allows for a rescaling of self-propulsion velocities in the active host system. These adjustments account for the fact that ABP-MCT generally overestimates the tendency for kinetic arrest. We confirm in the simulations a peculiar feature of the transient and stationary dynamical density-correlation functions regarding their lack of symmetry under time reversal, demonstrating the nonequilibrium nature of the system and how it manifests itself in the theory.

From the detailed balance-like relation, we propose a measure, $K^{*}$, of a nonequilibrium steady-state (NESS) distance from equilibrium. We investigate in particular the NESS of a particle confined in a time-dependent harmonic potential of constant stiffness but with an ON-OFF state following a telegraph process. Experimental results coupled with simulations show that $K^{*}$ increases at slow switching rates (far from equilibrium) and approaches to zero at equilibrium conditions. Thus, the steady-state distribution together with $K^{*}$ fully characterizes a NESS.

A rotational sliceform (RS) forms a stiff, ringlike array of intersecting planar slices. Removing a few slices and disconnecting the ends of an RS enables the incomplete array to be collapsed scissorlike into a compact stack; it can be expanded smoothly as far as the original incomplete configuration, but not beyond. Its structured architecture, coupled to apparent mechanistic motion and a natural self-locking ability, expresses equivalently a novel deployable metamaterial, and we set out to determine its natural limits of motion for symmetrical and asymmetrical RS architectures. We first reconceptualize the RS as an array of plane-faced pyramidal cells bounded by rigid slices of zero thickness. The minimum articulation range from all cells is shown to set an upper bound on the range of motion of an incomplete RS, specifically, that symmetrical architectures can collapse fully while asymmetrical cannot and that expansion always stops at the design configuration. We also find that planar rotation of slices is not possible without distorting the original intersections. Each slice is then permitted to kink out-of-plane while preserving the initial geometry of each cell, in order to marshal compatible rotations of now compliant slices. Our analysis then reliably captures the deployment features: the minimum collapsed state, the degree of slice deformation as they rotate, and the limit of expansion.

Many of us have the experience of inflating balloons and twisting them into different shapes and animals. Snapping the balloon into two separate compartments is a necessary step that bears resemblance to the pinch-off phenomenon when a water droplet detaches from the faucet. In addition to testing whether balloons exhibit the properties of self-similarity and memory effect that are often associated with the latter event, we determine their phase diagram by experiments. It turns out that a common party balloon does not just snap, but can assume five more shapes, i.e., straight, necking, wrinkled, helix, and supercoil, depending on the twist angle and ratio of its length and diameter. Moreover, history also matters due to their prominent hysteresis. One may shift the phase boundary and/or reshuffle the phases by untwisting or lengthening the balloon at different twist angle and initial length. A heuristic minimal model is provided to obtain analytic expressions for the phase boundaries.

As porous media play an essential role in a variety of industrial applications, it is essential to understand their physical properties. Nowadays, the super-dimensional (SD) reconstruction algorithm is used to stochastically reconstruct a three-dimensional (3D) structure of porous media from a given two-dimensional image. This algorithm exhibits superiority in accuracy compared with classical algorithms because it learns information from the real 3D structure. However, owing to the short development time of the SD algorithm, it also has some limitations, such as inexact porosity characterization, long run time, blocking artifacts, and suboptimal accuracy that may be improved. To mitigate these limitations, this study presents the design of a special template that contains two parts of data (i.e., adjacent blocks and a central block); the proposed method matches adjacent blocks during reconstruction and assigns the matched central block to the area to be reconstructed. Furthermore, we design two important mechanisms during reconstruction: one for block matching and the other for porosity control. To verify the effectiveness of the proposed method compared with an existing SD method, both methods were tested on silica particle material and three homogeneous sandstones with different porosities; meanwhile, we compared the proposed method with a multipoint statistics method and a simulated annealing method. The reconstructed results were then compared with the target both visually and quantitatively. The experimental results indicate that the proposed method can overcome the aforementioned limitations and further improve the accuracy of existing methods. This method achieved 4–6 speedup factor compared with the traditional SD method.

An improved high-speed NMR profiling method is introduced that enables us to measure capillary action in thin, nontransparent porous media. Liquid profiles can be measured as fast as 10 ms with a spatial resolution of $14.5\mu \text{m}$. Capillary absorption of microliter-sized droplets into nontransparent nylon-6,6 porous membranes was studied. In nylon-6,6 sharp fronts obeying a one-dimensional Darcy model were observed. In paper the uptake did not follow this model: retardation in uptake, broader fluid fronts, and swelling was observed. These measurements nicely demonstrate the possibilities and versatility of the developed method.

The charging of an open quantum battery is investigated where the charger and the quantum battery interact with a common environment. At zero temperature, the stored energy of the battery is optimal as the charger and the quantum battery share the same coupling strength $(g_C=g_B)$. By contrast, in the presence of the quantum jump-based feedback control, the energy stored in the battery can be greatly enhanced for different couplings $(g_C>g_B)$. Considering the feasibility of the experiment, a model of Rydberg quantum battery is proposed with cascade-type atoms interacting with a dissipative optical cavity. The effective coupling strength between the charger (quantum battery) and the cavity field is hence adjustable and one can make the battery close to perfect excitation. The adverse factors of charging quantum batteries such as time delay for feedback, finite temperature, and spontaneous emission of Rydberg atoms are also discussed, and the result shows that the quantum battery is still able to retain a satisfactory energy storage effect.

We present exact results for the classical version of the out-of-time-order commutator (OTOC) for a family of power-law models consisting of $N$ particles in one dimension and confined by an external harmonic potential. These particles are interacting via power-law interaction of the form $\propto {\sum }_{\begin{array}{c}i,j=1(i\ne j)\end{array}}^N{|x_i-x_j|}^{-k}\forall k>1$ where $x_i$ is the position of the $i\text{th}$ particle. We present numerical results for the OTOC for finite $N$ at low temperatures and short enough times so that the system is well approximated by the linearized dynamics around the many-body ground state. In the large-$N$ limit, we compute the ground-state dispersion relation in the absence of external harmonic potential exactly and use it to arrive at analytical results for OTOC. We find excellent agreement between our analytical results and the numerics. We further obtain analytical results in the limit where only linear and leading nonlinear (in momentum) terms in the dispersion relation are included. The resulting OTOC is in agreement with numerics in the vicinity of the edge of the “light cone.” We find remarkably distinct features in OTOC below and above $k=3$ in terms of going from non-Airy behavior ($1<k<3$) to an Airy universality class ($k>3$). We present certain additional rich features for the case $k=2$ that stem from the underlying integrability of the Calogero-Moser model. We present a field theory approach that also assists in understanding certain aspects of OTOC such as the sound speed. Our findings are a step forward towards a more general understanding of the spatiotemporal spread of perturbations in long-range interacting systems.

We focus on the derivation of a general position-dependent effective diffusion coefficient to describe two-dimensional (2D) diffusion in a narrow and smoothly asymmetric channel of varying width under a transverse gravitational external field, a generalization of the symmetric channel case using the projection method introduced earlier by Kalinay and Percus [P. Kalinay and J. K. Percus, J. Chem. Phys. 122, 204701 (2005)]. To this end, we project the 2D Smoluchowski equation into an effective one-dimensional generalized Fick-Jacobs equation in the presence of constant force in the transverse direction. The expression for the diffusion coefficient given in Eq. (34) is our main result. This expression is a more general effective diffusion coefficient for narrow 2D channels in the presence of constant transverse force, which contains the well-known previous results for a symmetric channel obtained by Kalinay, as well as the limiting cases when the transverse gravitational external field goes to zero and infinity. Finally, we show that diffusivity can be described by the interpolation formula proposed by Kalinay, $D_0/{[1+(1/4)w^{\prime 2}\left(x\right)]}^{-\eta }$, where spatial confinement, asymmetry, and the presence of a constant transverse force can be encoded in $\eta$, which is a function of channel width ($w$), channel centerline, and transverse force. The interpolation formula also reduces to well-known previous results, namely, those obtained by Reguera and Rubi [D. Reguera and J. M. Rubi, Phys. Rev. E 64, 061106 (2001)] and by Kalinay [P. Kalinay, Phys. Rev. E 84, 011118 (2011)].

$A+B\to C$ reaction fronts describe a wide variety of natural and engineered dynamics, according to the specific nature of reactants and product. Recent works have shown that the properties of such reaction fronts depend on the system geometry, by focusing on one-dimensional plug flow radial injection. Here, we extend the theoretical formulation to radial deformation in two-dimensional systems. Specifically, we study the effect of a Poiseuille advective velocity profile on $A+B\to C$ fronts when A is injected radially into B at a constant flow rate in a confined axisymmetric system consisting of two parallel impermeable plates separated by a thin gap. We analyze the front dynamics by computing the temporal evolution of the average over the gap of the front position, the maximum production rate, and the front width. We further quantify the effects of the nonuniform flow on the total amount of product, as well as on its radial concentration profile. Through analytical and numerical analyses, we identify three distinct temporal regimes, namely (i) the early-time regime where the front dynamics is independent of the reaction, (ii) the transient regime where the front properties result from the interplay of reaction, diffusion that smooths the concentration gradients and advection, which stretches the spatial distribution of the chemicals, and (iii) the long-time regime where Taylor dispersion occurs and the system becomes equivalent to the one-dimensional plug flow case.

Undular bores, or dispersive shock waves, are nonstationary waves propagating as oscillatory transitions between two basic states, in which the oscillatory structure gradually expands and grows in amplitude with distance traveled. In this work we report an important mechanism of generation of nonlinear dispersive shock waves in solids. We demonstrate, using high-speed pointwise photoelasticity, the generation of undular bores in solid (polymethylmethacrylate) prestrained bars by natural and induced tensile fracture. For the distances relevant to our experiments, the viscoelastic extended Korteweg–de Vries equation is shown to provide very good agreement with the key observed experimental features for suitable choice of material parameters, while some local features at the front of the bore are also captured reasonably well by the linearization near the nonzero prestrain level. The experimental and theoretical approaches presented open avenues and analytical tools for the study and applications of dispersive shock waves in solids.

We study the low-temperature phase equilibria of a fluid confined in an open capillary slit formed by two parallel walls separated by a distance $L$ which are in contact with a reservoir of gas. The top wall of the capillary is of finite length $H$ while the bottom wall is considered of macroscopic extent. This system shows rich phase equilibria arising from the competition between two different types of capillary condensation, corner filling, and meniscus depinning transitions depending on the value of the aspect ratio $a=L/H$ and divides into three regimes: For long capillaries, with $a<2/\pi$, the condensation is of type I involving menisci which are pinned at the top edges at the ends of the capillary. For intermediate capillaries, with $2/\pi <a<1$, depending on the value of the contact angle the condensation may be of type I or of type II, in which the menisci overspill into the reservoir and there is no pinning. For short capillaries, with $a>1$, condensation is always of type II. In all regimes, capillary condensation is completely suppressed for sufficiently large contact angles which is determined explicitly. For long and intermediate capillaries, we show that there is an additional continuous phase transition in the condensed liquid-like phase, associated with the depinning of each meniscus as they round the upper open edges of the slit. Meniscus depinning is third-order for complete wetting and second-order for partial wetting. Detailed scaling theories are developed for these transitions and phase boundaries which connect with the theories of wedge (corner) filling and wetting encompassing interfacial fluctuation effects and the direct influence of intermolecular forces. We test several of our predictions using a fully microscopic density functional theory which allows us to study the two types of capillary condensation and its suppression at the molecular level for different aspect ratios and contact angles.

Pore structures and gas transport properties in porous separators for polymer electrolyte fuel cells are evaluated both experimentally and through simulations. In the experiments, the gas permeabilities of two porous samples, a conventional sample and one with low electrical resistivity, are measured by a capillary flow porometer, and the pore size distributions are evaluated with mercury porosimetry. Local pore structures are directly observed with micro-x-ray computed tomography (CT). In the simulations, the effective diffusion coefficients of oxygen and the air permeability in porous samples are calculated using random walk Monte Carlo simulations and computational fluid dynamics (CFD) simulations, respectively, based on the x-ray CT images. The calculated porosities and air permeabilities of the porous samples are in good agreement with the experimental values. The simulation results also show that the in-plane permeability is twice the through-plane permeability in the conventional sample, whereas it is slightly higher in the low-resistivity sample. The results of this study show that CFD simulation based on micro-x-ray CT images makes it possible to evaluate anisotropic gas permeabilities in anisotropic porous media.