BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Shengfeng Cheng

Shengfeng Cheng contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mtrl-scicond-mat.softphysics.comp-phBiological PhysicsBiomoleculescond-mat.mes-hallphysics.chem-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2021arXiv

Coarse-Grained Molecular Dynamics Modeling of A Branched Polyetherimide

A coarse-grained model is developed to allow large-scale molecular dynamics (MD) simulations of a branched polyetherimide derived from two backbone monomers [4,4'-bisphenol A dianhydride (BPADA) and m-phenylenediamine (MPD)], a chain terminator [phthalic anhydride (PA)], and a branching agent [tris[4-(4-aminophenoxy)phenyl] ethane (TAPE)]. An atomistic model is first built for the branched polyetherimide. A systematic protocol based on chemistry-informed grouping of atoms, derivation of bond and angle interactions by direct Boltzmann inversion, and parameterization of nonbonded interactions by potential of mean force (PMF) calculations via gas-phase MD simulations of atomic group pairs, is used to construct the coarse-grained model. A six-pair geometry, with one atomic group at the center and six replicates of the other atomic group placed surrounding the central group in a NaCl structure, has been demonstrated to significantly speed up the PMF calculations and partially capture the many-body aspect of the PMFs. Furthermore, we propose a correction term to the PMFs that can make the resulting coarse-grained model transferable temperature-wise, by enabling the model to capture the thermal expansion property of the polymer. The coarse-grained model has been applied to explore the mechanical, structural, and rheological properties of the branched polyetherimide.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Hydration, Ion Distribution, and Ionic Network Formation in Sulfonated Poly(arylene ether sulfones)

We use molecular dynamics simulations to probe hydration, ion spacing, and cation-anion interaction in two sulfonated polysulfones with different ion distributions along the polymer backbone. At room temperature, these polymers remain in a glassy state even with water contents more than 10%. At the equilibrium water uptake, the ions exhibit a similar level of hydration as they would in their saturated aqueous solution. The framework of Manning's limiting law for counterion condensation is used to examine ionic interactions in the simulated polysulfones. The dielectric constant that the ions experience can be well approximated by a volume-weighted average of the dielectric constants of the polymer backbone and water. Our results show that a reasonable estimate of the average inter-ion distance, b, is obtained by using the distance where the sulfonate-sulfonate coordination number reaches 1. The spacing of the sulfonate ions along the polysulfone backbone plays a role in determining their spatial distribution inside the hydrated polymer. As a result, the value of b is slightly larger for polymers where the sulfonate ions are more evenly spaced along the backbone, which is consistent with experimental evidence. The simulations reveal that the sulfonate ions and sodium counterions form fibrillar aggregates at water contents below the equilibrium water uptake. Such extensive ionic aggregates are expected to facilitate ion transport in sulfonated polysulfone membranes, without the need for long-range chain motion as in the case of traditional rubbery ionic polymers. Our estimates for the dielectric constant and b are used in conjunction with Manning's theory to estimate the fraction of counterions condensed to the fixed ions. The prediction of Manning's theory agrees well with the simulation result.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Determination of Glass Transition Temperature of Polyimides from Atomistic Molecular Dynamics Simulations and Machine-Learning Algorithms

Glass transition temperature ($T_{\text{g}}$) plays an important role in controlling the mechanical and thermal properties of a polymer. Polyimides are an important category of polymers with wide applications because of their superior heat resistance and mechanical strength. The capability of predicting $T_{\text{g}}$ for a polyimide $a~priori$ is therefore highly desirable in order to expedite the design and discovery of new polyimide polymers with targeted properties and applications. Here we explore three different approaches to either compute $T_{\text{g}}$ for a polyimide via all-atom molecular dynamics (MD) simulations or predict $T_{\text{g}}$ via a mathematical model generated by using machine-learning algorithms to analyze existing data collected from literature. Our simulations reveal that $T_{\text{g}}$ can be determined from examining the diffusion coefficient of simple gas molecules in a polyimide as a function of temperature and the results are comparable to those derived from data on polymer density versus temperature and actually closer to the available experimental data. Furthermore, the predictive model of $T_{\text{g}}$ derived with machine-learning algorithms can be used to estimate $T_{\text{g}}$ successfully within an uncertainty of about 20 degrees, even for polyimides yet to be synthesized experimentally.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Capillary Adhesion at the Nanometer Scale

Molecular dynamics simulations are used to study the capillary adhesion from a nonvolatile liquid meniscus between a spherical tip and a flat substrate. The atomic structure of the tip, the tip radius, the contact angles of the liquid on the two surfaces, and the volume of the liquid bridge are varied. The capillary force between the tip and substrate is calculated as a function of their separation h. The force agrees with continuum predictions for h down to ~ 5 to 10nm. At smaller h, the force tends to be less attractive than predicted and has strong oscillations. This oscillatory component of the capillary force is completely missed in the continuum theory, which only includes contributions from the surface tension around the circumference of the meniscus and the pressure difference over the cross section of the meniscus. The oscillation is found to be due to molecular layering of the liquid confined in the narrow gap between the tip and substrate. This effect is most pronounced for large tip radii and/or smooth surfaces. The other two components considered by the continuum theory are also identified. The surface tension term, as well as the meniscus shape, is accurately described by the continuum prediction for h down to ~ 1nm, but the capillary pressure term is always more positive than the corresponding continuum result. This shift in the capillary pressure reduces the average adhesion by a factor as large as 2 from its continuum value and is found to be due to an anisotropy in the pressure tensor. The cross-sectional component is consistent with the capillary pressure predicted by the continuum theory (i.e., the Young-Laplace equation), but the normal pressure that determines the capillary force is always more positive than the continuum counterpart.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Molecular Dynamics Simulations of Evaporation-Induced Nanoparticle Assembly

While evaporating solvent is a widely used technique to assemble nano-sized objects into desired superstructures, there has been limited work on how the assembled structures are affected by the physical aspects of the process. We present large scale molecular dynamics simulations of the evaporation-induced assembly of nanoparticles suspended in a liquid that evaporates in a controlled fashion. The quality of the nanoparticle crystal formed just below the liquid/vapor interface is found to be better at relatively slower evaporation rates, as less defects and grain boundaries appear. This trend is understood as the result of the competition between the accumulation and diffusion times of nanoparticles at the liquid/vapor interface. When the former is smaller, nanoparticles are deposited so fast at the interface that they do not have sufficient time to arrange through diffusion, which leads to the prevalence of defects and grain boundaries. Our results have important implications in understanding assembly of nanoparticles and colloids in non-equilibrium liquid environments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Self-assembly of Chiral Tubules

The efficient and controlled assembly of complex structures from macromolecular building blocks is a critical open question in both biological systems and nanoscience. Using molecular dynamics simulations we study the self-assembly of tubular structures from model macromolecular monomers with multiple binding sites on their surfaces [Cheng et al., Soft Matter 8, 5666-5678 (2012)]. In this work we add chirality to the model monomer and a lock-and-key interaction. The self-assembly of free monomers into tubules yields a pitch value that often does not match the chirality of the monomer (including achiral monomers). We show that this mismatch occurs because of a twist deformation that brings the lateral interaction sites into alignment when the tubule pitch differs from the monomer chirality. The energy cost for this deformation is small as the energy distributions substantially overlap for small differences in the pitch and chirality. In order to control the tubule pitch by preventing the twist deformation, the interaction between the vertical surfaces must be increased without resulting in kinetically trapped structures. For this purpose, we employ the lock-and-key interactions and obtain good control of the self-assembled tubule pitch. These results explain some fundamental features of microtubules. The vertical interaction strength is larger than the lateral in microtubules because this yields a controlled assembly of tubules with the proper pitch. We also generally find that the control of the assembly into tubules is difficult, which explains the wide range of pitch and protofilament number observed in microtubule assembly.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Self-assembly of artificial microtubules

Understanding the complex self-assembly of biomacromolecules is a major outstanding question. Microtubules are one example of a biopolymer that possesses characteristics quite distinct from standard synthetic polymers that are derived from its hierarchical structure. In order to understand how to design and build artificial polymers that possess features similar to those of microtubules, we have initially studied the self-assembly of model monomers into a tubule geometry. Our model monomer has a wedge shape with lateral and vertical binding sites that are designed to form tubules. We used molecular dynamics simulations to study the assembly process for a range of binding site interaction strengths. In addition to determining the optimal regime for obtaining tubules, we have calculated a diagram of the structures that form over a wide range of interaction strengths. Unexpectedly, we find that the helical tubules form, even though the monomer geometry is designed for nonhelical tubules. We present the detailed dynamics of the tubule self-assembly process and show that the interaction strengths must be in a limited range to allow rearrangement within clusters. We extended previous theoretical methods to treat our system and to calculate the boundaries between different structures in the diagram.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Structure and Diffusion of Nanoparticle Monolayers Floating at Liquid/Vapor Interfaces: A Molecular Dynamics Study

Large-scale molecular dynamics simulations are used to simulate a layer of nanoparticles diffusing on the surface of a liquid. Both a low viscosity liquid, represented by Lennard-Jones monomers, and a high viscosity liquid, represented by linear homopolymers, are studied. The organization and diffusion of the nanoparticles are analyzed as the nanoparticle density and the contact angle between the nanoparticles and liquid are varied. When the interaction between the nanoparticles and liquid is reduced the contact angle increases and the nanoparticles ride higher on the liquid surface, which enables them to diffuse faster. In this case the short range order is also reduced as seen in the pair correlation function. For the polymeric liquids, the out-of-layer fluctuation is suppressed and the short range order is slightly enhanced. However, the diffusion becomes much slower and the mean square displacement even shows sub-linear time dependence at large times. The relation between diffusion coefficient and viscosity is found to deviate from that in bulk diffusion. Results are compared to simulations of the identical nanoparticles in 2-dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Evaporation of Lennard-Jones Fluids

Evaporation and condensation at a liquid/vapor interface are ubiquitous interphase mass and energy transfer phenomena that are still not well understood. We have carried out large scale molecular dynamics simulations of Lennard-Jones (LJ) fluids composed of monomers, dimers, or trimers to investigate these processes with molecular detail. For LJ monomers in contact with a vacuum, the evaporation rate is found to be very high with significant evaporative cooling and an accompanying density gradient in the liquid domain near the liquid/vapor interface. Increasing the chain length to just dimers significantly reduces the evaporation rate. We confirm that mechanical equilibrium plays a key role in determining the evaporation rate and the density and temperature profiles across the liquid/vapor interface. The velocity distributions of evaporated molecules and the evaporation and condensation coefficients are measured and compared to the predictions of an existing model based on kinetic theory of gases. Our results indicate that for both monatomic and polyatomic molecules, the evaporation and condensation coefficients are equal when systems are not far from equilibrium and smaller than one, and decrease with increasing temperature. For the same reduced temperature $T/T_c$, where $T_c$ is the critical temperature, these two coefficients are higher for LJ dimers and trimers than for monomers, in contrast to the traditional viewpoint that they are close to unity for monatomic molecules and decrease for polyatomic molecules. Furthermore, data for the two coefficients collapse onto a master curve when plotted against a translational length ratio between the liquid and vapor phase.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mechanical response of a self avoiding membrane: fold collisions and the birth of conical singularities

An elastic membrane that is forced to reside in a container smaller than its natural size will deform and, upon further volume reduction, eventually crumple. The crumpled state is characterized by the localization of energy in a complex network of highly deformed crescent-like regions joined by line ridges. Previous studies have focused on the onset of the crumpled state by analyzing the mechanical response and stability of a conical dislocation, while others have simulated the highly packed regime neglecting the importance of the connectivity of the membrane. Here we show, through a combination of experiments, numerical simulations, and analytic approach, that the emergence of new regions of high stretching is a generic outcome when a self avoiding membrane is subject to a severe geometrical constraint. We demonstrate that, at moderate packing fraction, interlayer interactions produce a response equivalent to the one of a thicker membrane that has the shape of the deformed one. Evidence is found that friction plays a key role stabilizing the folded structures.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Defining Contact at the Atomic Scale

Molecular dynamics simulations are used to study different definitions of contact at the atomic scale. The roles of temperature, adhesive interactions and atomic structure are studied for simple geometries. An elastic, crystalline substrate contacts a rigid, atomically flat surface or a spherical tip. The rigid surface is formed from a commensurate or incommensurate crystal or an amorphous solid. Spherical tips are made by bending crystalline planes or removing material outside a sphere. In continuum theory the fraction of atomically flat surfaces that is in contact rises sharply from zero to unity when a load is applied. This simple behavior is surprisingly difficult to reproduce with atomic scale definitions of contact. Due to thermal fluctuations, the number of atoms making contact at any instant rises linearly with load over a wide range of loads. Pressures comparable to the ideal hardness are needed to achieve full contact at typical temperatures. A simple harmonic mean-field theory provides a quantitative description of this behavior and explains why the instantaneous forces on atoms have a universal exponential form. Contact areas are also obtained by counting the number of atoms with a time-averaged repulsive force. For adhesive interactions, the resulting area is nearly independent of temperature and averaging interval, but usually rises from zero to unity over a range of pressures that is comparable to the ideal hardness. The only exception is the case of two identical commensurate surfaces. For nonadhesive surfaces, the mean pressure is repulsive if there is any contact during the averaging interval $Δt$. The associated area is very sensitive to $Δt$ and grows monotonically. Similar complications are encountered in defining contact areas for spherical tips.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Contact and Friction of Nano-Asperities: Effects of Adsorbed Monolayers

Molecular dynamics simulations are used to study contact between a rigid, nonadhesive, spherical tip with radius of order 30nm and a flat elastic substrate covered with a fluid monolayer of adsorbed chain molecules. Previous studies of bare surfaces showed that the atomic scale deviations from a sphere that are present on any tip constructed from discrete atoms lead to significant deviations from continuum theory and dramatic variability in friction forces. Introducing an adsorbed monolayer leads to larger deviations from continuum theory, but decreases the variations between tips with different atomic structure. Although the film is fluid, it remains in the contact and behaves qualitatively like a thin elastic coating except for certain tips at high loads. Measures of the contact area based on the moments or outer limits of the pressure distribution and on counting contacting atoms are compared. The number of tip atoms making contact in a time interval grows as a power of the interval when the film is present and logarithmically with the interval for bare surfaces. Friction is measured by displacing the tip at a constant velocity or pulling the tip with a spring. Both static and kinetic friction rise linearly with load at small loads. Transitions in the state of the film lead to nonlinear behavior at large loads. The friction is less clearly correlated with contact area than load.

0 citations0 reviewsNo code linkedOpen access