Local stress fields are routinely computed from molecular dynamics trajectories to understand the structure and mechanical properties of lipid bilayers. These calculations can be systematically understood with the Irving-Kirkwood-Noll theory. In identifying the stress tensor, a crucial step is the decomposition of the forces on the particles into pairwise contributions.However, such a decomposition is not unique in general, leading to an ambiguity in the definition of the stress tensor, particularly for multibody potentials. Furthermore, a theoretical treatment of constraints in local stress calculations has been lacking. Here, we present a new implementation of local stress calculations that systematically treats constraints and considers a privileged decomposition, the central force decomposition, that leads to a symmetric stress tensor by construction. We focus on biomembranes, although the methodology presented here is widely applicable. Our results show that some unphysical behavior obtained with previous implementations, e.g. non-constant normal stress profiles along an isotropic bilayer in equilibrium, is a consequence of an improper treatment of constraints. Furthermore, other valid force decompositions produce significantly different stress profiles, particularly in the presence of * To whom correspondence should be addressed † Universitat Politècnica de Catalunya-BarcelonaTech ‡ Contributed equally to this work 1 dihedral potentials. Our methodology reveals the striking effect of unsaturations on the bilayer mechanics, missed by previous stress calculation implementations.
Fluid deformable surfaces are ubiquitous in cell and tissue biology, including lipid bilayers, the actomyosin cortex, or epithelial cell sheets. These interfaces exhibit a complex interplay between elasticity, low Reynolds number interfacial hydrodynamics, chemistry, and geometry, and govern important biological processes such as cellular traffic, division, migration, or tissue morphogenesis. To address the modelling challenges posed by this class of problems, in which interfacial phenomena tightly interact with the shape and dynamics of the surface, we develop a general continuum mechanics and computational framework for fluid deformable surfaces. The dual solid-fluid nature of fluid deformable surfaces challenges classical Lagrangian or Eulerian descriptions of deforming bodies. Here, we extend the notion of Arbitrarily Lagrangian-Eulerian (ALE) formulations, wellestablished for bulk media, to deforming surfaces. To systematically develop models for fluid deformable surfaces, which consistently treat all couplings between fields and geometry, we follow a nonlinear Onsager formalism according to which the dynamics minimize a Rayleighian functional where dissipation, power input and energy release rate compete. Finally, we propose new computational methods, which build on Onsager's formalism and our ALE formulation, to deal with the resulting stiff system of higher-order of partial differential equations. We apply our theoretical and computational methodology to classical models for lipid bilayers and the cell cortex. The methods developed here allow us to formulate/simulate these models for the first time in their full three-dimensional generality, accounting for finite curvatures and finite shape changes.
The microscopic stress field provides a unique connection between atomistic simulations and mechanics at the nanoscale. However, its definition remains ambiguous. Rather than a mere theoretical preoccupation, we show that this fact acutely manifests itself in local stress calculations of defective graphene, lipid bilayers, and fibrous proteins. We find that popular definitions of the microscopic stress violate the continuum statements of mechanical equilibrium, and we propose an unambiguous and physically sound definition.
We revisit the derivation of the microscopic stress, linking the statistical mechanics of particle systems and continuum mechanics. The starting point in our geometric derivation is the Doyle-Ericksen formula, which states that the Cauchy stress tensor is the derivative of the free-energy with respect to the ambient metric tensor and which follows from a covariance argument. Thus, our approach to define the microscopic stress tensor does not rely on the statement of balance of linear momentum as in the classical Irving-Kirkwood-Noll approach. Nevertheless, the resulting stress tensor satisfies balance of linear and angular momentum. Furthermore, our approach removes the ambiguity in the definition of the microscopic stress in the presence of multibody interactions by naturally suggesting a canonical and physically motivated force decomposition into pairwise terms, a key ingredient in this theory. As a result, our approach provides objective expressions to compute a microscopic stress for a system in equilibrium and for force-fields expanded into multibody interactions of arbitrarily high order. We illustrate the proposed
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.