BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Gary S. Grest

Gary S. Grest appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softcond-mat.mtrl-scicond-mat.mes-hallphysics.flu-dynBiological Physicschem-phcond-mat.stat-mechmtrl-th

Catalog footprint

What is connected

33works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

33 published item(s)

preprint2026arXiv

Reverse segregation and self-organization in inclined chute flows of bidisperse granular mixtures

In the usual segregation scenario for stable inclined chute flows of bidisperse mixtures of fine and coarse spherical particles, coarse particles rise toward the free surface, forming a coarse-rich region atop the flowing pile. Beyond a threshold coarse-to-fine diameter ratio of approximately 4, conversely, the weight of the coarse particles exceeds the segregation driving forces, causing individual coarse particles to sink within the pile and producing a reversed segregation state. However, an understanding of the collective evolution of the pile structure is still lacking when the particle diameter ratio exceeds 4 {\textit{and}} the coarse particle mass fraction is appreciable. To explore this broadly bidisperse limit, we perform discrete element method simulations considering mean particle diameter ratios of up to 8 and coarse particle mass fractions spanning 0.1 to 0.9. The steady-state flow profiles reveal several intriguing behaviors that depend on the diameter ratio and mass fraction. These include a previously identified transition from usual to reverse segregation and a newfound tendency to self-organize into alternating coarse- and fine-rich particle layers stacked along the shear gradient direction, with layer thickness dictated by the coarse particle diameter. A fuller understanding of segregation at this scale could pave the way for enhanced mixing or demixing techniques at the commercial scale.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Shape Elasticity in Colloidal Bent-Core Liquid Crystals

Curved particles have been shown to stabilize a range of states with unique order in dense suspensions of colloidal bent core liquid crystals. The shape of the colloidal rods encourages the formation of curved director fields. However, states of constant bend cannot uniformly fill either two or three dimensional Euclidean space and are therefore geometrically frustrated. As a result, curved rods are forced to couple their preference for bend with additional twist and splay deformations, giving rise to twist-bend and splay-bend states of nematic and smectic order. In this article, we study the effect of rod curvature on these diverse states of liquid crystalline order using molecular dynamics simulations of a bonded particle model of curved rods with tunable shape elasticity. Focusing on the case of intermediately curved rods, we find that curved rods go through a sequence of isotropic, nematic twist-bend and smectic splay-bend ordering as the density is increased from the dilute limit, in agreement with previous studies of rigid rods. As the rods become more elastic, the critical concentration separating these phases is shifted to higher density. Lastly, we find that flexibility weakens the first-order phase transition separating the isotropic and nematic twist-bend phases.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Flow and Arrest in Stressed Granular Materials

Flowing granular materials often abruptly arrest if not driven by sufficient applied stresses. Such abrupt cessation of motion can be economically expensive in industrial materials handling and processing, and is significantly consequential in intermittent geophysical phenomena such as landslides and earthquakes. Using discrete element simulations, we calculate states of steady flow and arrest for granular materials under the conditions of constant applied pressure and shear stress, which are also most relevant in practice. Here the material can dilate or compact, and flow or arrest, in response to the applied stress. Our simulations highlight that under external stress, the intrinsic response of granular materials is characterized by uniquely-defined steady states of flow or arrest, which are highly sensitive to interparticle friction. While the flowing states can be equivalently characterized by volume fraction, coordination number or internal stress ratio, to characterize the states of shear arrest, one needs to also consider the structural anisotropy in the contact network. We highlight the role of dilation in the flow-arrest transition, and discuss our findings in the context of rheological transitions in granular materials.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fluctuations and power-law scaling of dry, frictionless granular rheology near the hard-particle limit

