Microstructural evolution and reverse flow in shear-banding of entangled polymer melts.

The temporal and spatial evolution of shear banding under startup of shear flow was simulated for highly entangled, linear, monodisperse polyethylene melts of differing molecular weight, C750H1502, C1200H2402, and C3000H6002, using a high-fidelity coarse-grained dissipative particle dynamics method. It was determined that shear stress was dominated by segmental orientation of entangled strands at low shear rates, but at a critical shear rate below the reciprocal of the Rouse time, flow-induced disentanglement resulted in the onset of chain tumbling that reduced the average degree of orientation, leading to a regime of decreasing shear stress, with a commensurate onset of increasing average chain extension imposed by the strong flow kinematics that ultimately drove the steady-state shear stress higher. During startup of shear flow, shear band development began immediately after the maximum in the first normal stress difference, where distinct fast and slow bands formed. The slow bands consisted of relatively entangled and coiled molecules, whereas the fast bands consisted of more disentangled and extended chains that experienced quasiperiodic rotation/retraction cycles. The simulation results often exhibited a generation of temporary reverse flow, in which the local fluid velocity was temporarily opposite to that of the bulk flow direction, at the onset of the shear-banding phenomena; this effect was consistent with earlier experiments and theoretical results. The physical mechanism for the generation of reverse flow during shear-band formation was investigated and found to be related to the recoil of the molecules comprising the slow band. Overall, the phenomenon of shear banding appeared to arise due to flow-induced disentanglement from orientational ordering and segmental stretching that affected individual chains to different degrees, ultimately resulting in regions of relatively coiled and entangled chains that evolved into a slow band, whereas the locally disentangled chains, experiencing quasiperiodic stretch-rotation cycles, formed a fast band. The transitional period resulted in a kinematic instability that generated the temporary reverse-flow phenomenon.

A method for calculating the nonequilibrium entropy of a flowing polymer melt via atomistic simulation.

Nonequilibrium thermodynamics as applied to polymeric liquids is limited by the inability to quantify the configurational entropy. There is no known experimental method to determine it rigorously. Theoretically, entropy is based entirely on the configurational microstate of the material, but for polymer liquids, the number of available configurations is immense and covers long length scales associated with the chain-like nature of the constituent molecules. In principle, however, it should be possible to calculate the entropy from a realistic molecular dynamics simulation that contains positional data for each atomic unit making up the polymer macromolecules. However, there are two challenges in calculating the entropy from an atomistic simulation: it is necessary to relate atomic positions to configurational mesostates, depending on the degree of coarse-graining assumed (if any), and then to entropy, and considerable computational resources are required to determine the three-dimensional probability distribution functions of the configurational mesostates. In this study, a method was developed to calculate nonequilibrium entropy using 3d probability distributions for a linear, entangled polyethylene melt undergoing steady-state shear and elongational flow. An approximate equation expressed in terms of second moments of the 3d distributions was also examined, which turned out to provide almost identical values of entropy as the fully 3d distributions at the mesoscopic level associated with the end-to-end vector of the polymer chains.

A theory for the coexistence of coiled and stretched configurational phases in the extensional flow of entangled polymer melts.

It has recently been demonstrated via nonequilibrium molecular dynamics (NEMD) simulation [M. H. Nafar Sefiddashti, B. J. Edwards, and B. Khomami, J. Chem. Phys. 148, 141103 (2018); Phys. Rev. Lett. 121, 247802 (2018)] that the extensional flow of entangled polymer melts can engender, within a definite strain-rate regime [expressed in terms of the Deborah number (De) based on the Rouse time], the coexistence of separate domains consisting primarily of either coiled or stretched chain-like macromolecules. This flow-induced phase separation results in bimodal configurational distributions, where transitions of individual molecules between the coiled and stretched states occur very slowly by hopping over an apparent energy activation barrier. We demonstrate that the qualitative aspects of this phenomenon can be described via the single-mode Rolie-Poly model including Convective Constraint Release (CCR) and finite extensibility of the chain-like macromolecules. This analysis reveals the physical mechanism for the configurational coexistence, namely, the nonlinear rate of change of the average entropic restoring force of a given entangled chain with extension. Under conditions of significant flow-induced disentanglement, the rate of change of the effective restoring force initially decreases with extension (effective spring softening) and then increases (hardens) as the maximum chain length is approached. When balanced by flow-induced chain stretching, we find that there can be two configuration states within the same De regime, as covered by the NEMD simulations; therefore, a region of conformational coexistence can indeed exist. However, we demonstrate that this coexistence of configurational microstates is only possible when the magnitude of the CCR parameters is consistent with the rate of flow-induced disentanglement, as observed in the NEMD simulations.

