We consider the influence of active speed fluctuations on the dynamics of a $d$-dimensional active Brownian particle performing a persistent stochastic motion. The speed fluctuation brings about a dynamical anisotropy even in the absence of shape anisotropy. We use the Laplace transform of the Fokker-Planck equation to obtain exact expressions for time-dependent dynamical moments. Our results agree with direct numerical simulations and show several dynamical crossovers determined by the activity, persistence, and speed fluctuation. The dynamical anisotropy leads to a subdiffusive scaling in the parallel component of displacement fluctuation at intermediate times. The kurtosis remains positive at short times determined by the speed fluctuation, crossing over to a negative minimum at intermediate times governed by the persistence before vanishing asymptotically. The probability distribution of particle displacement obtained from numerical simulations in two dimensions shows two crossovers between compact and extended trajectories via two bimodal distributions at intervening times. While the speed fluctuation dominates the first crossover, the second crossover is controlled by persistence like in the wormlike chain model of semiflexible polymers.

We consider one-dimensional systems of all-to-all harmonically coupled particles with arbitrary masses, subject to two Langevin thermal baths. The couplings correspond to the mean-field limit of long-range interactions. Additionally, the particles can be subject to a harmonic on-site potential to break momentum conservation. Using the nonequilibrium Green's operator formalism, we calculate the transmittance, the heat flow, and local temperatures for arbitrary configurations of masses. For identical masses, we show analytically that the heat flux decays with the system size $N$ as $1/N$ regardless of the conservation or not of the momentum and of the introduction or not of a Kac factor. These results describe, in good agreement, the thermal behavior of systems with small heterogeneity in the masses.

Active chiral viscoelastic materials exhibit elastic responses perpendicular to the applied stresses, referred to as odd elasticity. We use a covariant formulation of viscoelasticity combined with an entropy production analysis to show that odd elasticity is not only present in active systems but also in broad classes of passive chiral viscoelastic fluids. In addition, we demonstrate that linear viscoelastic chiral solids require activity in order to manifest odd elastic responses. To model the phenomenon of passive odd viscoelasticity we propose a chiral extension of Jeffreys model. We apply our covariant formalism in order to derive the dispersion relations of hydrodynamic modes and obtain clear imprints of odd viscoelastic behavior.

Acoustic radiation force and torque arising from wave scattering are able to translate and rotate matter without contact. However, the existing research mainly focused on manipulating simple symmetrical geometries, neglecting the significance of geometric features. For the nonspherical geometries, the shape of the object strongly affects its scattering properties, and thus the radiation force and torque as well as the acoustophoretic process. Here, we develop an analytical framework to calculate the radiation force and torque exerted on the sound-hard or the sound-soft, axisymmetric particles excited by a user-customized transducer array based on a conformal transformation approach, capturing the significance of the geometric features. The derivation framework is established under the computation coordinate system (CCS), whereas the particle is assumed to be static. For the dynamic processes, the rotation of particle is converted as the opposite rotation of transducer array, achieved by employing a rotation transformation to tune the incident driving field in the CCS. Later, the obtained radiation force and torque in the CCS should be transformed back to the observation coordinate system for force and torque analysis. The radiation force and torque exerted on particles with different orientations are validated by comparing the full three-dimensional numerical solution in different phase distributions. It is found that the proposed method presents superior computational accuracy, high geometric adaptivity, and good robustness to various geometric features, while the computational efficiency is more than 100 times higher than that of the full numerical method. Furthermore, it is found that the dynamic trajectories of particles with different geometric features are completely different, indicating that the geometric features can be a potential degree of freedom to tune acoustophoretic process. The ability to predict the acoustophoretic process of nonspherical particles above a user-customized transducer array has improved our understanding of the effect of shape asymmetry, which can also be used to verify the effectiveness of acoustic tweezers in manipulating nonspherical objects.

Simulation-based inference enables learning the parameters of a model even when its likelihood cannot be computed in practice. One class of methods uses data simulated with different parameters to infer models of the likelihood-to-evidence ratio, or equivalently the posterior function. Here we frame the inference task as an estimation of an energy function parametrized with an artificial neural network. We present an intuitive approach, named MINIMALIST, in which the optimal model of the likelihood-to-evidence ratio is found by maximizing the likelihood of simulated data. Within this framework, the connection between the task of simulation-based inference and mutual information maximization is clear, and we show how several known methods of posterior estimation relate to alternative lower bounds to mutual information. These distinct objective functions aim at the same optimal energy form and therefore can be directly benchmarked. We compare their accuracy in the inference of model parameters, focusing on four dynamical systems that encompass common challenges in time series analysis: dynamics driven by multiplicative noise, nonlinear interactions, chaotic behavior, and high-dimensional parameter space.