The flow of frictionless granular particles is studied with stress-controlled discrete element modeling simulations for systems varying in size from 300 to 100,000 particles. The volume fraction and shear stress ratio $μ$ are relatively insensitive to system size fo a wide range of inertial numbers $I$. Second-order effects in strain rate, such as second normal stress differences, require large system sizes to accurately extract meaningful results, notably a non-monotonic dependence in the first normal stress difference with strain rate. The first-order rheological response represented by the $μ(I)$ relationship works well at describing the lower-order aspects of the rheology, except near the quasi-static limit of these stress-controlled flows. The pressure is varied over five decades, and a pressure dependence of the coordination number is observed, which is not captured by the inertial number. Large fluctuations observed for small systems $N\le$ 1,000 near the quasi-static limit can lead to arrest of flow resulting in challenges to fitting the data to rheological relationships. The inertial number is also insufficient for capturing the pressure-dependent behavior of property fluctuations. Fluctuations in the flow and microstructural properties are measured in both the quasi-static and inertial regimes, including shear stress, pressure, strain rate, normal stress differences, volume fraction, coordination number and contact fabric anisotropy. The fluctuations in flow properties scale self-similarly with pressure and system size. A transition in the scaling of fluctuations of stress properties and contact fabric anisotropy are measured and proposed as a quantitative identification of the transition from inertial to quasi-static flow.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Large-scale frictionless jamming with power-law particle size distributions

Due to significant computational expense, discrete element method simulations of jammed packings of size-dispersed spheres with size ratios greater than 1:10 have remained elusive, limiting the correspondence between simulations and real-world granular materials with large size dispersity. Invoking a recently developed neighbor binning algorithm, we generate mechanically-stable jammed packings of frictionless spheres with power-law size distributions containing up to nearly four million particles with size ratios up to 1:100. By systematically varying the width and exponent of the underlying power laws, we analyze the role of particle size distributions on the structure of jammed packings. The densest packings are obtained for size distributions that balance the relative abundance of large-large/intermediate and small-small particle contacts. Although the proportion of rattler particles and mean coordination number strongly depend on the size distribution, the mean coordination of non-rattler particles attains the frictionless isostatic value of six in all cases. The size distribution of non-rattler particles that participate in the load-bearing network exhibits no dependence on the width of the total particle size distribution beyond a critical particle size for low-magnitude exponent power laws. This signifies that only particles with sizes greater than the critical particle size contribute to the mechanical stability. However, for high-magnitude exponent power laws, all particle sizes participate in the mechanical stability of the packing.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Shear Is Not Always Simple: Rate-Dependent Effects of Flow Type on Granular Rheology

Despite there being an infinite variety of types of flow, most rheological studies focus on a single type such as simple shear. Using discrete element simulations, we explore bulk granular systems in a wide range of flow types at large strains and characterize invariants of the stress tensor for different inertial numbers and interparticle friction coefficients. We identify a strong dependence on the type of flow which grows with increasing inertial number or friction. Standard models of yielding, repurposed to describe the dependence of the stress on flow type in steady-state flow and at finite rates, are compared with data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Super-stretchable Elastomer from Cross-linked Ring Polymers

The stretchability of polymeric materials is critical to many applications such as stretchable electronics and soft robotics, yet the stretchability of conventional cross-linked linear polymers is limited by the entanglements between polymer chains. We show using molecular dynamics simulations that cross-linked ring polymers are significantly more stretchable than cross-linked linear polymers. Compared to linear polymers, the entanglements between ring polymers do not act as effective cross-links. As a result, the stretchability of cross-linked ring polymers is determined by the maximum extension of polymer strands between cross-links, rather than between entanglements as in cross-linked linear polymers. The more compact conformation of ring polymers before deformation also contributes to the increase in stretchability.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Evolution of Internal Granular Structure at the Flow-Arrest Transition

The evolution of the internal granular structure in shear-arrested and shear-flowing states of granular materials is characterized using fabric tensors as descriptors of the internal contact and force networks. When a dilute system of frictional grains is subjected to a constant pressure and shear stress, the bulk stress ratio is well-predicted from the anisotropy of its contact and force networks during transient flow prior to steady shear flow or shear arrest. Although the onset of shear arrest is a stochastic process, the fabric tensors upon arrest are distributed around nearly equal contributions of force and contact network anisotropy to the bulk stress ratio. The distribution becomes seemingly narrower with increasing system size. The anisotropy of the contact network in shear-arrested states is reminiscent of the fabric anisotropy observed in shear-jammed packings.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Granular packings with sliding, rolling and twisting friction

