Equilibrium molecular dynamics evaluation of the solid-liquid friction coefficient: role of timescales

Solid-liquid friction plays a key role in nanofluidic systems. Yet, despite decades of method development to quantify solid-liquid friction using molecular dynamics (MD) simulations, an accurate and widely applicable method is still missing. Here, we propose a method to quantify the solid-liquid friction coefficient (FC) from equilibrium MD simulations of a liquid confined between parallel solid walls. In this method, the FC is evaluated by fitting the Green-Kubo (GK) integral of the S-L shear force autocorrelation for the range of time scales where the GK integral slowly decays with time. The fitting function was derived based on the analytical solution considering the hydrodynamic equations in our previous work [H. Oga et al., Phys. Rev. Research 3, L032019 (2021)], assuming that the timescales related to the friction kernel and to the bulk viscous dissipation can be separated. By comparing the results with those of other equilibrium MD-based methods and those of non-equilibrium MD for a Lennard-Jones liquid between flat crystalline walls with different wettability, we show that the FC is extracted with excellent accuracy by the present method, even in wettability regimes where other methods become innacurate. We then show that the method is also applicable to grooved solid walls, for which the GK integral displays a complex behavior at short times. Overall, the present method extracts efficiently the FC for various systems, with easy implementation and low computational cost.Comment: 22 pages, 6 Figure

Equilibrium molecular dynamics evaluation of the solid-liquid friction coefficient: Role of timescales

Solid-liquid friction plays a key role in nanofluidic systems. Following the pioneering work of Bocquet and Barrat, who proposed to extract the friction coefficient (FC) from the plateau of the Green-Kubo (GK) integral of the solid-liquid shear force autocorrelation, the so-called plateau problem has been identified when applying the method to finite-sized molecular dynamics simulations, e.g., with a liquid confined between parallel solid walls. A variety of approaches have been developed to overcome this problem. Here, we propose another method that is easy to implement, makes no assumptions about the time dependence of the friction kernel, does not require the hydrodynamic system width as an input, and is applicable to a wide range of interfaces. In this method, the FC is evaluated by fitting the GK integral for the timescale range where it slowly decays with time. The fitting function was derived based on an analytical solution of the hydrodynamics equations [Oga et al., Phys. Rev. Res. 3, L032019 (2021)], assuming that the timescales related to the friction kernel and the bulk viscous dissipation can be separated. By comparing the results with those of other GK-based methods and non-equilibrium molecular dynamics, we show that the FC is extracted with excellent accuracy by the present method, even in wettability regimes where other GK-based methods suffer from the plateau problem. Finally, the method is also applicable to grooved solid walls, where the GK integral displays complex behavior at short times.Oga H., Omori T., Joly L., et al. Equilibrium molecular dynamics evaluation of the solid-liquid friction coefficient: Role of timescales, Journal of Chemical Physics, 159(2), 024701, 14 July 2023, © 2023 American Chemical Society. https://doi.org/10.1063/5.0155628

Shear force measurement of the hydrodynamic wall position in molecular dynamics

Flows in nanofluidic systems are strongly affected by liquid-solid slip, which is quantified by the slip length and by the position where the slip boundary condition applies. Here, we show that the viscosity, slip length, and hydrodynamic wall position (HWP) can be accurately determined from a single molecular dynamics (MD) simulation of a Poiseuille flow, after identifying a relation between the HWP and the wall shear stress in that configuration. From this relation, we deduce that in gravity-driven flows, the HWP identifies with the Gibbs dividing plane of the liquid-vacuum density profile. Simulations of a generic Lennard-Jones liquid confined between parallel frozen walls show that the HWP for a pressure-driven flow is also close to the Gibbs dividing plane (measured at equilibrium), which therefore provides an inexpensive estimate of the HWP, going beyond the common practice of assuming a given position for the hydrodynamic wall. For instance, we show that the HWP depends on the wettability of the surface, an effect usually neglected in MD studies of liquid-solid slip. Overall, the method introduced in this article is simple, fast, and accurate and could be applied to a large variety of systems of interest for nanofluidic applications.Cecilia Herrero, Takeshi Omori, Yasutaka Yamaguchi, and Laurent Joly, "Shear force measurement of the hydrodynamic wall position in molecular dynamics", The Journal of Chemical Physics 151, 041103 (2019) https://doi.org/10.1063/1.5111966

Angiomyofibroblastoma of the vulva: a large pedunculated mass formation.

Angiomyofibroblastoma is a rare, usually small benign mesenchymal tumor that occurs in vulvar lesions of premenopausal women. A case of angiomyofibroblastoma that arose as a unique pedunculated and particularly large mass in the left vulva of a 48-year-old woman is presented herein. The patient had been aware of a gradually enlarged mass of 7 years duration without any other gynecological symptoms or signs. The maximum dimension of the tumor measured 11 cm. The resected tumor was well circumscribed with a bulging and glistening cut surface. Histological examination revealed an admixture of irregularly distributed hypercellular and hypocellular areas with spindled, plump spindled, or plasmacytoid stromal cells and abundant venular or capillary-sized vessels. Stromal cells characteristically cluster around delicate vessels within an edematous to collagenous matrix. In the present case, intralesional adipose tissue was present throughout the tumor. There was no significant nuclear atypia, and mitotic figures were very sparse. There was little stromal mucin throughout the tumor. Immunohistochemically, the stromal cells were characterized by strong reactivity for vimentin and CD34, with focal reactivity for desmin and alpha smooth muscle actin. Both estrogen and progesterone receptors were diffusely expressed in the stromal cells. These histological findings are consistent with angiomyofibroblastoma and support the hypothesis that angiomyofibroblastoma originates from perivascular stem cells with a capacity for myofibroblastic and fatty differentiation

Transplantation of a human induced pluripotent stem cell-derived airway epithelial cell sheet into the middle ear of rats

[Introduction] Early postoperative regeneration of the middle ear mucosa is essential for the prevention of postoperative refractory otitis media and recurrent cholesteatoma. As a means for intractable otitis media management, we focused on human induced pluripotent stem cell (hiPSC)-derived airway epithelial cells (AECs), which have been used in upper airway mucosal regeneration and transplantation therapy. In this study, we transplanted hiPSC-derived AECs into the middle ear of immunodeficient rats. [Methods] Following the preparation of AEC sheets from hiPSCs, the bilateral middle ear mucosa of X-linked severe combined immunodeficient rats was scraped, and the AEC sheets were transplanted in the ears unilaterally. [Results] Human nuclear antigen (HNA)-positive ciliated cells were observed on the transplanted side of the middle ear cavity surface in three of six rats in the 1-week postoperative group and in three of eight rats in the 2-week postoperative group. No HNA-positive cells were found on the control side. The percentage of HNA-positive ciliated cells in the transplanted areas increased in the 2-week postoperative group compared with the 1-week group, suggesting survival of hiPSC-derived AECs. Additionally, HNA-positive ciliated cells were mainly located at sites where the original ciliated cells were localized. Immunohistochemical analysis showed that the transplanted AECs contained cytokeratin 5- and mucin 5AC-positive cells, indicating that both basal cells and goblet cells had regenerated within the middle ear cavity. [Conclusions] The results of this study are an important first step in the establishment of a novel transplantation therapy for chronic otitis media