Source author record

Mark O. Robbins

Mark O. Robbins appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

20 published item(s)

preprint2020arXiv

Effect of flow-induced molecular alignment on welding and strength of polymer interfaces

Structures formed by fused filament fabrication are often substantially weaker than those made with conventional techniques, and fail at the welds between successive layers. One factor that may influence strength is flow-induced alignment of deposited material. Recent work suggests that alignment reduces the entanglement density and thus should accelerate welding by diffusion. Here, coarse-grained molecular simulations are used to test the effect of molecular alignment on diffusion and weld strength. While standard measures show a decrease of the entanglement density with alignment, there is no change in the rate of diffusion normal to the interface or the rate of formation of entanglements across the interface. The time for chain reorientation also remains equal to the equilibrium disentanglement time $τ_d$. Despite this, simulations of mechanical tests show that welds formed from aligned states are weaker until several $τ_d$. This is not because the weld itself is weaker, but because aligned material near the weld is weaker than unaligned material. The maximum shear strength and tensile fracture energy of welded systems are the same as bulk systems with the same alignment.

preprint2016arXiv

Contact area of rough spheres: Large scale simulations and simple scaling laws

We use molecular simulations to study the nonadhesive and adhesive atomic-scale contact of rough spheres with radii ranging from nanometers to micrometers over more than ten orders of magnitude in applied normal load. At the lowest loads, the interfacial mechanics is governed by the contact mechanics of the first asperity that touches. The dependence of contact area on normal force becomes linear at intermediate loads and crosses over to Hertzian at the largest loads. By combining theories for the limiting cases of nominally flat rough surfaces and smooth spheres, we provide parameter-free analytical expressions for contact area over the whole range of loads. Our results establish a range of validity for common approximations that neglect curvature or roughness in modeling objects on scales from atomic force microscope tips to ball bearings.

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.

preprint2016arXiv

Nanocapillary adhesion between parallel plates

Molecular dynamics simulations are used to study capillary adhesion from a nanometer scale liquid bridge between two parallel flat solid surfaces. The capillary force and the meniscus shape of the bridge are computed as the separation between the solid surfaces is varied. Macroscopic theory predicts the meniscus shape and the contribution of liquid/vapor interfacial tension to the capillary force quite accurately for separations as small as 2 or 3 molecular diameters (1-2nm). However the total capillary force differs in sign and magnitude from macroscopic theory for separations less than about 5nm (8-10 diameters) because of molecular layering that is not included in macroscopic theory. For these small separations, the pressure tensor in the fluid becomes anisotropic. The components in the plane of the surface vary smoothly and are consistent with theory based on the macroscopic surface tension. Capillary adhesion is affected by only the perpendicular component, which has strong oscillations as the molecular layering changes.

preprint2015arXiv

Chain Ends and the Ultimate Strength of Polyethylene Fibers

We use large scale molecular dynamics (MD) simulations to determine the tensile yield mechanism of orthorhombic polyethylene (PE) crystals with finite chains spanning $10^2-10^4$ carbons in length. We find the yield stress $σ_y$ saturates for long chains at 6.3 GPa, agreeing well with experiments. We show chains do not break but always yield by slip, after nucleation of 1D dislocations at chain ends. Dislocations are accurately described by a Frenkel-Kontorova model parametrized by the mechanical properties of an ideal crystal. We compute a dislocation core size $ξ\approx25$Å and determine the high and low strain rate limits of $σ_y$. Our results suggest characterizing the 1D dislocations of polymer crystals as an efficient method for numerically predicting the ultimate tensile strength of aligned fibers.

preprint2015arXiv

The effect of inertia on sheared disordered solids: Critical scaling of avalanches in two and three dimensions

Molecular dynamics simulations with varying damping are used to examine the effects of inertia and spatial dimension on sheared disordered solids in the athermal, quasistatic limit. In all cases the distribution of avalanche sizes follows a power law over at least three orders of magnitude in dissipated energy or stress drop. Scaling exponents are determined using finite-size scaling for systems with thousands to millions of particles. Three distinct universality classes are identified corresponding to overdamped and underdamped limits, as well as a crossover damping that separates the two regimes. For each universality class, the exponent describing the avalanche distributions is the same in two and three dimensions. The spatial extent of plastic damage is proportional to the energy dissipated in an avalanche. Both rise much more rapidly with system size in the underdamped limit where inertia is important. Inertia also lowers the mean energy of configurations sampled by the system and leads to an excess of large events like that seen in earthquake distributions for individual faults. The distribution of stress values during shear narrows to zero with increasing system size and may provide useful information about the size of elemental events in experimental systems. For overdamped and crossover systems the stress variation scales inversely with the square root of the system size. For underdamped systems the variation is determined by the size of the largest events.

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.

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.

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$.

preprint2013arXiv

Finite-size scaling in the interfacial stiffness of rough elastic contacts