Intuition tells us that a rolling or spinning sphere will eventually stop due to the presence of friction and other dissipative interactions. The resistance to rolling and spinning/twisting torque that stops a sphere also changes the microstructure of a granular packing of frictional spheres by increasing the number of constraints on the degrees of freedom of motion. We perform discrete element modeling simulations to construct sphere packings implementing a range of frictional constraints under a pressure-controlled protocol. Mechanically stable packings are achievable at volume fractions and average coordination numbers as low as 0.53 and 2.5, respectively, when the particles experience high resistance to sliding, rolling and twisting. Only when the particle model includes rolling and twisting friction, were experimental volume fractions reproduced.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Topological Linking Drives Anomalous Thickening of Ring Polymers In Weak Extensional Flows

Molecular dynamics simulations confirm recent extensional flow experiments showing ring polymer melts exhibit strong extension-rate thickening of the viscosity at Weissenberg numbers $Wi<<1$. Thickening coincides with the extreme elongation of a minority population of rings that grows with $Wi$. The large susceptibility of some rings to extend is due to a flow-driven formation of topological links that connect multiple rings into supramolecular chains. Links form spontaneously with a longer delay at lower $Wi$ and are pulled tight and stabilized by the flow. Once linked, these composite objects experience larger drag forces than individual rings, driving their strong elongation. The fraction of linked rings generated by flow depends non-monotonically on $Wi$, increasing to a maximum when $Wi\sim1$ before rapidly decreasing when the strain rate approaches the relaxation rate of the smallest ring loops $\sim 1/τ_e$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Crazing of Nanocomposites with Polymer-Tethered Nanoparticles

The crazing behavior of polymer nanocomposites formed by blending polymer grafted nanoparticles with an entangled polymer melt is studied by molecular dynamics simulations. We focus on the three key differences in the crazing behavior of a composite relative to the pure homopolymer matrix, namely, a lower yield stress, a smaller extension ratio and a grafted chain length dependent failure stress. The yield behavior is found to be mostly controlled by the local nanoparticle-grafted polymer interfacial energy, with the grafted polymer-polymer matrix interfacial structure being of little to no relevance. Increasing the attraction between nanoparticle core and the grafted polymer inhibits void nucleation and leads to a higher yield stress. In the craze growth regime, the presence of grafted chain sections of 100 monomers alters the mechanical response of composite samples, giving rise to smaller extension ratios and higher drawing stresses than for the homopolymer matrix. The dominant failure mechanism of composite samples depends strongly on the length of the grafted chains, with disentanglement being the dominant mechanism for short chains, while bond breaking is the failure mode for chain lengths greater than 10Ne, where Ne is the entanglement length.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dispersing Nanoparticles in a Polymer Film via Solvent Evaporation

Large scale molecular dynamics simulations are used to study the dispersion of nanoparticles (NPs) in a polymer film during solvent evaporation. As the solvent evaporates, a dense polymer-rich skin layer forms at the liquid/vapor interface, which is either NP rich or poor depending on the strength of the NP/polymer interaction. When the NPs are strongly wet by the polymer, the NPs accumulate at the interface and form layers. However when the NPs are only partially wet by the polymer, most NPs are uniformly distributed in the bulk of the polymer film with the dense skin layer serving as a barrier to prevent the NPs from moving to the interface. Our results point to a possible route to employ less favorable NP/polymer interactions and fast solvent evaporation to uniformly disperse NPs in a polymer film, contrary to the common belief that strong NP/polymer attractions are needed to make NPs well dispersed in polymer nanocomposites.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Polymer Entanglement Dynamics: Role of Attractive Interactions

The coupled dynamics of entangled polymers which span a broad time and length scales govern the unique viscoelastic properties of polymers. To follow chain mobility by numerical simulations from the intermediate Rouse and reptation regimes to the late time diffusive regime, highly coarse grained models with purely repulsive interactions between monomers are widely used since they are computationally the most efficient. Here using large scale molecular dynamics simulations the effect of including the attractive interaction between monomers on the dynamics of entangled polymers melts is explored for the first time over wide temperature range. Attractive interactions has little effect on the local packing for all temperatures T and on the chain mobility for T higher than about twice the glass transition Tg. These results, across a broad range of molecular weight, show that to study the dynamics of entangled polymers melts the interactions can be treated as pure repulsive, confirming a posteriori the validity of previous studies and opening the way to new large scale numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Resolving Dynamic Properties of Polymers through Coarse-Grained Computational Studies