Effects of chain length and polydispersity on shear banding in simple shear flow of polymeric melts.

The characteristics of shear banding were investigated in entangled, polydisperse, linear polymer melts under steady-state and startup conditions of simple shear flow. This virtual experimentation was conducted using course-grained nonequilibrium dissipative particle dynamics simulations expressed in terms of a force-field representation that faithfully models the atomistic system dynamics. We examined melts with two mean molecular bead numbers of Nn = 2 50 and 400 and polydispersity indexes of 1.0, 1.025, and 1.05. The wide range of relaxation timescales in the polydisperse melts decreased the nonmonotonic character of the steady-state shear stress vs. shear rate profile compared to a monodisperse linear melt. The polydispersity level required to observe a stress plateau in the shear stress profile at intermediate shear rates was correlated with the nominal entanglement density. Startup of shear flow simulations revealed the development of spatial inhomogeneities and dynamic instabilities in polydisperse fluids containing both monotonic and nonmonotonic shear stress flow curves. Although the shape and duration of instabilities were found to be correlated with the monotonicity of the shear stress profile, the onset and underlying mechanism leading to the formation of shear bands were generally universal. The simulations revealed that perturbations arose soon after the occurrence of a large stress overshoot under startup conditions, and that banded structures stemmed from local reorientation and subsequent deconstruction of the entanglement network. Furthermore, data indicated that the inception of strain localization occurred at shear rates near the reciprocal of the Rouse characteristic timescale, [small gamma, Greek, dot above] > τR-1. Transient shear banding was observed in shorter chain melts undergoing startup of shear flow in which instabilities arose after the appearance of a stress overshoot. These instabilities eventually decayed, but only long after the stresses had attained their steady-state values. The longer chain melt exhibited a shear band structure that remained indefinitely after the stresses had attained steady state.

Configurational Microphase Separation in Elongational Flow of an Entangled Polymer Liquid.

Manufacturing of plastics is typically performed via flow processing of a molten polymeric fluid. Until recently, conventional knowledge has maintained that the deformation of the constituent molecules under flow is homogeneous and obeys Gaussian statistics. In this study via virtual experimentation, an entangled polyethylene melt subjected to planar elongational flow displays an unanticipated microphase separation into a heterogeneous liquid composed of regions of either highly stretched or tightly coiled macromolecules, thus providing a natural realization of a biphasic coil-stretch transition.

Communication: A coil-stretch transition in planar elongational flow of an entangled polymeric melt.

Virtual experimentation of atomistic entangled polyethylene melts undergoing planar elongational flow revealed an amazingly detailed depiction of individual macromolecular dynamics and the resulting effect on bistable configurational states. A clear coil-stretch transition was evident, in much the same form as first envisioned by de Gennes for dilute solutions of high polymers, resulting in an associated hysteresis in the configurational flow profile over the range of strain rates predicted by theory. Simulations conducted at steady state revealed bimodal distribution functions, in which equilibrium configurational states were simultaneously populated by relatively coiled and stretched molecules which could transition from one conformational mode to the other over a relatively long time scale at critical values of strain rates. The implication of such behavior points to a double-well conformational free energy potential with an activation barrier between the two configurational minima.

Improvement of the Tribological Properties of a Lithium-Based Grease by Addition of Graphene.

