When living cells are observed at rest on a flat substrate, they can typically exhibit a rounded (symmetric) or an elongated (polarized) shape. Although the cells are apparently at rest, the active stress generated by the molecular motors continuously stretches and drifts the actin network, the cytoskeleton of the cell. In this paper we theoretically compare the energy stored and dissipated in this active system in two geometric configurations of interest: symmetric and polarized. We find that the stored energy is larger for a radially symmetric cell at low activation regime, while the polar configuration has larger strain energy when the active stress is beyond a critical threshold. Conversely, the dissipation of energy in a symmetric cell is always larger than that of a nonsymmetric one. By a combination of symmetry arguments and competition between surface and bulk stress, we argue that radial symmetry is an energetically expensive metastable state that provides access to an infinite number of lower-energy states, the polarized configurations.

The aim of the present paper is to provide a quantitative prediction of the stretch-induced crystallization in natural rubber, the exclusive reason for its history-dependent thermomechanical features. A constitutive model based on a micromechanism inspired molecular chain approach is formulated within the context of the thermodynamic framework. The molecular configuration of the partially crystallized single chain is analyzed and calculated by means of some statistical mechanical methods. The random thermal oscillation of the crystal orientation, considered as a continuous random variable, is treated by means of a representative angle. The physical expression of the chain free energy is derived according to a two-step strategy by separating crystallization and stretching. This strategy ensures that the stretch-induced part of the thermodynamic crystallization force is null at the initial instant and allows, without any additional constraint, the formulation of a simple linear relationship for the crystallinity evolution law. The model contains very few physically interpretable material constants to simulate the complex mechanism: two chain-scale constants, one crystallinity kinetics constant, three thermodynamic constants related to the newly formed crystallites, and a function controlling the crystal orientation with respect to the chain. The model is used to discuss some important aspects of the micromechanism and the macroresponse under the equilibrium state and the nonequilibrium state involved during stretching and recovery, and continuous relaxation.

We study a two-dimensional model for interacting colloidal particles which displays spontaneous clustering. Within this model we investigate the competition between the pinning to a periodic corrugation potential and a sideways constant pulling force which would promote a sliding state. For a few sample particle densities and amplitudes of the periodic corrugation potential we investigate the depinning from the statically pinned to the dynamically sliding regime. This sliding state exhibits the competition between a dynamics where entire clusters are pulled from a minimum to the next and a dynamics where single colloids or smaller groups leave a cluster and move across the corrugation energy barrier to join the next cluster downstream in the force direction. Both kinds of sliding states can occur either coherently across the entire sample or asynchronously: the two regimes result in different average mobilities. Finite temperature tends to destroy separate sliding regimes, generating a smoother dependence of the mobility on the driving force.

The transition to the chaos of the air flow between two vertical plates maintained at different temperatures is studied in the Boussinesq approximation. After the first bifurcation at critical Rayleigh number ${\mathrm{Ra}}_c$, the flow consists of two-dimensional (2D) corotating rolls. The stability of the 2D rolls is examined, confronting linear predictions with nonlinear integration. In all cases the 2D rolls are destabilized in the spanwise direction. Efficient linear stability analysis based on an Arnoldi method shows competition between two eigenmodes, corresponding to different spanwise wavelengths and different types of roll distortion. Nonlinear integration shows that the lower-wave-number mode is always dominant. A partial route to chaos is established through the nonlinear simulations. The flow becomes temporally chaotic for $\mathrm{Ra}=1.05{\mathrm{Ra}}_c$, but remains characterized by the spatial patterns identified by linear stability analysis. This highlights the complementary role of linear stability analysis and nonlinear simulation.

It is well recognized that there exist additional cubic terms of velocity in the lattice Boltzmann (LB) model based on the standard lattice. In this work, elimination of these cubic terms in the pseudopotential LB model for multiphase flow is investigated, where the force term and density gradient are considered. By retaining high-order ($\ge 3$) Hermite terms in the equilibrium distribution function and the discrete force term, as well as introducing correction terms in the LB equation, the additional cubic terms of velocity are entirely eliminated. With this technique, the computational simplicity of the pseudopotential LB model is well maintained. Numerical tests, including stationary and moving flat and circular interface problems, are carried out to show the effects of such cubic terms on the simulation of multiphase flow. It is found that the elimination of additional cubic terms is beneficial to reduce the numerical error, especially when the velocity is relatively large. Numerical results also suggest that these cubic terms mainly take effect in the interfacial region and that the density-gradient-related cubic terms are more important than the other cubic terms for multiphase flow.