Coupled length and time scales determine the dynamic behavior of polymers and underlie their unique viscoelastic properties. To resolve the long-time dynamics it is imperative to determine which time and length scales must be correctly modeled. Here we probe the degree of coarse graining required to simultaneously retain significant atomistic details and access large length and time scales. The degree of coarse graining in turn sets the minimum length scale instrumental in defining polymer properties and dynamics. Using linear polyethylene as a model system, we probe how coarse graining scale affects the measured dynamics. Iterative Boltzmann inversion is used to derive coarse-grained potentials with 2-6 methylene groups per coarse-grained bead from a fully atomistic melt simulation. We show that atomistic detail is critical to capturing large scale dynamics. Using these models we simulate polyethylene melts for times over 500 μs to study the viscoelastic properties of well-entangled polymer chains.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Diffusion in jammed particle packs

Using random walk simulations we explore diffusive transport through monodisperse sphere packings over a range of packing fractions, $ϕ$, in the vicinity of the jamming transition at $ϕ_{c}$. Various diffusion properties are computed over several orders of magnitude in both time and packing pressure. Two well-separated regimes of normal, "Fickian" diffusion, where the mean squared displacement is linear in time, are observed. The first corresponds to diffusion inside individual spheres, while the latter is the long-time bulk diffusion. The intermediate anomalous diffusion regime and the long-time value of the diffusion coefficient are both shown to be controlled by particle contacts, which in turn depend on proximity to $ϕ_{c}$. The time required to recover normal diffusion $t^*$ scales as $(ϕ-ϕ_c)^{-0.5}$ and the long-time diffusivity $D_\infty \sim (ϕ- ϕ_c)^{0.5}$, or $D_\infty \sim 1/t^*$. It is shown that the distribution of mean first passage times associated with the escape of random walkers between neighboring particles controls both $t^*$ and $D_\infty$ in the limit $ϕ\rightarrow ϕ_c$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Coating thickness and coverage effects on the forces between silica nanoparticles in water

The structure and interactions of coated silica nanoparticles have been studied in water using molecular dynamics simulations. For 5 nm diameter amorphous silica nanoparticles we studied the effects of varying the chain length and grafting density of polyethylene oxide (PEO) on the nanoparticle coating's shape and on nanoparticle-nanoparticle effective forces. For short ligands of length $n=6$ and $n=20$ repeat units, the coatings are radially symmetric while for longer chains ($n=100$) the coatings are highly anisotropic. This anisotropy appears to be governed primarily by chain length, with coverage playing a secondary role. For the largest chain lengths considered, the strongly anisotropic shape makes fitting to a simple radial force model impossible. For shorter ligands, where the coatings are isotropic, we found that the force between pairs of nanoparticles is purely repulsive and can be fit to the form $(R/2r_\text{core}-1)^{-b}$ where $R$ is the separation between the center of the nanoparticles, $r_\text{core}$ is the radius of the silica core, and $b$ is measured to be between 2.3 and 4.1.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Healing of polymer interfaces: Interfacial dynamics, entanglements, and strength

Self-healing of polymer films often takes place as the molecules diffuse across a damaged region, above their melting temperature. Using molecular dynamics simulations we probe the healing of polymer films and compare the results with those for thermal welding of homopolymer slabs. The two processes differ in their interfacial structure since damage leads to increased polydispersity and more short chains. A polymer sample was cut into two separate films that were then held together in the melt state. The recovery of the damaged film was followed as time elapsed and polymer molecules diffused across the interface. The mass uptake and formation of entanglements, as obtained from primitive path analysis, are extracted and correlated with the interfacial strength obtained from shear simulations. We find that the interdiffusion is significantly faster in the damaged film compared to welding because of the presence of short chains. Though interfacial entanglements increase more rapidly for the damaged films, a large fraction of them are near chain ends. As a result, the interfacial strength of the healing film increases more slowly than for welding. For both healing and welding, the interfacial strength saturates as the bulk entanglement density is recovered across the interface. However, the saturation strength of the damaged film is below the bulk strength for the polymer sample. At saturation, cut chains remain near the healing interface. They are less entangled and thus mechanically weaken the interface. When the interfacial strength saturates, the number of interfacial entanglements scales with the corresponding bulk entanglement density. Chain stiffness increases the density of entanglements, which increases the interfacial strength. Our results show that a few entanglements across the interface are sufficient to resist interfacial chain pullout and enhance the mechanical strength.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Tensile Fracture of Welded Polymer Interfaces: Miscibility, Entanglements and Crazing

