2025-07-11 | | Total: 12
Periodic Boundary Conditions (PBC) introduce well-known lattice artifacts. We present a novel Pseudoperiodic Spherical Boundary Condition (SBC) that is perfectly isotropic. Through detailed comparative simulations, we demonstrate that SBC eliminates the structural and dynamic anisotropy inherent to PBC. For the crowded systems where these artifacts are most prominent, our method is also computationally more efficient than standard Minimum Image Convention implementations. This establishes SBC as a powerful, high-fidelity alternative for simulating isotropic matter.
How internal forces are transduced into motion through soft, fluid membranes remains a fundamental question in the study of active systems. To investigate this coupling, we develop a minimal system consisting of a single ferromagnetic particle encapsulated within a lipid vesicle with controlled membrane composition and phase behavior. An external rotating magnetic field actuates the particle, which rotates and translates along the inner membrane leaflet. This motion generates local slip in the membrane; near a substrate, the slip creates a shear gradient across the lubrication gap that propels the vesicle forward. Propulsion is intermittent and strongest when the particle moves near the vesicle bottom, where stress transmission is most effective. We find that the coupling between internal flows and vesicle motion is highly sensitive to membrane elasticity, excess area, and phase coexistence. Local membrane deformation and flow dissipate part of the stress, limiting the efficiency of force transduction. Additionally, membrane fluctuations and external boundaries reduce particle mobility, and in phase-separated membranes, line tension at domain boundaries deflects the particle and gradually reorients membrane structure. These results demonstrate that lipid membranes not only transmit internal stresses but also remodel themselves in response, actively shaping the dynamics of force transduction and motion in active systems.
The main objective of this work is to shed light on the effect of fiber plasticity on the macroscopic response and domain formation in soft biological composites. This goal is pursued by analyzing the plane-strain response of two-phase laminates. In the context of this problem, the effect of fiber plasticity is accounted for by allowing the elastically stiffer layers (``fiber'' phase) to also yield plastically and by taking the soft layers (``matrix'' phase) to be purely elastic solids. The analysis is carried out at finite elastic and plastic strains, but it is restricted to unidirectional, non-monotonic loading paths, applied by initially increasing the macroscopic stretch along the direction of the layers up to a prescribed maximum value and then decreasing the same stretch down to a minimum value. A simple expression is derived for the critical conditions at which the homogenized behavior of the laminate loses strong ellipticity for the first time along the loading path. The relevance of this result stems from the fact that the loss of macroscopic ellipticity of these composites is known to coincide with the onset of bifurcations of the long-wavelength type. It follows from this result that, just like hyperelastic laminates, elastoplastic laminates may lose macroscopic ellipticity whenever their incremental strength in shear perpendicular to the layers vanishes for the first time. For situations in which loss of macroscopic ellipticity does take place, a corresponding post-bifurcation solution for the homogenized behavior of the laminate is computed. The deformed state of the material described by this solution is characterized by twin lamellar domains that are formed at a length scale much larger than the width of the original, microscopic layers, but still much smaller than the overall dimensions of the macroscopic specimen under consideration.
A central and extensively debated question in glass physics concerns whether a single, growing lengthscale fundamentally controls glassy dynamics, particularly in systems lacking obvious structural motifs or medium range crystalline order (MRCO). In this work, we investigate structural and dynamical lengthscales in supercooled liquids using the Kob Andersen binary Lennard Jones (KALJ) model in two compositions: 80:20 and 60:40. We compute the dynamical lengthscale from displacement displacement correlation functions and observe a consistent growth as temperature decreases. To explore the static counterpart, we use a structural order parameter (SOP) based on the mean field caging potential. While this SOP is known to predict short time dynamics effectively, its bare correlation function reveals minimal spatial growth. Motivated by recent findings that long time dynamics reflect collective rearrangements, we perform spatial coarse-graining of the SOP and identify an optimal lengthscale $L_{max}$ that maximises structure dynamics correlation. We show that the structural correlation length derived from SOP coarse-grained over $L_{max}$ exhibits clear growth with cooling and closely tracks the dynamical lengthscale, especially for A particles in the 80:20 mixture and for both A and B particles in the 60:40 system. Our results reconcile the previously observed absence of static length growth in MRCO-free models like KALJ by highlighting the necessity of intermediate range structural descriptors. Furthermore, we find that the particles with larger structural length growth also correspond to species with latent crystallisation tendencies, suggesting a possible link between structural order, dynamics, and incipient crystallisation.
Elastomeric pressure-sensitive adhesives (PSAs) form adhesive bonds under light pressure. Liquid crystal elastomers (LCEs) are exciting PSA candidates as they can impart both anisotropy and temperature-dependence to adhesion, but the full potential of their anisotropic adhesion is unexplored. Here, identical side-chain LCEs, produced as transparent isotropic or nematic films are investigated; the latter aligned in homeotropic or planar geometries. Their room-temperature adhesion, determined through a 90-degree peel test, is consistent with theoretical predictions and strongest in a planar geometry (peeled parallel to the director) with adhesive force per unit length of 0.67 Nmm-1. In contrast, adhesion of the planar perpendicular, isotropic and homeotropic films is 62.5%, 38.5% and 23.0% lower, respectively. The surface contribution to adhesion is identical for all films, confirming that the variation in adhesion is determined solely by the bulk LCE alignment controlled during film preparation. A temperature-dependent adhesion factor is determined from 0 Celsius to 80 Celsius using dynamic mechanical analysis, and found to be in excellent agreement with the peel data at room temperature. Molecular relaxations active above the glass transition temperature are dominant in determining LCE adhesion. The results show that side-chain LCEs can function as transparent, tunable, broad-temperature smart PSAs
The emergence of hydrodynamic bend instabilities in ordered suspensions of active particles is widely observed across diverse living and synthetic systems, and is considered to be governed by dipolar active stresses generated by the self-propelled particles. Here, using linear stability analyses and numerical simulations, we show that a hydrodynamic bend instability can emerge in the absence of any dipolar active stress and solely due to the self-propulsion force acting on polar active units suspended in an incompressible fluid confined to a substrate. Specifically, we show analytically, and confirm in simulations, that a uniformly ordered state develops bend instability above a critical self-propulsion force. Numerical simulations show that a further increase in the self-propulsion strength leads the system towards a disorderly flow state. The results offer a new route for development of hydrodynamic instabilities in two-dimensional self-propelled materials that are in contact with a substrate, with wide implications in layers of orientationally ordered cells and synthetic active particles.
We study the thermodynamics of continuous transformations between two canonical, cubic network phases of block copolymer melts: double-gyroid, an equilibrium morphology for many systems; and double-diamond, often thought to be a close competitor. We use a strong-segregation approach to compute the free energy of double network morphologies as a function of two structural parameters that convert between two limiting cubic cases; a tetragonal stretch of the unit cell in combination with fusion of pairs of trihedal gyroid nodes into tetrahedral diamond nodes. For the simplest case of conformationally symmetric diblock melts, we find that cubic double-diamond sits at an unstable saddle point that is continuously deformable into the lower free energy gyroid, as well as a second metastable, tetragonal network composed by trihedral nodes. We confirm the broad instability of double-diamond at finite segregation using self-consistent field studies and further show that it derives directly from the entropic free energy cost of chain packing in the tubular domains of tetrahedral nodes. Correspondingly, we demonstrate two factors that quench the entropic cost of packing in the tubular domain -- homopolymer blending and elastic asymmetry between the blocks -- promoting double-diamond to a metastable state by way of free energy barrier that separates it from double-gyroid.
We report the combined results of experiments and molecular dynamics simulations conducted to investigate superheating phenomena in vertically vibrated granular matter. Specifically, we consider a system of cubic particles densely packed in a square-lattice array and subjected to different shaking strengths, denoted by Gamma, approaching a critical value Gamma_c. Below Gamma_c, the superheated crystalline array remains indefinitely stable. Above Gamma_c, it transitions progressively into a granular liquid-like state over a Gamma-dependent timescale tau. We show that while an increase in frictional contacts significantly prolongs the lifetime of the superheated crystalline state, it does not play a major role in the scaling laws governing how that lifetime depends on shaking strength. Our findings also show that the transition from the superheated solid to the liquid state in the vibrated system is well described by a Kolmogorov-Johnson-Mehl-Avrami (KJMA) equation, which is commonly used to model phase transformations in thermal systems.
Accurate computational modeling of damage and fracture remains a central challenge in solid mechanics. The finite element method (FEM) is widely used for numerical modeling of fracture problems; however, classical damage models without gradient regularization yield mesh-dependent and usually inaccurate predictions. The use of gradient damage with FEM improves numerical robustness but introduces significant mathematical and numerical implementation complexities. Physics-informed neural networks (PINNs) can encode the governing partial differential equations, boundary conditions, and constitutive models into the loss functions, offering a new method for fracture modeling. Prior applications of PINNs have been limited to small-strain problems and have incorporated gradient damage formulation without a critical evaluation of its necessity. Since PINNs in their basic form are meshless, this work presents a PINN framework for modeling fracture in elastomers undergoing large deformation without the gradient damage formulation. The PINN implementation here does not require training data and utilizes the collocation method to formulate physics-informed loss functions. We have validated the PINN's predictions for various defect configurations using benchmark solutions obtained from FEM with gradient damage formulation. The crack paths obtained using the PINN are approximately insensitive to the collocation point distribution. This study offers new insights into the feasibility of using PINNs without gradient damage and suggests a simplified and efficient computational modeling strategy for fracture problems. The performance of the PINN has been evaluated through systematic variations in key neural network parameters to provide guidance for future applications. The results provide motivation for extending PINN-based approaches to a broader class of materials and damage models in mechanics.
A long and narrow channel confines disks of two sizes. The disks are randomly agitated in a widened channel under moderate pressure, then jammed according to a tunable protocol. We present exact results that characterize jammed macrostates (volume, entropy, jamming patterns). The analysis divides jammed disk sequences into overlapping tiles out of which statistically interacting quasiparticles are constructed. The fractions of small and large disks are controlled by a chemical potential adapted to configurational statistics of granular matter. The results show regimes for the energy parameters (determined by the jamming protocol) that either enhance or suppress the mixing of disk sizes. Size segregation or size alternation driven by steric forces alone are manifestations of a broken symmetry.
We describe the formation of short-wavelength mesophases in a two-dimensional core-softened particle system. By proposing a series of specific ansatz for each relevant phase, we performed a variational analysis to obtain the ground-state phase diagram. Our results reveal a variety of cluster lattice phases with distinct cluster orientations, alongside traditional two-dimensional Bravais lattices such as square, triangular, oblique, and rectangular structures, as well as other non-Bravais arrangements including honeycomb and kagome phases. We characterize in detail the ground-state phase transitions and identify coexistence regions between competing phases, capturing both first-order and continuous transitions. In addition, we highlight the crucial role of the competing length scales introduced by the hard-core repulsion in shaping the rich landscape of mesophases, emphasizing the interplay between intra-cluster structure and inter-cluster organization. This study provides a systematic framework that could support future investigations of classical thermal melting behavior or quantum phase transitions in similar cluster-forming systems.
We report the flow dynamics of two shear-banding wormlike micellar solutions with distinct surfactant chemistries in a Taylor-Couette (TC) setup following a startup shear. The solutions, formulated with CTAB/NaSal and CPyCl/NaSal, exhibit comparable bulk rheology and equilibrium microstructural properties. By varying the TC gap size, we systematically examine elasticity number over a range of 128000-4470000, while flow curvature spans from 0.022 to 0.171. Under a step shear into the stress plateau, both solutions exhibit a pronounced stress overshoot that intensifies with increasing elasticity number and flow curvature, followed by the development of growing flow heterogeneities. Beyond a critical threshold of elasticity number and curvature, the CTAB/NaSal solution exhibits transient flow reversal, whereas the CPyCl/NaSal solution, despite developing similar heterogeneities, does not undergo flow reversal under any tested conditions. We hypothesize that this difference arises from the stronger electrostatic repulsions between CPyCl head groups due to their larger size, coupled with the higher dissociation rate of Cl counterions. Additionally, the quasi-steady velocity profiles are significantly influenced by elasticity number and flow curvature. Wall slip at the outer cylinder exhibits a non-monotonic dependence on these parameters: it is negligible at low values of elasticity number and flow curvature, peaks at intermediate values, and diminishes at high values. Our findings highlight how surfactant chemistry, elasticity number, and flow curvature control shear-banding flows in wormlike micelles.