We study explosive percolation (EP) on an Erdös-Rényi network for product rule (PR) and sum rule (SR). Initially, it was claimed that EP describes discontinuous phase transition; now it is well accepted as a probabilistic model for thermal continuous phase transition (CPT). However, no model for CPT is complete unless we know how to relate its observable quantities with those of thermal CPT. To this end, we define entropy and specific heat, redefine susceptibility, and show that they behave exactly like their thermal counterparts. We obtain the critical exponents $\nu ,\alpha ,\beta$, and $\gamma$ numerically and find that both PR and SR belong to the same universality class and obey Rushbrooke inequality.

A method is introduced to measure the entanglement entropy using a wavelet analysis. Using this method, the two-dimensional Haar wavelet transform of a configuration of Fortuin-Kasteleyn (FK) clusters is performed. The configuration represents a direct snapshot of spin-spin correlations since spin degrees of freedom are traced out in FK representation. A snapshot of FK clusters loses image information at each coarse-graining process by the wavelet transform. It is shown that the loss of image information measures the entanglement entropy in the Potts model.

We investigate in detail the process of formation of the multipeak low-temperature structure in the behavior of the specific heat capacity in frustrated magnetic systems in the framework of the exactly solvable antiferromagnetic spin-$1/2$ Ising model with the multisite interaction in the presence of the external magnetic field on the kagome-like Husimi lattice. The behavior of the entropy of the model is studied and exact values of the residual entropies of all ground states are found. It is shown that the multipeak structure in the behavior of the specific heat capacity is related to the formation of the multilevel hierarchical ordering in the system of all ground states of the model. Direct relation between the maximal number of peaks in the specific heat capacity behavior and the number of independent interactions in studied frustrated magnetic system is identified. The mechanism of the formation of the multipeak structure in the specific heat capacity is described and studied in detail, and it is generalized to frustrated magnetic systems with arbitrary numbers of independent interactions.

The simplest Klein-Gordon (KG) breather is a compacton on a string subject to the force of gravity in a frictionless V-shaped trough. Its dynamics, spectrum, and energy are discussed and it is compared to sine-Gordon breathers. A generalization of this problem consists of a charged string subject to the electrostatic force of two semi-infinite coplanar charged planes separated by a gap of constant width. For motion in the midplane between these planes, the string's displacement $u(x,t)$ satisfies the nonlinear KG equation $({\partial }_t^2-{\partial }_x^2)u=-{\text{tan}}^{-1}u$ in dimensionless form. Simulations of this equation reveal long-lived, breatherlike states, or “pseudobreathers,” which preserve shape and speed to high accuracy when Lorentz-transformed to simulate collisions.

In linear time-invariant systems acoustic reciprocity holds by the Onsager-Casimir principle of microscopic reversibility, and it can be broken only by odd external biases, nonlinearities, or time-dependent properties. Recently it was shown that one-dimensional lattices composed of a finite number of identical nonlinear cells with internal scale hierarchy and asymmetry exhibit nonreciprocity both locally and globally. Considering a single cell composed of a large scale nonlinearly coupled to a small scale, local dynamic nonreciprocity corresponds to vibration energy transfer from the large to the small scale, but absence of energy transfer (and localization) from the small to the large scale. This has been recently proven both theoretically and experimentally. Then, considering the entire lattice, global acoustic nonreciprocity has been recently proven theoretically, corresponding to preferential energy transfer within the lattice under transient excitation applied at one of its boundaries, and absence of similar energy transfer (and localization) when the excitation is applied at its other boundary. This work provides experimental validation of the global acoustic nonreciprocity with a one-dimensional asymmetric lattice composed of three cells, with each cell incorporating nonlinearly coupled large and small scales. Due to the intentional asymmetry of the lattice, low impulsive excitations applied to one of its boundaries result in wave transmission through the lattice, whereas when the same excitations are applied to the other end, they lead in energy localization at the boundary and absence of wave transmission. This global nonreciprocity depends critically on energy (i.e., the intensity of the applied impulses), and reduced-order models recover the nonreciprocal acoustics and clarify the nonlinear mechanism generating nonreciprocity in this system.

We investigate the response of quasiperiodically driven nonlinear systems exhibiting strange nonchaotic attractors (SNAs) to deterministic input signals. We show that if one uses two square waves in an aperiodic manner as input to a quasiperiodically driven double-well Duffing oscillator system, the response of the system can produce logical output controlled by such a forcing. Changing the threshold or biasing of the system changes the output to another logic operation and memory latch. The interplay of nonlinearity and quasiperiodic forcing yields logical behavior, and the emergent outcome of such a system is a logic gate. It is further shown that the logical behaviors persist even for an experimental noise floor. Thus the SNA turns out to be an efficient tool for computation.

