Mechanically bonded fabrics account for a significant portion of nonwoven products, and serve many niche areas of nonwoven manufacturing. Such fabrics are characterized by layers of disordered fibrous webs, but we lack an understanding of how such microstructures determine bulk material response. Here we numerically determine the linear shear response of needle-punched fabrics modeled as cross-linked sheets of two-dimensional (2D) Mikado networks. We systematically vary the intra-sheet fiber density, inter-sheet separation distance, and direction of shear, and quantify the macroscopic shear modulus alongside the degree of affinity and energy partition. For shear parallel to the sheets, the response is dominated by intrasheet fibers and follows known trends for 2D Mikado networks. By contrast, shears perpendicular to the sheets induce a softer response dominated by either intrasheet or intersheet fibers depending on a quadratic relation between sheet separation and fiber density. These basic trends are reproduced and elucidated by a simple scaling argument that we provide. We discuss the implications of our findings in the context of real nonwoven fabrics.

A large-density-ratio and tunable-viscosity-ratio multicomponent multiphase pseudopotential lattice Boltzmann model is used to study the dissolution process of a bubble under pressure. The multi-relaxation-time collision operator, exact-difference-method external force scheme, and scaling coefficient $k$ are applied to ensure the numerical stability of the model. The influence of $k$ in the equation of state (EOS) and intermolecule interaction strength on the stationary bubble evolution process are discussed, and the effect of $k$ on thermodynamic consistency is also analyzed. The results indicate that adjusting the scaling coefficient in the EOS changes the surface tension and interface thickness, and that the gas-liquid interface width $w$ is proportional to $1/\sqrt{k}$. Considering the effect of $k$ on the surface tension, interface thickness, and thermodynamic consistency, the scaling coefficient should be between 0.6 and 1. Furthermore, the dissolution process of a single bubble under pressure is studied using the developed model, and it is found that the dissolution mass and concentration of dissolved gas increase linearly with increases in the pressure difference, and that the concentration of dissolved gas is proportional to the gas pressure after the fluid system reaches equilibrium. These results are consistent with Henry's law.

We report the experimental observation of Faraday waves on soft gels. These were obtained using agarose in a mechanically vibrated cylindrical container. Low driving frequencies induce subharmonic standing waves with spatial structure that conforms to the geometry of the container. We report the experimental observation of the first 15 resonant Faraday wave modes that can be defined by the mode number $(n,\ell )$ pair. We also characterize the shape of the instability tongue and show the complex dependence upon material properties can be understood as an elastocapillary effect.

Transfer entropy (TE) is widely used in time-series analysis to detect causal couplings between temporally evolving objects. As a coupling strength quantifier, the TE alone often seems insufficient, raising the question of its further interpretations. Here the TE is related to dynamical causal effects (DCEs) which quantify long-term responses of a coupling recipient to variations in a coupling source or in a coupling itself: Detailed relationships are established for a paradigmatic stochastic dynamical system of bidirectionally coupled linear overdamped oscillators, their practical applications and possible extensions are discussed. It is shown that two widely used versions of the TE (original and infinite-history) can become qualitatively distinct, diverging to different long-term DCEs.

In this paper, we examine the virial- and the potential-energy correlation for quasireal model systems. This correlation constitutes the framework of the theory of the isomorph in the liquid phase diagram commonly examined using simple liquids. Interestingly, our results show that for the systems characterized by structural anisotropy and flexible bonds, the instantaneous values of total virial and total potential energy are entirely uncorrelated. It is due to the presence of the intramolecular interactions because the contributions to the virial and potential energy resulting from the intermolecular interactions still exhibit strong linear dependence. Interestingly, in contrast to the results reported for simple liquids, the slope of the mentioned linear dependence is different than the values of the density scaling exponent. However, our findings show that for quasireal materials, the slope of dependence between the virial and potential energy (resulting from the intermolecular interactions) strongly depends on the interval of intermolecular distances that are taken into account. Consequently, the value of the slope of the discussed relationship, which enables satisfactory density scaling, can be obtained. Interestingly, this conclusion is supported by the results obtained for analogous systems without intermolecular attraction, for which the value the slope of the virial–potential-energy correlation is independent of considered intermolecular distances, directly corresponds to the exponent of the intermolecular repulsion, and finally leads to accurate density scaling.

We study an agent-based model of animals marking their territory and evading adversarial territory in one dimension with respect to the distribution of the size of the resulting territories. In particular, we use sophisticated sampling methods to determine it over a large part of territory sizes, including atypically small and large configurations, which occur with probability of less than ${10}^{-30}$. We find hints for the validity of a large deviation principle, the shape of the rate function for the right tail of the distribution, and insight into the structure of atypical realizations.