Graphene with layer number less than seven, prepared by a mechanical exfoliation method, was used as a friction-reducing additive to a lithium-based grease. The graphene was characterized via AFM, TEM, and Raman spectroscopy. The as-prepared graphene had few defects according to the characterization analysis and appeared to be composed primarily of sheets averaging 1-4 atomic layers. When the graphene was added to a lithium-based grease, the lubrication and antiwear properties of the grease were improved, as quantified by friction coefficient and wear scar measurements. The weld point of the lithium-based grease increased proportionally with graphene loading. At the maximum graphene loading tested (2 wt%), the weld point was 1.6 times that of the pure lithium-based grease. Hence the mechanical properties of the graphene sheets played an important role in improving the tribological properties of the grease.

The effect of particle size on the morphology and thermodynamics of diblock copolymer/tethered-particle membranes.

A combination of self-consistent field theory and density functional theory was used to examine the effect of particle size on the stable, 3-dimensional equilibrium morphologies formed by diblock copolymers with a tethered nanoparticle attached either between the two blocks or at the end of one of the blocks. Particle size was varied between one and four tenths of the radius of gyration of the diblock polymer chain for neutral particles as well as those either favoring or disfavoring segments of the copolymer blocks. Phase diagrams were constructed and analyzed in terms of thermodynamic diagrams to understand the physics associated with the molecular-level self-assembly processes. Typical morphologies were observed, such as lamellar, spheroidal, cylindrical, gyroidal, and perforated lamellar, with the primary concentration region of the tethered particles being influenced heavily by particle size and tethering location, strength of the particle-segment energetic interactions, chain length, and copolymer radius of gyration. The effect of the simulation box size on the observed morphology and system thermodynamics was also investigated, indicating possible effects of confinement upon the system self-assembly processes.

Block copolymer morphology formation on topographically complex surfaces: a self-consistent field theoretical study.

A self-consistent field theoretic study is performed to study morphological development of lamellae-forming diblock copolymers on substrates with a well-defined roughness, modeled as trenches of varying depth and width engraved into the substrates. There are three possible lamellar orientations observed: horizontal lamellae, vertical lamellae that are parallel to the trench direction, and vertical lamellae that are perpendicular to the trench direction. Which of these three morphologies formed depends upon the trench width and surface affinity; however, trench depth has a relatively insignificant effect on the morphological development. Therefore, tuning trench width, but not trench depth, should allow for a reduction of the morphological defect density in directed self-assembly of lamellar morphology of diblock copolymers.

A macroscopic model of proton transport through the membrane-ionomer interface of a polymer electrolyte membrane fuel cell.

The membrane-ionomer interface is the critical interlink of the electrodes and catalyst to the polymer electrolyte membrane (PEM); together forming the membrane electrode assembly in current state-of-the-art PEM fuel cells. In this paper, proton conduction through the interface is investigated to understand its effect on the performance of a PEM fuel cell. The water containing domains at this interface were modeled as cylindrical pores/channels with the anionic groups (i.e., -SO(3)(-)) assumed to be fixed on the pore wall. The interactions of each species with all other species and an applied external field were examined. Molecular-based interaction potential energies were computed in a small test element of the pore and were scaled up in terms of macroscopic variables. Evolution equations of the density and momentum of the species (water molecules and hydronium ions) were derived within a framework of nonequilibrium thermodynamics. The resulting evolution equations for the species were solved analytically using an order-of-magnitude analysis to obtain an expression for the proton conductivity. Results show that the conductivity increases with increasing water content and pore radius, and strongly depends on the separation distance between the sulfonate groups and their distribution on the pore wall. It was also determined that the conductivity of two similar pores of different radii in series is limited by the pore with the smaller radius.

A self-consistent field study of diblock copolymer/charged particle system morphologies for nanofiltration membranes.

A combination of self-consistent field theory and density functional theory was used to examine the stable, 3-dimensional equilibrium morphologies formed by diblock copolymers with a tethered nanoparticle attached either between the two blocks or at the end of one of the blocks. Both neutral and interacting particles were examined, with and without favorable/unfavorable energetic potentials between the particles and the block segments. The phase diagrams of the various systems were constructed, allowing the identification of three types of ordered mesophases composed of lamellae, hexagonally packed cylinders, and spheroids. In particular, we examined the conditions under which the mesophases could be generated wherein the tethered particles were primarily located within the interface between the two blocks of the copolymer. Key factors influencing these properties were determined to be the particle position along the diblock chain, the interaction potentials of the blocks and particles, the block copolymer composition, and molecular weight of the copolymer.