Training an artificial neural network involves an optimization process over the landscape defined by the cost (loss) as a function of the network parameters. We explore these landscapes using optimization tools developed for potential energy landscapes in molecular science. The number of local minima and transition states (saddle points of index one), as well as the ratio of transition states to minima, grow rapidly with the number of nodes in the network. There is also a strong dependence on the regularization parameter, with the landscape becoming more convex (fewer minima) as the regularization term increases. We demonstrate that in our formulation, stationary points for networks with $N_h$ hidden nodes, including the minimal network required to fit the XOR data, are also stationary points for networks with $N_h+1$ hidden nodes when all the weights involving the additional node are zero. Hence, smaller networks trained on XOR data are embedded in the landscapes of larger networks. Our results clarify certain aspects of the classification and sensitivity (to perturbations in the input data) of minima and saddle points for this system, and may provide insight into dropout and network compression.

Cardiac alternans is a disturbance in heart rhythm that is linked to the onset of lethal cardiac arrhythmias. Mechanical perturbation control has been recently used to suppress alternans in cardiac tissue of relevant size. In this control strategy, cardiac tissue mechanics are perturbed via active tension generated by the heart's electrical activity, which alters the tissue's electric wave profile through mechanoelectric coupling. We analyze the effects of mechanical perturbation on the dynamics of a map model that couples the membrane voltage and active tension systems at the cellular level. Therefore, a two-dimensional iterative map of the heart beat-to-beat dynamics is introduced, and a stability analysis of the system of coupled maps is performed in the presence of a mechanical perturbation algorithm. To this end, a bidirectional coupling between the membrane voltage and active tension systems in a single cardiac cell is provided, and a discrete form of the proposed control algorithm, that can be incorporated in the coupled maps, is derived. In addition, a realistic electromechanical model of cardiac tissue is employed to explore the feasibility of suppressing alternans at cellular and tissue levels. Electrical activity is represented in two detailed ionic models, the Luo-Rudy 1 and the Fox models, while two active contractile tension models, namely a smooth variant of the Nash-Panfilov model and the Niederer-Hunter-Smith model, are used to represent mechanical activity in the heart. The Mooney-Rivlin passive elasticity model is employed to describe passive mechanical behavior of the myocardium.

We numerically examine ballistic active disks driven through a random obstacle array. Formation of a pinned or clogged state occurs at much lower obstacle densities for the active disks than for passive disks. As a function of obstacle density, we identify several distinct phases including a depinned fluctuating cluster state, a pinned single-cluster or jammed state, a pinned multicluster state, a pinned gel state, and a pinned disordered state. At lower active disk densities, a drifting uniform liquid forms in the absence of obstacles, but when even a small number of obstacles are introduced, the disks organize into a pinned phase-separated cluster state in which clusters nucleate around the obstacles, similar to a wetting phenomenon. We examine how the depinning threshold changes as a function of disk or obstacle density and find a crossover from a collectively pinned cluster state to a disordered plastic depinning transition as a function of increasing obstacle density. We compare this to the behavior of nonballistic active particles and show that as we vary the activity from completely passive to completely ballistic, a clogged phase-separated state appears in both the active and passive limits, while for intermediate activity, a readily flowing liquid state appears and there is an optimal activity level that maximizes the flux through the sample.

The present paper investigates the mean impact force exerted by a granular mass flowing down an incline and impacting a rigid wall of semi-infinite height. First, this granular flow-wall interaction problem is modeled by numerical simulations based on the discrete element method (DEM). These DEM simulations allow computing the depth-averaged quantities—thickness, velocity, and density—of the incoming flow and the resulting mean force on the rigid wall. Second, that problem is described by a simple analytic solution based on a depth-averaged approach for a traveling compressible shock wave, whose volume is assumed to shrink into a singular surface, and which coexists with a dead zone. It is shown that the dead-zone dynamics and the mean force on the wall computed from DEM can be reproduced reasonably well by the analytic solution proposed over a wide range of slope angle of the incline. These results are obtained by feeding the analytic solution with the thickness, the depth-averaged velocity, and the density averaged over a certain distance along the incline rather than flow quantities taken at a singular section before the jump, thus showing that the assumption of a shock wave volume shrinking into a singular surface is questionable. The finite length of the traveling wave upstream of the grains piling against the wall must be considered. The sensitivity of the model prediction to that sampling length remains complicated, however, which highlights the need of further investigation about the properties and the internal structure of the propagating granular wave.