Large-scale molecular simulations are performed to investigate tensile failure of polymer interfaces as a function of welding time $t$. Changes in the tensile stress, mode of failure and interfacial fracture energy $G_I$ are correlated to changes in the interfacial entanglements as determined from Primitive Path Analysis. Bulk polymers fail through craze formation, followed by craze breakdown through chain scission. At small $t$ welded interfaces are not strong enough to support craze formation and fail at small strains through chain pullout at the interface. Once chains have formed an average of about one entanglement across the interface, a stable craze is formed throughout the sample. The failure stress of the craze rises with welding time and the mode of craze breakdown changes from chain pullout to chain scission as the interface approaches bulk strength. The interfacial fracture energy $G_I$ is calculated by coupling the simulation results to a continuum fracture mechanics model. As in experiment, $G_I$ increases as $t^{1/2}$ before saturating at the average bulk fracture energy $G_b$. As in previous simulations of shear strength, saturation coincides with the recovery of the bulk entanglement density. Before saturation, $G_I$ is proportional to the areal density of interfacial entanglements. Immiscibiltiy limits interdiffusion and thus suppresses entanglements at the interface. Even small degrees of immisciblity reduce interfacial entanglements enough that failure occurs by chain pullout and $G_I \ll G_b$.

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

Structure and Strength at Immiscible Polymer Interfaces

Thermal welding of polymer-polymer interfaces is important for integrating polymeric elements into devices. When two different polymers are joined, the strength of the weld depends critically on the degree of immiscibility. We perform large-scale molecular dynamics simulations of the structure-strength relation at immiscible polymer interfaces. Our simulations show that immiscibility arrests interdiffusion and limits the equilibrium interfacial width. Even for weakly immiscible films, the narrow interface is unable to transfer stress upon deformation as effectively as the bulk material, and chain pullout at the interface becomes the dominant failure mechanism. This greatly reduces the interfacial strength. The weak response of immiscible interfaces is shown to arise from an insufficient density of entanglements across the interface. We demonstrate that there is a threshold interfacial width below which no significant entanglements can form between opposite sides to strengthen the interface.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Polymer Welding: Strength Through Entanglements

Large-scale simulations of thermal welding of polymers are performed to investigate the rise of mechanical strength at the polymer-polymer interface with the welding time. The welding process is in the core of integrating polymeric elements into devices as well as in thermal induced healing of polymers; processes that require development of interfacial strength equal to that of the bulk. Our simulations show that the interfacial strength saturates at the bulk shear strength much before polymers diffuse by their radius of gyration. Along with the strength increase, the dominant failure mode changes from chain pullout at the interface to chain scission as in the bulk. Formation of sufficient entanglements across the interface, which we track using a Primitive Path Analysis is required to arrest catastrophic chain pullout at the interface. The bulk response is not fully recovered until the density of entanglements at the interface reaches the bulk value. Moreover, the increase of interfacial strength before saturation is proportional to the number of interfacial entanglements between chains from opposite sides.

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

Molecular Dynamics Simulation Study of Nonconcatenated Ring Polymers in a Melt: I. Statics

Molecular dynamics simulations were conducted to investigate the structural properties of melts of nonconcatenated ring polymers and compared to melts of linear polymers. The longest rings were composed of N=1600 monomers per chain which corresponds to roughly 57 entanglement lengths for comparable linear polymers. For the rings, the radius of gyration squared was found to scale as N to the 4/5 power for an intermediate regime and N to the 2/3 power for the larger rings indicating an overall conformation of a crumpled globule. However, almost all beads of the rings are "surface beads" interacting with beads of other rings, a result also in agreement with a primitive path analysis performed in the following paper (DOI: 10.1063/1.3587138). Details of the internal conformational properties of the ring and linear polymers as well as their packing are analyzed and compared to current theoretical models.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Molecular Dynamics Simulation Study of Nonconcatenated Ring Polymers in a Melt: II. Dynamics