Granular rheology is experimentally investigated in a vertical Couette-Poiseuille-like channel flow of photoelastic disks, where an erodible bed is sheared intermittently by an upward-moving shear band and a gravity-induced reverse flow. The shear band conforms to the existing nonlocal Eyring-like rheology but the bed exhibits discontinuous shear thickening from the Bagnold inertial regime near the band-bed interface to the Herschel-Bulkley plastic regime near the static wall. This newly discovered bed rheology is rate dependent and is associated with the fragility of the contact networks indicated by the statistics of local stress states inferred from the material photoelastic responses.

We present the solution of the weak noise theory (WNT) for the Kardar-Parisi-Zhang equation in one dimension at short time for flat initial condition (IC). The nonlinear hydrodynamic equations of the WNT are solved analytically through a connection to the Zakharov-Shabat (ZS) system using its classical integrability. This approach is based on a recently developed Fredholm determinant framework previously applied to the droplet IC. The flat IC provides the case for a nonvanishing boundary condition of the ZS system and yields a richer solitonic structure comprising the appearance of multiple branches of the Lambert function. As a byproduct, we obtain the explicit solution of the WNT for the Brownian IC, which undergoes a dynamical phase transition. We elucidate its mechanism by showing that the related spontaneous breaking of the spatial symmetry arises from the interplay between two solitons with different rapidities.

We study the question of universality in the two-dimensional spin-1 Baxter-Wu model in the presence of a crystal field $\mathrm{\Delta }$. We employ extensive numerical simulations of two types, providing us with complementary results: Wang-Landau sampling at fixed values of $\mathrm{\Delta }$ and a parallelized variant of the multicanonical approach performed at constant temperature $T$. A detailed finite-size scaling analysis in the regime of second-order phase transitions in the $(\mathrm{\Delta },T)$ phase diagram indicates that the transition belongs to the universality class of the four-state Potts model. Previous controversies with respect to the nature of the transition are discussed and attributed to the presence of strong finite-size effects, especially as one approaches the pentacritical point of the model.

We investigate the phase transitions of the $q$-state Brownian Potts model in two dimensions (2D) comprising Potts spins that diffuse like Brownian particles and interact ferromagnetically with other spins within a fixed distance. With extensive Monte Carlo simulations we find a continuous phase transition from a paramagnetic to a ferromagnetic phase even for $q>4$. This is in sharp contrast to the existence of a discontinuous phase transition in the equilibrium $q$-state Potts model in 2D with $q>4$. We present detailed numerical evidence for a continuous phase transition and argue that diffusion generated dynamical positional disorder suppresses phase coexistence leading to a continuous transition.

We investigate the unsteady heat (energy) transport in an infinite mass-in-mass chain with a given initial temperature profile. The chain consists of two sublattices: the $\beta$-Fermi-Pasta-Ulam-Tsingou (FPUT) chain and oscillators (of a different mass) connected to each FPUT particle. Initial conditions are such that initial kinetic temperatures of the FPUT particles and the oscillators are equal. Using the harmonic theory, we analytically describe evolution of these two temperatures in the ballistic regime. In particular, we derive a closed-form fundamental solution and solution for a sinusoidal initial temperature profile in the case when the oscillators are significantly lighter than the FPUT particles. The harmonic theory predicts that during the heat transfer the temperatures of sublattices are significantly different, while initially and finally (at large times) they are equal. This may look like an artifact of the harmonic approximation, but we show that it is not the case. Two distinct temperatures are also observed in the anharmonic case, even when the heat transport regime is no longer quasiballistic. We show that the value of the nonlinearity coefficient required to equalize the temperatures strongly depends on the particle mass ratio. If the oscillators are much lighter than the FPUT particles, then a fairly strong nonlinearity is required for the equalization.

We investigate restricted diffusion in a bounded domain towards a small partially reactive target in three- and higher-dimensional spaces. We propose a simple explicit approximation for the principal eigenvalue of the Laplace operator with mixed Robin-Neumann boundary conditions. This approximation involves the harmonic capacity and the surface area of the target, the volume of the confining domain, the diffusion coefficient, and the reactivity. The accuracy of the approximation is checked by using a finite-elements method. The proposed approximation determines also the mean first-reaction time, the long-time decay of the survival probability, and the overall reaction rate on that target. We identify the relevant lengthscale of the target, which determines its trapping capacity, and we investigate its relation to the target shape. In particular, we study the effect of target anisotropy on the principal eigenvalue by computing the harmonic capacity of prolate and oblate spheroids in various space dimensions. Some implications of these results in chemical physics and biophysics are briefly discussed.