Molecular Processes Leading to Shear Banding in Entangled Polymeric Solutions.

The temporal and spatial evolution of shear banding during startup and steady-state shear flow was studied for solutions of entangled, linear, monodisperse polyethylene C3000H6002 dissolved in hexadecane and benzene solvents. A high-fidelity coarse-grained dissipative particle dynamics method was developed and evaluated based on previous NEMD simulations of similar solutions. The polymeric contribution to shear stress exhibited a monotonically increasing flow curve with a broad stress plateau at intermediate shear rates. For startup shear flow, transient shear banding was observed at applied shear rates within the steady-state shear stress plateau. Shear bands were generated at strain values where the first normal stress difference exhibited a maximum, with lifetimes persisting for up to several hundred strain units. During the lifetime of the shear bands, an inhomogeneous concentration distribution was evident within the system, with higher polymer concentration in the slow bands at low effective shear rate; i.e., γ˙<τR-1, and vice versa at high shear rate. At low values of applied shear rate, a reverse flow phenomenon was observed in the hexadecane solution, which resulted from elastic recoil of the molecules within the slow band. In all cases, the shear bands dissipated at high strains and the system attained steady-state behavior, with a uniform, linear velocity profile across the simulation cell and a homogeneous concentration.

Atomistic Simulation of Flow-Induced Microphase Separation and Crystallization of an Entangled Polyethylene Melt Undergoing Uniaxial Elongational Flow and the Role of Kuhn Segment Extension.

Atomistic simulations of the linear, entangled polyethylene C1000H2002 melt undergoing steady-state and startup conditions of uniaxial elongational flow (UEF) over a wide range of flow strength were performed using a united-atom model for the atomic interactions between the methylene groups constituting the polymer macromolecules. Rheological, topological, and microstructural properties of these nonequilibrium viscoelastic materials were computed as functions of strain rate, focusing on regions of flow strength where flow-induced phase separation and flow-induced crystallization were evident. Results of the UEF simulations were compared with those of prior simulations of planar elongational flow, which revealed that uniaxial and planar flows exhibited essentially a universal behavior, although over strain rate ranges that were not completely equivalent. At intermediate flow strength, a purely configurational microphase separation was evident that manifested as a bicontinuous phase composed of regions of highly stretched molecules that enmeshed spheroidal domains of relatively coiled chains. At high flow strength, a flow-induced crystallization (FIC) occurred, producing a semicrystalline material possessing a high degree of crystallinity and primarily a monoclinic lattice structure. This FIC phase formed at a temperature (450 K) high above the quiescent melting point (≈400 K) and remained stable after cessation of flow for temperature at or below 435 K. Careful examination of the Kuhn segments constituting the polymer chains revealed that the FIC phase only formed once the Kuhn segments had become essentially fully extended under the UEF flow field. Thermodynamic properties such as the heat of fusion and heat capacity were estimated from the simulations and found to compare favorably with experimental values.

Morphology tailoring of thin film block copolymers on patterned substrates.

It is well known that chemically patterned substrates can direct the assembly of adsorbed layers or thin films of block copolymers. For a cylinder-forming diblock copolymer on periodically spot-patterned substrates, the morphology of the block copolymer follows the pattern at the substrate; however, with different periodic spacing and spot size of the pattern, novel morphologies can be created. Specifically, we have demonstrated that new morphologies that are absent in the bulk system can be tailored by judiciously varying the mismatch between the width of the pattern and the periodic spacing of the bulk block copolymer, the top surface affinity, and spot size. New morphologies can thus be achieved, such as honeycomb and ring structures, which do not appear in the bulk system. These results demonstrate a promising strategy for fabrication of new nanostructures from chemically patterned substrates.

Why do medical graduates choose rural careers?