Molecular dynamics simulations were conducted to investigate the dynamic properties of melts of nonconcatenated ring polymers and compared to melts of linear polymers. The longest rings were composed of N=1600 monomers per chain which corresponds to roughly 57 entanglement lengths for comparable linear polymers. The ring melts were found to diffuse faster than their linear counterparts, with the self-diffusion coefficient for both architectures scaling as approximately N to the -2.4 power for large N. The mean-square displacement of the center-of-mass of the rings follows a sub-diffusive behavior for times and distances beyond the mean-square gyration radius, neither compatible with the Rouse nor the reptation model. The rings relax stress much faster than linear polymers and the zero-shear viscosity was found to vary as approximately N to the 1.4 power which is much weaker than the N to the 3.4 power of linear chains, not matching any commonly known model for polymer dynamics when compared to the observed mean-square displacements. These findings are discussed in view of the conformational properties of the rings presented in the preceding paper (DOI: 10.1063/1.3587137).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Rheology of Ring Polymer Melts: From Linear Contaminants to Ring/Linear Blends

Ring polymers remain a major challenge to our current understanding of polymer dynamics. Experimental results are difficult to interpret because of the uncertainty in the purity and dispersity of the sample. Using both equilibrium and non-equilibrium molecular dynamics simulations we have systematically investigated the structure, dynamics and rheology of perfectly controlled ring/linear polymer blends with chains of such length and flexibility that the number of entanglements is up to about 14 per chain, which is comparable to experimental systems examined in the literature. The smallest concentration at which linear contaminants increase the zero-shear viscosity of a ring polymer melt of these chain lengths by 10% is approximately one-fifth of their overlap concentration. When the two architectures are present in equal amounts the viscosity of the blend is approximately twice as large as that of the pure linear melt. At this concentration the diffusion coefficient of the rings is found to decrease dramatically, while the static and dynamic properties of the linear polymers are mostly unaffected. Our results are supported by a primitive path analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Effect of Particle Shape and Charge on Bulk Rheology of Nanoparticle Suspensions

The rheology of nanoparticle suspensions for nanoparticles of various shapes with equal mass is studied using molecular dynamics simulations. The equilibrium structure and the response to imposed shear are analyzed for suspensions of spheres, rods, plates, and jacks in an explicit solvent for both charged and uncharged nanoparticles. For the volume fraction studied, ?$ϕ_{vf}=0.075$, the uncharged systems are all in their isotropic phase and the viscosity is only weakly dependent on shape for spheres, rods, and plate whereas for the jacks the viscosity is an order of magnitude larger than for the other three shapes. The introduction of charge increases the viscosity for all four nanoparticle shapes with the increase being the largest for rods and plates. The presence of a repulsive charge between the particles decreases the amount of stress reduction that can be achieved by particle reorientation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

End Grafted Polymer Nanoparticles in a Polymeric Matrix: Effect of Coverage and Curvature

It has recently been proposed that the miscibility of nanoparticles with a polymer matrix can be controlled by grafting polymer chains to the nanoparticle surface. As a first step to study this situation, we have used molecular dynamics simulations on a single nanoparticle of radius R ($4σ\le$R$\le 16σ$, where $σ$ is the diameter of a polymer monomer) grafted with chains of length 500 in a polymer melt of chains of length 1000. The grafting density $Σ$ was varied between $0.04$-$0.32$ chains/$σ^2$. To facilitate equilibration a Monte Carlo double-bridging algorithm is applied - new bonds are formed across a pair of chains, creating two new chains each substantially different from the original. For the long brush chains studied here, the structure of the brush assumes its large particle limit even for $R$ as small as 8$σ$, which is consistent with recent experimental findings. We study autophobic dewetting of the melt from the brush as a function of increasing $Σ$. Even these long brush and matrix chains of length $6$ and $12$ $N_e$, respectively, (the entanglement length is $N_e \sim 85$) give somewhat ambiguous results for the interfacial width, showing that studies of two or more nanoparticles are necessary to properly understand these miscibility issues. Entanglement between the brush and melt chains were identified using the primitive path analysis. We find that the number of entanglements between the brush and melt chains scale simply with the product of the local monomer densities of brush and melt chains.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stress Relaxation in Entangled Polymer Melts