We present a method to renormalize stochastic differential equations subjected to multiplicative noise. The method is based on the widely used concept of effective potential in high-energy physics and has already been successfully applied to the renormalization of stochastic differential equations subjected to additive noise. We derive a general formula for the one-loop effective potential of a single ordinary stochastic differential equation (with arbitrary interaction terms) subjected to multiplicative Gaussian noise (provided the noise satisfies a certain normalization condition). To illustrate the usefulness (and limitations) of the method, we use the effective potential to renormalize a toy chemical model based on a simplified Gray-Scott reaction. In particular, we use it to compute the scale dependence of the toy model's parameters (in perturbation theory) when subjected to a Gaussian power-law noise with short time correlations.

Random walks process on networks plays a fundamental role in understanding the importance of nodes and the similarity of them, which has been widely applied in PageRank, information retrieval, and community detection, etc. An individual's memory has been proved to be crucial to affect network evolution and dynamical processes unfolding on the network. In this work, we study the random-walk process on an extended activity-driven network model by taking account of an individual's memory. We analyze how an individual's memory affects random-walk process unfolding on the network when the timescales of the processes of the random walk and the network evolution are comparable. Under the constraints of long-time evolution, we derive analytical solutions for the distribution of walkers at the stationary state and the mean first-passage time of the random-walk process. We find that, compared with the memoryless activity-driven model, an individual's memory enhances the activity fluctuation and leads to the formation of small clusters of mutual contacts with high activity nodes, which reduces a node's capability of gathering walkers, especially for the nodes with large activity, and memory also delays the mean first-passage time. The results on real networks also support the theoretical analysis and numerical results with artificial networks.

We consider the network constraints on the bounds of the assortativity coefficient, which aims to quantify the tendency of nodes with the same attribute values to be connected. The assortativity coefficient can be considered as the Pearson's correlation coefficient of node metadata values across network edges and lies in the interval $[-1,1]$. However, properties of the network, such as degree distribution and the distribution of node metadata values, place constraints upon the attainable values of the assortativity coefficient. This is important as a particular value of assortativity may say as much about the network topology as about how the metadata are distributed over the network—a fact often overlooked in literature where the interpretation tends to focus simply on the propensity of similar nodes to link to each other, without any regard on the constraints posed by the topology. In this paper we quantify the effect that the topology has on the assortativity coefficient in the case of binary node metadata. Specifically, we look at the effect that the degree distribution, or the full topology, and the proportion of each metadata value has on the extremal values of the assortativity coefficient. We provide the means for obtaining bounds on the extremal values of assortativity for different settings and demonstrate that under certain conditions the maximum and minimum values of assortativity are severely limited, which may present issues in interpretation when these bounds are not considered.

Quantifying the influence of microscopic details on the dynamics of development of the overall structure of a filamentous network is important in a number of biologically relevant contexts, but it is not obvious what order parameters can be used to adequately describe this complex process. In this paper we investigated the role of multivalent actin-binding proteins (ABPs) in reorganizing actin filaments into higher-order complex networks via a computer model of semiflexible filaments. We characterize the importance of local connectivity among actin filaments, as well as the global features of actomyosin networks. We first map the networks into local graph representations and then, using principles from network-theory order parameters, combine properties from these representations to gain insight into the heterogeneous morphologies of actomyosin networks at a global level. We find that ABPs with a valency greater than 2 promote filament bundles and large filament clusters to a much greater extent than bivalent multilinkers. We also show that active myosinlike motor proteins promote the formation of dendritic branches from a stalk of actin bundles. Our work motivates future studies to embrace network theory as a tool to characterize complex morphologies of actomyosin detected by experiments, leading to a quantitative understanding of the role of ABPs in manipulating the self-assembly of actin filaments into unique architectures that underlie the structural scaffold of a cell relating to its mobility and shape.

Tunicates are small invertebrates which possess a unique ability to reverse flow in their hearts. Scientists have debated various theories regarding how and why flow reversals occur. Here we explore the electrophysiological basis for reversals by simulating action potential propagation in an idealized model of the tubelike tunicate heart. Using asymptotic formulas for action potential duration and conduction velocity, we propose tunicate-specific parameters for a two-current ionic model of the action potential. Then, using a kinematic model, we derive analytical criteria for reversals to occur. These criteria inform subsequent numerical simulations of action potential propagation in a fiber paced at both ends. In particular, we explore the role that variability of pacemaker firing rates plays in generating reversals, and we identify various favorable conditions for triggering retrograde propagation. Our analytical framework extends to other species; for instance, it can be used to model competition between the sinoatrial node and abnormal ectopic foci in human heart tissue.