We use a microscopically motivated generalized Langevin equation (GLE) approach to link the vibrational density of states (VDOS) to the dielectric response of orientational glasses (OGs). The dielectric function calculated based on the GLE is compared with experimental data for the paradigmatic case of two OGs: freon-112 and freon-113, around and just above $T_g$. The memory function is related to the integral of the VDOS times a spectral coupling function $\gamma \left({\omega }_p\right)$, which tells the degree of dynamical coupling between molecular degrees of freedom at different eigenfrequencies. The comparative analysis of the two freons reveals that the appearance of a secondary $\beta$ relaxation in freon-112 is due to cooperative dynamical coupling in the regime of mesoscopic motions caused by stronger anharmonicity (absent in freon-113) and is associated with the comparatively lower boson peak in the VDOS. The proposed framework brings together all the key aspects of glassy physics (VDOS with the boson peak, dynamical heterogeneity, dissipation, and anharmonicity) into a single model.

The Harada-Sasa equality elegantly connects the energy dissipation rate of a moving object with its measurable violation of the Fluctuation-Dissipation Theorem (FDT). Although proven for Langevin processes, its validity remains unclear for discrete Markov systems whose forward and backward transition rates respond asymmetrically to external perturbation. A typical example is a motor protein called kinesin. Here we show generally that the FDT violation persists surprisingly in the high-frequency limit due to the asymmetry, resulting in a divergent FDT violation integral and thus a complete breakdown of the Harada-Sasa equality. A renormalized FDT violation integral still well predicts the dissipation rate when each discrete transition produces a small entropy in the environment. Our study also suggests a way to infer this perturbation asymmetry based on the measurable high-frequency-limit FDT violation.

This article has to do with the derivation and solution of the Fokker-Planck equation associated to the momentum-integrated Wigner function of a particle subjected to a harmonic external field in contact with an ohmic thermal bath of quantum harmonic oscillators. The strategy employed is a simplified version of the phenomenological approach of Schramm, Jung, and Grabert of interpreting the operators as c numbers to derive the quantum master equation arising from a twofold transformation of the Wigner function of the entire phase space. The statistical properties of the random noise comes from the integral functional theory of Grabert, Schramm, and Ingold. By means of a single Wigner transformation, a simpler equation than that mentioned before is found. The Wigner function reproduces the known results of the classical limit. This allowed us to rewrite the underdamped classical Langevin equation as a first-order stochastic differential equation with time-dependent drift and diffusion terms.

In a crowd, individuals make different motion choices such as “moving to destination,” “following another pedestrian,” and “making a detour.” For the sake of convenience, the three direction choices are respectively called destination direction, following direction, and detour direction in this paper. Here, it is found that the featured direction choices could be inspired by the shape characteristics of the Voronoi diagram. To be specific, in the Voronoi cell of a pedestrian, the direction to a Voronoi node is regarded as a potential “detour” direction and the direction perpendicular to a Voronoi link is regarded as a potential “following” direction. A pedestrian generally owns several alternative Voronoi nodes and Voronoi links in a Voronoi cell, and the optimal detour and following direction are determined by considering related factors such as deviation. Plus the destination direction which is directly pointing to the destination, the three basic direction choices are defined in a Voronoi cell. In order to evaluate the Voronoi diagram based basic directions, the empirical trajectory data in both uni- and bi-directional flow experiments are extracted. A time series method considering the step frequency is used to reduce the original trajectories' swaying phenomena which might disturb the recognition of actual forward direction. The deviations between the empirical velocity direction and the basic directions are investigated, and each velocity direction is classified into a basic direction or regarded as an inexplicable direction according to the deviations. The analysis results show that each basic direction could be a potential direction choice for a pedestrian. The combination of the three basic directions could cover most empirical velocity direction choices in both uni- and bi-directional flow experiments.

A fundamental property of complex networks is the tendency for edges to cluster. The extent of the clustering is typically quantified by the clustering coefficient, which is the probability that a length-2 path is closed, i.e., induces a triangle in the network. However, higher-order cliques beyond triangles are crucial to understanding complex networks, and the clustering behavior with respect to such higher-order network structures is not well understood. Here we introduce higher-order clustering coefficients that measure the closure probability of higher-order network cliques and provide a more comprehensive view of how the edges of complex networks cluster. Our higher-order clustering coefficients are a natural generalization of the traditional clustering coefficient. We derive several properties about higher-order clustering coefficients and analyze them under common random graph models. Finally, we use higher-order clustering coefficients to gain new insights into the structure of real-world networks from several domains.