The total elastic stiffness of two contacting bodies with a microscopically rough interface has an interfacial contribution K that is entirely attributable to surface roughness. A quantitative understanding of K is important because it can dominate the total mechanical response and because it is proportional to the interfacial contributions to electrical and thermal conductivity in continuum theory. Numerical simulations of the dependence of K on the applied squeezing pressure p are presented for nominally flat elastic solids with a range of surface roughnesses. Over a wide range of p, K rises linearly with p. Sublinear power-law scaling is observed at small p, but the simulations reveal that this is a finite-size effect. We derive accurate, analytical expressions for the exponents and prefactors of this low-pressure scaling of K by extending the contact mechanics theory of Persson to systems of finite size. In agreement with our simulations, these expressions show that the onset of the low-pressure scaling regime moves to lower pressure as the system size increases.

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.

preprint2013arXiv

When are rough surfaces sticky?

At the molecular scale there are strong attractive interactions between surfaces, yet few macroscopic surfaces are sticky. Extensive simulations of contact by adhesive surfaces with roughness on nanometer to micrometer scales are used to determine how roughness reduces the area where atoms contact and thus weakens adhesion. The material properties, adhesive strength and roughness parameters are varied by orders of magnitude. In all cases the area of atomic contact rises linearly with load, and the prefactor rises linearly with adhesive strength for weak interactions. Above a threshold adhesive strength, the prefactor changes sign, the surfaces become sticky and a finite force is required to separate them. A parameter-free analytic theory is presented that describes changes in these numerical results over up to five orders of magnitude in load. It relates the threshold strength to roughness and material properties, explaining why most macroscopic surfaces do not stick. The numerical results are qualitatively and quantitatively inconsistent with classical theories based on the Greenwood-Williamson approach that neglect the range of adhesion and do not include asperity interactions.

preprint2012arXiv

Avalanches in Strained Amorphous Solids: Does Inertia Destroy Critical Behavior?

Simulations are used to determine the effect of inertia on athermal shear of a two-dimensional binary Lennard-Jones glass. In the quasistatic limit, shear occurs through a series of rapid avalanches. The distribution of avalanches is analyzed using finite-size scaling with thousands to millions of particles. Inertia takes the system to a new underdamped universality class rather than driving the system away from criticality. Scaling exponents are determined for the underdamped and overdamped limits and a critical damping that separates the two regimes. Systems are in the overdamped universality class even when most vibrational modes are underdamped.

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.

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.

preprint2010arXiv

Nanoindentation of virus capsids in a molecular model

A molecular-level model is used to study the mechanical response of empty cowpea chlorotic mottle virus (CCMV) and cowpea mosaic virus (CPMV) capsids. The model is based on the native structure of the proteins that consitute the capsids and is described in terms of the C-alpha atoms. Nanoindentation by a large tip is modeled as compression between parallel plates. Plots of the compressive force versus plate separation for CCMV are qualitatively consistent with continuum models and experiments, showing an elastic region followed by an irreversible drop in force. The mechanical response of CPMV has not been studied, but the molecular model predicts an order of magnitude higher stiffness and a much shorter elastic region than for CCMV. These large changes result from small structural changes that increase the number of bonds by only 30% and would be difficult to capture in continuum models. Direct comparison of local deformations in continuum and molecular models of CCMV shows that the molecular model undergoes a gradual symmetry breaking rotation and accommodates more strain near the walls than the continuum model. The irreversible drop in force at small separations is associated with rupturing nearly all of the bonds between capsid proteins in the molecular model while a buckling transition is observed in continuum models.

preprint2010arXiv

Stiffness of Contacts Between Rough Surfaces

The effect of self-affine roughness on solid contact is examined with molecular dynamics and continuum calculations. The contact area and normal and lateral stiffnesses rise linearly with the applied load, and the load rises exponentially with decreasing separation between surfaces. Results for a wide range of roughnesses, system sizes and Poisson ratios can be collapsed using Persson's contact theory for continuous elastic media. The atomic scale response at the interface between solids has little affect on the area or normal stiffness, but can greatly reduce the lateral stiffness. The scaling of this effect with system size and roughness is discussed.

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.

preprint2009arXiv

Strain hardening in bidisperse polymer glasses: Separating the roles of chain orientation and interchain entanglement

The effects of entanglement and chain orientation on strain hardening in glassy polymers are separated by examining mixtures of chains with different lengths. Simulations show that the orientation of a molecule of a given chain length is the same in monodisperse systems and bidisperse mixtures, even when entangled and unentangled chains are mixed. In addition, the stress in mixtures is equal to the weighted average of the stresses in monodisperse systems. These results indicate that chains contribute independently to strain hardening, that chain orientation is determined by local interactions with the surrounding glass, and that entanglements play at most an indirect role in strain hardening in the range of strains typically studied. We discuss these results in the context of recent microscopic theories.