The event graph representation of temporal networks suggests that the connectivity of temporal structures can be mapped to a directed percolation problem. However, similarly to percolation theory on static networks, this mapping is valid under the approximation that the structure and interaction dynamics of the temporal network are determined by its local properties, and, otherwise, it is maximally random. We challenge these conditions and demonstrate the robustness of this mapping in case of more complicated systems. We systematically analyze random and regular network topologies and heterogeneous link-activation processes driven by bursty renewal or self-exciting processes using numerical simulation and finite-size scaling methods. We find that the critical percolation exponents characterizing the temporal network are not sensitive to many structural and dynamical network heterogeneities, while they recover known scaling exponents characterizing directed percolation on low-dimensional lattices. While it is not possible to demonstrate the validity of this mapping for all temporal network models, our results establish the first batch of evidence supporting the robustness of the scaling relationships in the limited-time reachability of temporal networks.

Cell membranes are composed of a great variety of protein and lipid species with distinct unperturbed hydrophobic thicknesses. To achieve hydrophobic matching, the lipid bilayer tends to deform around membrane proteins so as to match the protein hydrophobic thickness at bilayer-protein interfaces. Such protein-induced distortions of the lipid bilayer hydrophobic thickness incur a substantial energy cost that depends critically on the bilayer-protein hydrophobic mismatch, while distinct conformational states of membrane proteins often show distinct hydrophobic thicknesses. As a result, hydrophobic interactions between membrane proteins and lipids can yield a rich interplay of lipid-protein organization and transitions in protein conformational state. We combine here the membrane elasticity theory of protein-induced lipid bilayer thickness deformations with the Landau-Ginzburg theory of lipid domain formation to systematically explore the coupling between local lipid organization, lipid and protein hydrophobic thickness, and protein-induced lipid bilayer thickness deformations in membranes with heterogeneous lipid composition. We allow for a purely mechanical coupling of lipid and protein composition through the energetics of protein-induced lipid bilayer thickness deformations as well as a chemical coupling driven by preferential interactions between particular lipid and protein species. We find that the resulting lipid-protein organization can endow membrane proteins with diverse and controlled mechanical environments that, via protein-induced lipid bilayer thickness deformations, can strongly influence protein function. The theoretical approach employed here provides a general framework for the quantitative prediction of how membrane thickness deformations influence the joint organization and function of lipids and proteins in cell membranes.

This study proposes a multiscale model combining molecular simulation and the lattice Boltzmann method (LBM) to explore gas flow behaviors with multiple transport mechanisms in nanoporous media of shale matrix. The gas adsorption characteristics in shale nanopores are first investigated by molecular simulations, which are then integrated and upscaled into the LBM model through a local adsorption density parameter. In order to adapt to high Knudsen number and nanoporous shale matrix, a multiple-relaxation-time pore-scale LBM model with a regularization procedure is developed. The combination of bounce-back and full diffusive boundary condition is adopted to take account of gas slippage and surface diffusion induced by gas adsorption. Molecular simulation results at the atomic scale show that gas adsorption behaviors are greatly affected by the pressure and pore size of the shale organic nanopore. At the pore scale, the gas transport behaviors with multiple transport mechanisms in nanoporous shale matrix are explored by the developed multiscale model. Simulation results indicate that pressure exhibits more significant influences on the transport behaviors of shale gas than temperature does. Compared with porosity, the average pore size of nanoporous shale matrix plays a more significant role in determining the apparent permeability of gas transport. The roles of the gas adsorption layer and surface diffusion in shale gas transport are discussed. It is observed that under low pressure, the gas adsorption layer has a positive influence on gas transport in shale matrix due to the strong surface diffusion effect. The nanoporous structure with the anisotropy characteristic parallel to the flow direction can enhance gas transport in shale matrix. The obtained results may provide underlying and comprehensive understanding of gas flow behaviors considering multiple transport mechanisms in shale matrix. Also, the proposed multiscale model can be considered as a powerful tool to invesigate the multiscale and multiphysical flow behaviors in porous media.

Many datasets give partial information about an ordering or ranking by indicating which team won a game, which item a user prefers, or who infected whom. We define a continuous spin system whose Gibbs distribution is the posterior distribution on permutations, given a probabilistic model of these interactions. Using the cavity method, we derive a belief propagation algorithm that computes the marginal distribution of each node's position. In addition, the Bethe free energy lets us approximate the number of linear extensions of a partial order and perform model selection between competing probabilistic models, such as the Bradley-Terry-Luce model of noisy comparisons and its cousins.