Transient or sustained permeability transition pore (PTP) opening is important in normal physiology or cell death, respectively. These are closely linked to ${\mathrm{Ca}}^{2+}$ and reactive oxygen species (ROS). The entry of ${\mathrm{Ca}}^{2+}$ into mitochondria regulates ROS production, and both ${\mathrm{Ca}}^{2+}$ and ROS trigger PTP opening. In addition to this feedforward loop, there exist four feedback loops in the ${\mathrm{Ca}}^{2+}$-ROS-PTP system. ROS promotes ${\mathrm{Ca}}^{2+}$ entering (F1) and induces further ROS generation (F2), forming two positive feedback loops. PTP opening results in the efflux of ${\mathrm{Ca}}^{2+}$ (F3) and ROS (F4) from the mitochondria, forming two negative feedback loops. Owing to these complexities, we construct a mathematical model to dissect the roles of these feedback loops in the dynamics of PTP opening. The qualitative agreement between simulation results and recent experimental observations supports our hypothesis that under physiological conditions the PTP opens in an oscillatory state, while under pathological conditions it opens in a high steady state. We clarify that the negative feedback loops are responsible for producing oscillations, wherein F3 plays a more prominent role than F4; whereas the positive feedback loops are beneficial for maintaining oscillation robustness, wherein F1 has a more dominant role than F2. Furthermore, we manifest that the proper increase in negative feedback strength or decrease in positive feedback strength not only facilitates the occurrence of oscillations and thus protects the system against a high steady state, but also assists in lowering the oscillation peak. This study may provide potential therapeutic strategies in treating neurodegenerative diseases due to PTP dysfunction.

Fast shocks that form in optically thick media are mediated by Compton scattering and, if relativistic, pair creation. Since the radiation force acts primarily on electrons and positrons, the question arises of how the force is mediated to the ions which are the dominant carriers of the shock energy. It has been widely thought that a small charge separation induced by the radiation force generates an electric field inside the shock that decelerates the ions. In this paper we argue that, while this is true in subrelativistic shocks which are devoid of positrons, in relativistic radiation mediated shocks (RRMS), which are dominated by newly created $e^{+}{e}^{-}$ pairs, additional coupling is needed, owing to the opposite electric force acting on electrons and positrons. Specifically, we show that dissipation of the ions energy must involve collective plasma interactions. By constructing a multifluid model for RRMS that incorporates friction forces, we estimate that momentum transfer between electrons and positrons (and/or ions) via collective interactions on scales of tens to thousands of proton skin depths, depending on whether friction is effective only between $e^{+}{e}^{-}$ pairs or also between pairs and ions, is sufficient to couple all particles and radiation inside the shock into a single fluid. This leaves open the question of whether in relativistic RMS particles can effectively accelerate to high energies by scattering off plasma turbulence. Such acceleration might have important consequences for relativistic shock breakout signals.

A laser plasma wake field in a single-color femtosecond laser filament determines the acceleration of ionized electrons, which affects the intensity and bandwidth of the emitted terahertz wave and is important for understanding the fundamental nonlinear process of THz generation. Since the THz wave generated by a laser wake field is extremely small and easily hidden by other THz generation mechanisms, no method exists to measure this wake field directly. In this paper, a simple and stable method for determining the amplitude of the laser plasma wake field is presented. Based on the cancellation of a positive laser plasma wake field and an external negative electric field, the “zero point” of the intensity of the generated THz wave at some frequency can be used to determine the exact amplitude of the corresponding laser plasma wake field. This finding opens an avenue toward the clarification of ultrafast electronic dynamic processes in laser-induced plasmas.

Using high-precision Monte Carlo simulations based on a parallel version of the Wang-Landau algorithm and finite-size scaling techniques, we study the effect of quenched disorder in the crystal-field coupling of the Blume-Capel model on a square lattice. We mainly focus on the part of the phase diagram where the pure model undergoes a continuous transition, known to fall into the universality class of a pure Ising ferromagnet. A dedicated scaling analysis reveals concrete evidence in favor of the strong universality hypothesis with the presence of additional logarithmic corrections in the scaling of the specific heat. Our results are in agreement with an early real-space renormalization-group study of the model as well as a very recent numerical work where quenched randomness was introduced in the energy exchange coupling. Finally, by properly fine tuning the control parameters of the randomness distribution we also qualitatively investigate the part of the phase diagram where the pure model undergoes a first-order phase transition. For this region, preliminary evidence indicate a smoothing of the transition to second-order with the presence of strong scaling corrections.