We present an extensive set of simulation results for the stress relaxation in equilibrium and step-strained bead-spring polymer melts. The data allow us to explore the chain dynamics and the shear relaxation modulus, $G(t)$, into the plateau regime for chains with $Z=40$ entanglements and into the terminal relaxation regime for $Z=10$. Using the known (Rouse) mobility of unentangled chains and the melt entanglement length determined via the primitive path analysis of the microscopic topological state of our systems, we have performed parameter -free tests of several different tube models. We find excellent agreement for the Likhtman-McLeish theory using the double reptation approximation for constraint release, if we remove the contribution of high-frequency modes to contour length fluctuations of the primitive chain.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Forces between functionalized silica nanoparticles in solution

To prevent the flocculation and phase separation of nanoparticles in solution, nanoparticles are often functionalized with short chain surfactants. Here we present fully-atomistic molecular dynamics simulations which characterize how these functional coatings affect the interactions between nanoparticles and with the surrounding solvent. For 5 nm diameter silica nanoparticles coated with poly(ethylene oxide) (PEO) oligomers in water, we determined the hydrodynamic drag on two approaching nanoparticles moving through solvent and on a single nanoparticle as it approaches a planar surface. In most circumstances, acroscale fluid theory accurately predicts the drag on these nano-scale particles. Good agreement is seen with Brenner's analytical solutions for wall separations larger than the soft nanoparticle radius. For two approaching coated nanoparticles, the solvent-mediated (velocity-independent) and lubrication (velocity-dependent) forces are purely repulsive and do not exhibit force oscillations that are typical of uncoated rigid spheres.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Three-dimensional shear in granular flow

The evolution of granular shear flow is investigated as a function of height in a split-bottom Couette cell. Using particle tracking, magnetic-resonance imaging, and large-scale simulations we find a transition in the nature of the shear as a characteristic height $H^*$ is exceeded. Below $H^*$ there is a central stationary core; above $H^*$ we observe the onset of additional axial shear associated with torsional failure. Radial and axial shear profiles are qualitatively different: the radial extent is wide and increases with height while the axial width remains narrow and fixed.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Large-scale simulation of adhesion dynamics for end-graphed polymers

The adhesion between a polymer melt and substrate is studied in the presence of chemically attached chains on the substrate surface. Extensive molecular dynamics simulations have been carried out to study the effect of temperature, tethered chain areal density ($Σ$), tethered chain length ($N_{t}$), chain bending energy ($k_θ$) and tensile pull velocity ($v$) on the adhesive failure mechanisms of pullout and/or scission of the tethered chains. We observe a crossover from pure chain pullout to chain scission as $N_{t}$ is increased. Below the glass transition, the value of $N_{t}$ for which this crossover begins approaches the bulk entanglement length $N_{e}$. For the values of $N_{t}$ and $Σ$ used here, no crossover to crazing is observed.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Phase Coexistence of a Stockmayer Fluid in an Applied Field

We examine two aspects of Stockmayer fluids which consists of point dipoles that additionally interact via an attractive Lennard-Jones potential. We perform Monte Carlo simulations to examine the effect of an applied field on the liquid-gas phase coexistence and show that a magnetic fluid phase does exist in the absence of an applied field. As part of the search for the magnetic fluid phase, we perform Gibbs ensemble simulations to determine phase coexistence curves at large dipole moments, $μ$. The critical temperature is found to depend linearly on $μ^2$ for intermediate values of $μ$ beyond the initial nonlinear behavior near $μ=0$ and less than the $μ$ where no liquid-gas phase coexistence has been found. For phase coexistence in an applied field, the critical temperatures as a function of the applied field for two different $μ$ are mapped onto a single curve. The critical densities hardly change as a function of applied field. We also verify that in an applied field the liquid droplets within the two phase coexistence region become elongated in the direction of the field.

0 citations0 reviewsNo code linkedOpen access