Search for: All records

Far-from-equilibrium phenomena are critical to all natural and engineered systems, and essential to biological processes responsible for life. For over a century and a half, since Carnot, Clausius, Maxwell, Boltzmann, and Gibbs, among many others, laid the foundation for our understanding of equilibrium processes, scientists and engineers have dreamed of an analogous treatment of nonequilibrium systems. But despite tremendous efforts, a universal theory of nonequilibrium behavior akin to equilibrium statistical mechanics and thermodynamics has evaded description. Several methodologies have proved their ability to accurately describe complex nonequilibrium systems at the macroscopic scale, but their accuracy and predictive capacity is predicated on either phenomenological kinetic equations fit to microscopic data or on running concurrent simulations at the particle level. Instead, we provide a novel framework for deriving stand-alone macroscopic thermodynamic models directly from microscopic physics without fitting in overdamped Langevin systems. The only necessary ingredient is a functional form for a parameterized, approximate density of states, in analogy to the assumption of a uniform density of states in the equilibrium microcanonical ensemble. We highlight this framework’s effectiveness by deriving analytical approximations for evolving mechanical and thermodynamic quantities in a model of coiled-coil proteins and double-stranded DNA, thus producing, to the authors’ knowledge, the first derivation of the governing equations for a phase propagating system under general loading conditions without appeal to phenomenology. The generality of our treatment allows for application to any system described by Langevin dynamics with arbitrary interaction energies and external driving, including colloidal macromolecules, hydrogels, and biopolymers. 

“Viscosity is the most ubiquitous dissipative mechanical behavior” (Maugin, 1999). Despite its ubiquity, even for those systems where the mechanisms causing viscous and other forms of dissipation are known there are only a few quantitative models that extract the macroscopic rheological response from these microscopic mechanisms. One such mechanism is the stochastic breaking and forming of bonds which is present in polymer networks with transient cross-links, strong inter-layer bonding between graphene sheets, and sliding dry friction. In this paper we utilize a simple yet flexible model to show analytically how stochastic bonds can induce an array of rheological behaviors at the macroscale. We find that varying the bond interactions induces a Maxwell-type macroscopic material behavior with Newtonian viscosity, shear thinning, shear thickening, or solid like friction when subjected to shear at constant rates. When bond rupture is independent of the force applied, Newtonian viscosity is the predominant behavior. When bond breaking is accelerated by the applied force, a shear thinning response becomes most prevalent. Further connections of the macroscopic response to the interaction potential and rates of bonding and unbonding are illustrated through phase diagrams and analysis of limiting cases. Finally, we apply this model to polymer networks and to experimental data on “solid bridges” in polydisperse granular media. We imagine possible applications to material design through engineering bonds with specific interactions to bring about a desired macroscopic behavior. 

Self-assembly of proteins on lipid membranes underlies many important processes in cell biology, such as, exo- and endo-cytosis, assembly of viruses, etc. An attractive force that can cause self-assembly is mediated by membrane thickness interactions between proteins. The free energy profile associated with this attractive force is a result of the overlap of thickness deformation fields around the proteins which can be calculated from the solution of a boundary value problem. Yet, the time scales over which two inclusions coalesce has not been explored, even though the evolution of particle concentrations on membranes has been modeled using phase-field approaches. In this paper we compute this time scale as a function of the initial distance between two inclusions by viewing their coalescence as a first passage time problem. The mean first passage time is computed using Langevin dynamics and a partial differential equation, and both methods are found to be in excellent agreement. Inclusions of three different shapes are studied and it is found that for two inclusions separated by about hundred nanometers the time to coalescence is hundreds of milliseconds irrespective of shape. An efficient computation of the interaction energy of inclusions is central to our work. We compute it using a finite difference technique and show that our results are in excellent agreement with those from a previously proposed semi-analytical method based on Fourier–Bessel series. The computational strategies described in this paper could potentially lead to efficient methods to explore the kinetics of self-assembly of proteins on lipid membranes. 

Abstract The mechanism of the scoliotic curve development in healthy adolescents remains unknown in the field of orthopedic surgery. Variations in the sagittal curvature of the spine are believed to be a leading cause of scoliosis in this patient population. Here, we formulate the mechanics of S-shaped slender elastic rods as a model for pediatric spine under physiological loading. Second, applying inverse mechanics to clinical data of the subtypes of scoliotic spines, with characteristic 3D deformity, we determine the undeformed geometry of the spine before the induction of scoliosis. Our result successfully reproduces the clinical data of the deformed spine under varying loads, confirming that the prescoliotic sagittal curvature of the spine impacts the 3D loading that leads to scoliosis. 

The double-helical topology of DNA molecules observed at room temperature in the absence of any external loads can be disrupted by increasing the bath temperature or by applying tensile forces, leading to spontaneous strand separation known as DNA melting. Here, continuum mechanics of a 2D birod is combined with statistical mechanics to formulate a unified framework for studying both thermal melting and tensile force induced melting of double-stranded molecules: it predicts the variation of melting temperature with tensile load, provides a mechanics-based understanding of the cooperativity observed in melting transitions, and reveals an interplay between solution electrostatics and micromechanical deformations of DNA which manifests itself as an increase in the melting temperature with increasing ion concentration. This novel predictive framework sheds light on the micromechanical aspects of DNA melting and predicts trends that were observed experimentally or extracted phenomenologically using the Clayperon equation.