INTRODUCTION: This study is based on the metaphor of the 'rural pipeline' into medical practice. The four stages of the rural pipeline are: (1) contact between rural secondary schools and the medical profession; (2) selection of rural students into medical programs; (3) rural exposure during medical training; and (4) measures to address retention of the rural medical workforce. METHODS: Using the rural pipeline template we conducted a literature review, analysed the selection methods of Australian graduate entry medical schools and interviewed 17 interns about their medical career aspirations. LITERATURE REVIEW: The literature was reviewed to assess the effectiveness of selection practices to predict successful gradation and the impact of rural pipeline components on eventual rural practice. Undergraduate academic performance is the strongest predictor of medical course academic performance. The predictive power of interviews is modest. There are limited data on the predictive power of other measures of non-cognitive performance or the content of the undergraduate degree. Prior rural residence is the strongest predictor of choice of a rural career but extended rural exposure during medical training also has a significant impact. The most significant influencing factors are: professional support at national, state and local levels; career pathway opportunities; contentedness of the practitioner's spouse in rural communities; preparedness to adopt a rural lifestyle; educational opportunities for children; and proximity to extended family and social circle. Analysis of selection methods: Staff involved in student selection into 9 Australian graduate entry medical schools were interviewed. Four themes were identified: (1) rurality as a factor in student selection; (2) rurality as a factor in student selection interviews; (3) rural representation on student selection interview panels; (4) rural experience during the medical course. Interns' career intentions: Three themes were identified: (1) the efficacy of the rural pipeline; (2) community connectedness through the rural pipeline; (3) impediments to the effect of the rural pipeline, the most significant being a partner who was not committed to rural life CONCLUSION: Based on the literature review and interviews, 11 strategies are suggested to increase the number of graduates choosing a career in rural medicine, and one strategy for maintaining practitioners in rural health settings after graduation.

Explanation recruits comparison in a category-learning task.

Generating explanations can be highly effective in promoting category learning; however, the underlying mechanisms are not fully understood. We propose that engaging in explanation can recruit comparison processes, and that this in turn contributes to the effectiveness of explanation in supporting category learning. Three experiments evaluated the interplay between explanation and various comparison strategies in learning artificial categories. In Experiment 1, as expected, prompting participants to explain items' category membership led to (a) higher ratings of self-reported comparison processing and (b) increased likelihood of discovering a rule underlying category membership. Indeed, prompts to explain led to more self-reported comparison than did direct prompts to compare pairs of items. Experiment 2 showed that prompts to compare all members of a particular category ("group comparison") were more effective in supporting rule learning than were pairwise comparison prompts. Experiment 3 found that group comparison (as assessed by self-report) partially mediated the relationship between explanation and category learning. These results suggest that one way in which explanation benefits category learning is by inviting comparisons in the service of identifying broad patterns.

Individual Molecular Dynamics of an Entangled Polyethylene Melt Undergoing Steady Shear Flow: Steady-State and Transient Dynamics.

The startup and steady shear flow properties of an entangled, monodisperse polyethylene liquid (C1000H2002) were investigated via virtual experimentation using nonequilibrium molecular dynamics. The simulations revealed a multifaceted dynamical response of the liquid to the imposed flow field in which entanglement loss leading to individual molecular rotation plays a dominant role in dictating the bulk rheological response at intermediate and high shear rates. Under steady shear conditions, four regimes of flow behavior were evident. In the linear viscoelastic regime ( γ ˙ < τ d - 1 ), orientation of the reptation tube network dictates the rheological response. Within the second regime ( τ d - 1 < γ ˙ < τ R - 1 ), the tube segments begin to stretch mildly and the molecular entanglement network begins to relax as flow strength increases; however, the dominant relaxation mechanism in this region remains the orientation of the tube segments. In the third regime ( τ R - 1 < γ ˙ < τ e - 1 ), molecular disentangling accelerates and tube stretching dominates the response. Additionally, the rotation of molecules become a significant source of the overall dynamic response. In the fourth regime ( γ ˙ > τ e - 1 ), the entanglement network deteriorates such that some molecules become almost completely unraveled, and molecular tumbling becomes the dominant relaxation mechanism. The comparison of transient shear viscosity, η + , with the dynamic responses of key variables of the tube model, including the tube segmental orientation, S , and tube stretch, λ , revealed that the stress overshoot and undershoot in steady shear flow of entangled liquids are essentially originated and dynamically controlled by the S x y component of the tube orientation tensor, rather than the tube stretch, over a wide range of flow strengths.