Trusting in others and reciprocating that trust with trustworthy actions are crucial to successful and prosperous societies. The trust game has been widely used to quantitatively study trust and trustworthiness, involving a sequential exchange between an investor and a trustee. Deterministic evolutionary game theory predicts no trust and no trustworthiness, whereas the behavioral experiments with the one-shot anonymous trust game show that people substantially trust and respond trustworthily. To explain these discrepancies, previous works often turn to additional mechanisms, which are borrowed from other games such as the prisoner's dilemma. Although these mechanisms lead to the evolution of trust and trustworthiness to an extent, the optimal or the most common strategy often involves no trustworthiness. In this paper, we study the impact of asymmetric demographic parameters (e.g., different population sizes) on game dynamics of the trust game. We show that, in a weak-mutation limit, stochastic evolutionary dynamics with the asymmetric parameters can lead to the evolution of high trust and high trustworthiness without any additional mechanisms in well-mixed finite populations. Even full trust and near full trustworthiness can be the most common strategies. These results are qualitatively different from those of the previous works. Our results thereby demonstrate rich evolutionary dynamics of the asymmetric trust game.

Geometric confinement strongly influences the behavior of microparticles in liquid environments. However, to date, nonspherical particle behaviors close to confining boundaries, even as simple as planar walls, remain largely unexplored. Here, we measure the height distribution and orientation of colloidal dumbbells above walls by means of digital in-line holographic microscopy. We find that while larger dumbbells are oriented almost parallel to the wall, smaller dumbbells of the same material are surprisingly oriented at preferred angles. We determine the total height-dependent force acting on the dumbbells by considering gravitational effects and electrostatic particle-wall interactions. Our modeling reveals that at specific heights both net forces and torques on the dumbbells are simultaneously below the thermal force and energy, respectively, which makes the observed orientations possible. Our results highlight the rich near-wall dynamics of nonspherical particles and can further contribute to the development of quantitative frameworks for arbitrarily shaped microparticle dynamics in confinement.

Adsorption of asymmetric particles or molecules into monolayers is important for many biological and technologically relevant physical systems. In-plane ordering can drastically affect adsorption and phase behavior. In this work, a generalized van der Waals theory previously developed [M. V. Zonta and E. R. Soulé, Phys. Rev. E 100, 062703 (2019)] is used to calculated phase behavior and adsorption isotherms in a system of hard-core rodlike particles with in-plane nematic order, as a function of the model parameters (aspect ratio L/B, isotropic and anisotropic interaction parameters χ and ν, and adsorption constant $K_{\mathrm{ads}}$). For small L/B, isotropic-nematic and/or (depending on χ) isotropic liquid-gas coexistence is observed; as L/B increases, coexistence between two different nematic phases appears at low temperature, and liquid-gas equilibrium ceases to be observed for large enough L/B; this is understood considering that as aspect ratio increases, the range of stability of the nematic phase becomes larger. Adsorption isotherms are found to significantly deviate from Langmuir behavior, and are strongly affected by ordering and interactions (surface density in the adsorbed layer increases as interaction parameters and ordering increase). Phase coexistence is observed as discontinuous transitions in adsorption isotherms, where adsorption-desorption hysteresis cycles are possible.

An external load on a particle packing is distributed internally through a heterogeneous network of particle contacts. This contact force distribution determines the stability of the particle packing and the resulting structure. Here, we investigate the homogeneity of the contact force distribution in packings of highly nonconvex particles both in two-dimensional (2D) and three-dimensional (3D) packings. A recently developed discrete element method is used to model packings of nonconvex particles of varying sphericity. Our results establish that in 3D packings the distribution of the contact forces in the normal direction becomes increasingly heterogeneous with decreasing particle sphericity. However, in 2D packings the contact force distribution is independent of particle sphericity, indicating that results obtained in 2D packings cannot be extrapolated readily to 3D packings. Radial distribution functions show that the crystallinity in 3D packings decreases with decreasing particle sphericity. We link the decreasing homogeneity of the contact force distributions to the decreasing crystallinity of 3D packings. These findings are complementary to the previously observed link between the heterogeneity of the contact force distribution and a decreasing packing crystallinity due to an increasing polydispersity of spherical particles.

In the limit of small inertia, stratification, and advection of density, Ardekani and Stocker [Phys. Rev. Lett. 105, 084502 (2010)] derived the flow due to a point-force and force-dipole placed in a linearly density-stratified fluid. In this limit, these flows also represent the far-field flow due to a towed particle and a neutrally buoyant swimming organism in a stratified fluid. Here, we derive these two far-field flows in the limit of small inertia, stratification but at large advection of density. In both these limits, the flow in a stratified fluid decays rapidly and has closed streamlines but certain symmetries present at small advection are lost at large advection. To illustrate the application of these flows, we use them to calculate the drift induced by a towed drop and a swimming organism, as a means to quantify the mixing caused by them. The drift induced in a stratified fluid is less than that in the homogeneous fluid. A towed drop induces a large drift relative to its own volume at small advection while it induces at least an order of magnitude smaller drift at large advection. On the other hand, a swimming organism induces a large partial drift as compared with its own volume irrespective of the magnitude of advection, unless the stresslet exerted by the swimmer is small. These results are useful in understanding the stratification effects on the drift-based contributions to mixing.