We consider a detailed-balance-violating dynamics whose stationary state is a prescribed Boltzmann distribution. Such dynamics have been shown to be faster than any equilibrium counterpart. We quantify the gain in convergence speed for a system whose energy landscape displays one and then an infinite number of energy barriers. In the latter case, we work with the mean-field disordered $p$ spin and show that the convergence to equilibrium or to the nonergodic phase is accelerated during both the $\beta$- and $\alpha$-relaxation stages. An interpretation in terms of trajectories in phase space and of an accidental fluctuation-dissipation theorem is provided.

We introduce a hybrid projection scheme that combines linear Mori projection and conditional Zwanzig projection techniques and use it to derive a generalized Langevin equation (GLE) for a general interacting many-body system. The resulting GLE includes (i) explicitly the potential of mean force (PMF) that describes the equilibrium distribution of the system in the chosen space of reaction coordinates, (ii) a random force term that explicitly depends on the initial state of the system, and (iii) a memory friction contribution that splits into two parts: a part that is linear in the past reaction-coordinate velocity and a part that is in general nonlinear in the past reaction coordinates but does not depend on velocities. Our hybrid scheme thus combines all desirable properties of the Zwanzig and Mori projection schemes. The nonlinear memory friction contribution is shown to be related to correlations between the reaction-coordinate velocity and the random force. We present a numerical method to compute all parameters of our GLE, in particular the nonlinear memory friction function and the random force distribution, from a trajectory in reaction coordinate space. We apply our method on the dihedral-angle dynamics of a butane molecule in water obtained from atomistic molecular dynamics simulations. For this example, we demonstrate that nonlinear memory friction is present and that the random force exhibits significant non-Gaussian corrections. We also present the derivation of the GLE for multidimensional reaction coordinates that are general functions of all positions in the phase-space of the underlying many-body system; this corresponds to a systematic coarse-graining procedure that preserves not only the correct equilibrium behavior but also the correct dynamics of the coarse-grained system.

We develop a general method to calculate the exact time dependence of the cumulants of the position of a tracer particle in a dense lattice gas of hardcore particles. More precisely, we calculate the cumulant-generating function associated with the position of a tagged particle at arbitrary time, and at leading order in the density of vacancies on the lattice. In particular, our approach gives access to the short-time dynamics of the cumulants of the tracer position, a regime in which few results are known. The generality of our approach is demonstrated by showing that it goes beyond the case of a symmetric 1D random walk and covers the important situations of (1) a biased tracer, (2) comblike structures, and (3) $d$-dimensional situations.

Loosely speaking, the Mpemba effect appears when hotter systems cool sooner or, in a more abstract way, when systems further from equilibrium relax faster. In this paper, we investigate the Mpemba effect in a molecular gas with nonlinear drag, both analytically (by employing the tools of kinetic theory) and numerically (direct simulation Monte Carlo of the kinetic equation and event-driven molecular dynamics). The analysis is carried out via two alternative routes, recently considered in the literature: first, the kinetic or thermal route, in which the Mpemba effect is characterized by the crossing of the evolution curves of the kinetic temperature (average kinetic energy), and, second, the stochastic thermodynamics or entropic route, in which the Mpemba effect is characterized by the crossing of the distance to equilibrium in probability space. In general, a nonmutual correspondence between the thermal and entropic Mpemba effects is found, i.e., there may appear the thermal effect without its entropic counterpart or vice versa. Furthermore, a nontrivial overshoot with respect to equilibrium of the thermal relaxation makes it necessary to revise the usual definition of the thermal Mpemba effect, which is shown to be better described in terms of the relaxation of the local equilibrium distribution. Our theoretical framework, which involves an extended Sonine approximation in which not only the excess kurtosis but also the sixth cumulant is retained, gives an excellent account of the behavior observed in simulations.

The motion of particles along channels of finite width is known to be hindered by either the presence of energy barriers along the channel direction or by variations in the width of the channel in the transverse direction (rugged channel). Remarkably, when both features are present, they can interact to produce a counterintuitive result: adding energy barriers to a rugged channel can enhance the rate of diffusion along it. This is the result of competing energetic and entropic effects. Under the approximation of particles instantaneously in equilibrium in the transverse direction, one can tailor the energy barriers to the ruggedness to recover free diffusion. However, such fine-tuning and potentially restrictive approximations are not necessary to observe an enhanced rate of diffusion as we demonstrate by adding a range of (non-fine-tuned) energy barriers to a channel of sinusoidally varying curvature. Furthermore, this was observed to hold for systems with a finite characteristic timescale for motion in the transverse direction, thus, suggesting that the phenomenon lends itself to be exploited for practical applications.