Comparison of the hydration and diffusion of protons in perfluorosulfonic acid membranes with molecular dynamics simulations. scui@utk.edu.

Classical molecular dynamics (MD) simulations were performed to determine the hydrated morphology and hydronium ion diffusion coefficients in two different perfluorosulfonic acid (PFSA) membranes as functions of water content. The structural and transport properties of 1143 equivalent weight (EW) Nafion, with its relatively long perfluoroether side chains, are compared to the short-side-chain (SSC) PFSA ionomer at an EW of 977. The separation of the side chains was kept uniform in both ionomers consisting of -(CF 2) 15- units in the backbone, and the degree of hydration was varied from 5 to 20 weight % water. The MD simulations indicated that the distribution of water clusters is more dispersed in the SSC ionomer, which leads to a more connected water-channel network at the low water contents. This suggests that the SSC ionomer may be more inclined to form sample-spanning aqueous domains through which transport of water and protons may occur. The diffusion coefficients for both hydronium ions and water molecules were calculated at hydration levels of 4.4, 6.4, 9.6, and 12.8 H 2O/SO 3H for each ionomer. When compared to experimental proton diffusion coefficients, this suggests that as the water content is increased the contribution of proton hopping to the overall proton diffusion increases.

A quantum mechanical study of the decomposition of CF3OCF3 and CF3CF2OCF2CF3 in the presence of AlF3.

The effect of AlF3 on the decomposition of CF3OCF3 and CF3CF2OCF2CF3 is investigated using ab initio theory. Previous work by Pancansky et al. [Pacansky, J.; Waltman, R. J. J Fluorine Chem. 1997, 83, 41] showed that AlF3 significantly reduces the activation energy of the decomposition of CF3OCF3 due to the strong electrostatic interaction between the aluminum trifluoride and the reactant. In this work, a new transition-state structure and reaction mechanism have been identified for the decomposition of CF3OCF3 in the presence of AlF3. This new mechanism shows that AlF3 functions by accepting a fluorine atom from one carbon and simultaneously donating a fluorine atom to the other carbon. We show that the same pathway is obtained independently of the level of theory. The reaction rate, generated via statistical mechanics and transition-state theory, is 2-3 orders of magnitude higher for the new transition state when compared to that of the old one. The study was also performed for CF3CF2OCF2CF3 in order to ascertain the effect of chain length on the reaction mechanism and rate. We find that an analogous transition state, with lower activation energy, provides the lowest-energy path for decomposition of the longer chain.

In-plane and out-of-plane rotational motion of individual chain molecules in steady shear flow of polymer melts and solutions.

Recent nonequilibrium molecular dynamics (NEMD) simulations of mildly entangled C400H802 and moderately entangled C700H1402 linear polyethylene melts undergoing steady shear flow have revealed that several inconsistencies between theory and experiment could be rectified by consideration of the rotational motion of individual polymer chains that occurs at moderate to high flow strengths. In this study, we investigated the configurational dynamics of the individual molecular chains that allow these once-entangled, long-chain molecules to execute retraction/extension semi-periodic cycles in response to the imposed shear via NEMD simulations. Brownian dynamics simulations were also performed to extract dynamical and configurational information about the similar cycles of polymer chain behavior that occur in dilute solutions of macromolecular chain liquids dissolved in low molecular weight solvents. Results revealed that the configurational motions of the individual chains in both melt and solution were essentially the same and governed by a single timescale that scaled exponentially with the magnitude of the shear rate. This configurational motion contained both in-plane and out-of-plane components with respect to the flow-gradient plane, with the out-of-plane component playing a much larger role during the retraction phase of the cycle than during the extension phase. This was determined to be caused by the enhancement of the retraction motion by the out-of-plane entropic Brownian forces; however, these entropic forces were detrimental to the in-plane hydrodynamic diffusive forces during the extension phase of the cycle and were thus suppressed. Consequently, the configuration of a rotating chain was significantly more compact during the retraction stage than during the extension stage, wherein the latter phase most molecules were more preferentially distributed in the flow-gradient plane.