Structure of Strongly Adsorbed Polymer Systems: A Computer Simulation Study

The structure of very thin polymer films formed by strongly adsorbed macromolecules was studied by computer simulation. A coarse-grained model of strictly two-dimensional polymer systems was built, and its properties determined by an efficient Monte Carlo simulation algorithm. Properties of the model system were determined by means of Monte Carlo simulations with a sampling algorithm that combines Verdier-Stockmayer, pivot and reputation moves. The effects of temperature, chain length and polymer concentration on the macromolecular structure were investigated. It was shown that at low temperatures, the chain size increases with the concentration, that is, inversely with high temperatures. This behavior should be explained by the influence of inter-chain interactions.

Fractional Coupling of Primary and Johari-Goldstein Relaxations in a Model Polymer

A polymer model exhibiting heterogeneous Johari−Goldstein (JG) secondary relaxation is studied by extensive molecular-dynamics simulations of states with different temperature and pressure. Time−temperature−pressure superposition of the primary (segmental) relaxation is evidenced. The time scales of the primary and the JG relaxations are found to be highly correlated according to a power law. The finding agrees with key predictions of the Coupling Model (CM) accounting for the decay in a correlation function due to the relaxation and diffusion of interacting systems. Nonetheless, the exponent of the power law, even if it is found in the range predicted by CM (0<ξ<1), deviates from the expected one. It is suggested that the deviation could depend on the particular relaxation process involved in the correlation function and the heterogeneity of the JG process.

Molecular Insights into the Wall Slip Behavior of Pseudoplastic Polymer Melt in Nanochannels during Micro Injection Molding

Wall slip directly affects the molding quality of plastic parts by influencing the stability of the filling flow field during micro injection molding. The accurate modeling of wall slip in nanochannels has been a great challenge for pseudoplastic polymer melts. Here, an effective modeling method for polymer melt flow in nanochannels based on united-atom molecular dynamics simulations is presented. The effects of driving forces and wall–fluid interactions on the behavior of polyethylene melt under Poiseuille flow conditions were investigated by characterizing the slip velocity, dynamics information of the flow process, and spatial configuration parameters of molecular chains. The results indicated that the united-atom molecular dynamics model could better describe the pseudoplastic behavior in nanochannels than the commonly used finitely extensible nonlinear elastic (FENE) model. It was found that the slip velocity could be increased with increasing driving force and show completely opposite variation trends under different orders of magnitude of the wall–fluid interactions. The influence mechanism was interpreted by the density distribution and molecular chain structure parameters, including disentanglement and orientation, which also coincides with the change in the radius of gyration.

Non-Isothermal Free-Surface Viscous Flow of Polymer Melts in Pipe Extrusion Using an Open-Source Interface Tracking Finite Volume Method

Polymer extrudate swelling is a rheological phenomenon that occurs after polymer melt flow emerges at the die exit of extrusion equipment due to molecular stress relaxations and flow redistributions. Specifically, with the growing demand for large scale and high productivity, polymer pipes have recently been produced by extrusion. This study reports the development of a new incompressible non-isothermal finite volume method, based on the Arbitrary Lagrangian–Eulerian (ALE) formulation, to compute the viscous flow of polymer melts obeying the Herschel–Bulkley constitutive equation. The Papanastasiou-regularized version of the constitutive equation is employed. The influence of the temperature on the rheological behavior of the material is controlled by the Williams–Landel–Ferry (WLF) function. The new method is validated by comparing the extrudate swell ratio obtained for Bingham and Herschel–Bulkley flows (shear-thinning and shear-thickening) with reference data found in the scientific literature. Additionally, the essential flow characteristics including yield-stress, inertia and non-isothermal effects were investigated.

New Considerations on the Determination of the Apparent Shear Viscosity of Polymer Melt with Micro Capillary Dies

Capillary rheometers have been widely used for the rheological measurement of polymer melts. However, when micro capillary dies are used, the results are usually neither accurate nor consistent, even under the same measurement conditions. In this work, theoretical modeling and experimental studies were conducted for a more profound understanding of the mechanism by which the initial and boundary conditions influence the inaccuracy in the apparent shear viscosity determination with micro capillary dies (diameters: 500 μm, 200 μm, 100 μm). The results indicate that the amount of polymer initially in the barrel, the pre-compaction pressure and the capillary die diameter have a significant influence on the development of the micro scale inlet pressure, which directly determines the accuracy of the measurement at low and medium shear rates. The varying melt compressibility was confirmed to be the main factor directly related to the inaccuracy in the micro scale apparent shear viscosity determination. It is suggested that measures such as reducing the amount of polymer initially in the barrel and increasing the pre-compaction pressure could be used to reduce the measurement inaccuracy.

Retrieving Equivalent Shear Viscosity for Molten Polymers from 3-D Nonisothermal Capillary Flow Simulation

For highly viscous polymer melts, considerable fluid temperature rises produced by viscous heating can be a disturbing factor in viscosity measurements. By scrutinizing the experimental and simulated capillary pressure losses for polymeric liquids, we demonstrate the importance of applying a viscous heating correction to the shear viscosity, so as to correct for large errors introduced by the undesirable temperature rises. Specifically, on the basis of a theoretical derivation and 3-D nonisothermal flow simulation, an approach is developed for retrieving the equivalent shear viscosity in capillary rheometry, and we show that the shear viscosity can be evaluated by using the average fluid temperature at the wall, instead of the bulk temperature, as previously assumed. With the help of a viscous Cross model in analyzing the shear-dominated capillary flow, it is possible to extract the viscous heating contribution to capillary pressure loss, and the general validity of the methodology is assessed using the experiments on a series of thermoplastic melts, including polymers of amorphous, crystalline, and filler-reinforced types. The predictions of the viscous model based on the equivalent viscosity are found to be in good to excellent agreement with experimental pressure drops. For all the materials studied, a near material-independent scaling relation between the dimensionless temperature rise (Θ) and the Nahme number (Na) is found, Θ ~ Na^0.72, from which the fluid temperature rise due to viscous heating as well as the resultant viscosity change can be predicted.

Highlighting Thermo-Elastic Effects in Confined Fluids

The recent identification of a finite shear elasticity in mesoscopic fluids has motivated the search of other solid-like properties of liquids. We present an innovative thermal approach of liquids. We identify a dynamic thermo-elastic mesoscopic behavior by building the thermal image produced by different liquids upon applying a low frequency mechanical shear field. We selected three fluids: a low molecular weight polybutylacrylate (PBuA), polypropyleneglycol (PPG), and glycerol. We demonstrate that a part of the energy of the shear strain is converted in cold and hot shear bands varying synchronously with the applied shear field. This thermodynamic change suggests a coupling to shear elastic modes in agreement with the low frequency shear elasticity theoretically foreseen and experimentally demonstrated.

Multiscale modelling of heterogeneous fillers in polymer composites: the case of polyisoprene and carbon black

The dispersion of inorganic particles within polymeric materials is an extensively used method to enhance their mechanical properties. One of the major challenges in the simulation of polymer composites is to model the uneven surface of the fillers which strongly affects the dynamics of the adsorbed polymers and consequently the macroscopic mechanical properties of the final composite. Here we propose a new multiscale approach that, using experimental adsorption data, constructs the filler surface to statistically reproduce the surface defects. We use this approach to analyse the structure and dynamics of highly entangled polyisoprene melt in contact with different realistic carbon black samples. We show that the presence of the heterogeneous surface has a negligible influence on the structure of the polymer chains but a major effect on their dynamics and the surface wettability.

Bonding widths of Deposited Polymer Strands in Additive Manufacturing

In this study, we explore the deformation of a polymer extrudate upon the deposition on a build platform, to determine the bonding widths between stacked strands in fused-filament fabrication. The considered polymer melt has an extremely high viscosity, which dominates in its deformation. Mainly considering the viscous effect, we derive analytical expressions of the flat width, compressed depth, bonding width and cross-sectional profile of the filament in four special cases, which have different combinations of extrusion speed, print speed and nozzle height. We further validate the derived relations, using our experimental results on acrylonitrile butadiene styrene (ABS), as well as existing experimental and numerical results on ABS and polylactic acid (PLA). Compared with existing theoretical and numerical results, our derived analytic relations are simple, which need less calculations. They can be used to quickly predict the geometries of the deposited strands, including the bonding widths.

Thermo-mechanical processing of plant proteins using shear cell and high-moisture extrusion cooking

Consumption of plant-based meat analogues offers a way to reduce the environmental footprint of the human diet. High-moisture extrusion cooking (HMEC) and shear cell processing both rely on thermo-mechanical treatment of proteins to product fibrous meat-like products. However, the mechanisms underlying these processes are not well understood. In this review we discuss the effect of thermo-mechanical processing on the physicochemical properties and phase behavior of proteins and protein mixtures. The HMEC and shear cell processes are comparable in their basic unit operations, which are (1) mixing and hydration, (2) thermo-mechanical treatment, and (3) cooling. An often overlooked part of the extruder that could be crucial to fibrillation is the so-called breaker plate, which is situated between the barrel and die sections. We found a lack of consensus on the effect of heat on protein-protein interactions, and that the experimental tools to study protein-protein interactions are limited. The different mechanisms for structure formation proposed in literature all consider the deformation and alignment of the melt. However, the mechanisms differ in their underlying assumptions. Further investigation using novel and dedicated tools is required to fully understand these thermo-mechanical processes.

The Use of Scattering Data in the Study of the Molecular Organisation of Polymers in the Non-Crystalline State

Scattering data for polymers in the non-crystalline state, i.e., the glassy state or the molten state, may appear to contain little information. In this work, we review recent developments in the use of scattering data to evaluate in a quantitative manner the molecular organization of such polymer systems. The focus is on the local structure of chain segments, on the details of the chain conformation and on the imprint the inherent chemical connectivity has on this structure. We show the value of tightly coupling the scattering data to atomistic-level computer models. We show how quantitative information about the details of the chain conformation can be obtained directly using a model built from definitions of relatively few parameters. We show how scattering data may be supplemented with data from specific deuteration sites and used to obtain information hidden in the data. Finally, we show how we can exploit the reverse Monte Carlo approach to use the data to drive the convergence of the scattering calculated from a 3d atomistic-level model with the experimental data. We highlight the importance of the quality of the scattering data and the value in using broad Q scattering data obtained using neutrons. We illustrate these various methods with results drawn from a diverse range of polymers.

The Structure of Dipolar Polymer Brushes and Their Interaction in the Melt. Impact of Chain Stiffness

By using the numerical lattice Scheutjens-Fleer self-consistent field (SF-SCF) method we have studied the effect of the restricted flexibility of grafted chains on the structure and mutual interaction of two opposing planar conventional and A-type dipolar brushes. Brushes are immersed in the solvent consisting of chains similar to the grafted ones. The increase of the chain rigidity enhances the segregation of grafted chains in a A-type dipolar brush into two populations: backfolded chains with terminal monomers near the grafting surface and chains with the ends at the brush periphery. The fraction of backfolded chains grows by an increase of the Kuhn segment length. It is shown that two opposite A-type dipolar brushes from semi-rigid chains are attracted to each other at short distances. The attraction becomes more pronounced and begins at larger distances for more rigid chains with the same brush characteristics: polymerization degree, grafting density, and dipole moments of monomer units. This attraction is connected with the dipole-dipole interactions between chains of oncoming brushes with oppositely directed dipoles penetrating deeply into each other upon contact. This effect of the chain rigidity is opposite to that for conventional brushes without dipoles in the chains. For such brushes, an increase in the chain rigidity leads to the enhanced repulsion between them.

Violation of Stokes-Einstein and Stokes-Einstein-Debye relations in polymers at the gas-supercooled liquid coexistence

Molecular dynamics simulations are performed on a system of model linear polymers to look at the violations of Stokes-Einstein (SE) and Stokes-Einstein-Debye (SED) relations near the mode coupling theory transition temperatureT_cat three (one higher and two lower) densities. At low temperatures, both lower density systems show stable gas-supercooled-liquid coexistence whereas the higher density system is homogeneous. We show that monomer density relaxation exhibits SE violation for all three densities, whereas molecular density relaxation shows a weak violation of the SE relation nearT_cin both lower density systems. This study identifies disparity in monomer mobility and observation of jumplike motion in the typical monomer trajectories resulting in the SE violations. In addition to the SE violation, a weak SED violation is observed in the gas-supercooled-liquid coexisting domains of the lower densities. Both lower density systems also show a decoupling of translational and rotational dynamics in this polymer system.

Coincident Correlation between Vibrational Dynamics and Primary Relaxation of Polymers with Strong or Weak Johari-Goldstein Relaxation

The correlation between the vibrational dynamics, as sensed by the Debye-Waller factor, and the primary relaxation in the presence of secondary Johari-Goldstein (JG) relaxation, has been investigated through molecular dynamics simulations. Two melts of polymer chains with different bond length, resulting in rather different strength of the JG relaxation are studied. We focus on the bond-orientation correlation function, exhibiting higher JG sensitivity with respect to alternatives provided by torsional autocorrelation function and intermediate scattering function. We find that, even if changing the bond length alters both the strength and the relaxation time of the JG relaxation, it leaves unaffected the correlation between the vibrational dynamics and the primary relaxation. The finding is in harmony with previous studies reporting that numerical models not showing secondary relaxations exhibit striking agreement with experimental data of polymers also where the presence of JG relaxation is known.

Crystallisation in Melts of Short, Semi-Flexible Hard-Sphere Polymer Chains: The Role of the Non-Bonded Interaction Range

A melt of short semi-flexible polymers with hard-sphere-type non-bonded interaction undergoes a first-order crystallisation transition at lower density than a melt of hard-sphere monomers or a flexible hard-sphere chain. In contrast to the flexible hard-sphere chains, the semi-flexible ones have an intrinsic stiffness energy scale, which determines the natural temperature scale of the system. In this paper, we investigate the effect of weak additional non-bonded interaction on the phase transition temperature. We study the system using the stochastic approximation Monte Carlo (SAMC) method to estimate the micro-canonical entropy of the system. Since the density of states in the purely hard-sphere non-bonded interaction case already covers 5600 orders of magnitude, we consider the effect of weak interactions as a perturbation. In this case, the system undergoes the same ordering transition with a temperature shift non-uniformly depending on the additional interaction. Short-range attractions impede ordering of the melt of semi-flexible polymers and decrease the transition temperature, whereas relatively long-range attractions assist ordering and shift the transition temperature to higher values, whereas weak repulsive interactions demonstrate an opposite effect on the transition temperature.

Assessment of the Tumbling-Snake Model against Linear and Nonlinear Rheological Data of Bidisperse Polymer Blends

We have recently solved the tumbling-snake model for concentrated polymer solutions and entangled melts in the academic case of a monodisperse sample. Here, we extend these studies and provide the stationary solutions of the tumbling-snake model both analytically, for small shear rates, and via Brownian dynamics simulations, for a bidisperse sample over a wide range of shear rates and model parameters. We further show that the tumbling-snake model bears the necessary capacity to compare well with available linear and non-linear rheological data for bidisperse systems. This capacity is added to the already documented ability of the model to accurately predict the shear rheology of monodisperse systems.

Tumbling-Snake Model for Polymeric Liquids Subjected to Biaxial Elongational Flows with a Focus on Planar Elongation

We have recently solved the tumbling-snake model for concentrated polymer solutions and entangled melts in the presence of both steady-state and transient shear and uniaxial elongational flows, supplemented by a variable link tension coefficient. Here, we provide the transient and stationary solutions of the tumbling-snake model under biaxial elongation both analytically, for small and large elongation rates, and via Brownian dynamics simulations, for the case of planar elongational flow over a wide range of rates, times, and the model parameters. We show that both the steady-state and transient first planar viscosity predictions are similar to their uniaxial counterparts, in accord with recent experimental data. The second planar viscosity seems to behave in all aspects similarly to the shear viscosity, if shear rate is replaced by elongation rate.

Surface Segregation of Cyclic Chains in Binary Melts of Thin Polymer Films: The Influence of Constituent Concentration

We carry out extensive molecular dynamics simulations of thin films of bead-spring models of binary mixtures composed of cyclic and linear polymer chains. We study the equilibrium behavior of the polymer chains for two very different chain lengths, which resemble short (10-mers) and long (100-mers) chains, at different concentrations of the binary mixture. We clearly show how the concentration variable affects the enrichment of either of the two polymer species at the interface, and also how the chain length influences this process.

Monitoring Nanoscale Deformations in a Drawn Polymer Melt with Single-Molecule Fluorescence Polarization Microscopy

Elongating a polymer melt causes polymer segments to align and polymer coils to deform along the drawing direction. Despite the importance of this molecular response for understanding the viscoelastic properties and relaxation behavior of polymeric materials, studies on the single-molecule level are rare and were not performed in real time. Here we use single-molecule fluorescence polarization microscopy for monitoring the position and orientation of single fluorescent perylene diimide molecules embedded in a free-standing thin film of a polymethyl acrylate (PMA) melt with a time resolution of 500 ms during the film drawing and the subsequent stress relaxation period. The orientation distribution of the perylene diimide molecules is quantitatively described with a model of rod-like objects embedded in a uniaxially elongated matrix. The orientation of the fluorescent probe molecules is directly coupled to the local deformation of the PMA melt, which we derive from the distances between individual dye molecules. In turn, the fluorescence polarization monitors the shape deformation of the polymer coils on a length scale of 5 nm. During stress relaxation, the coil shape relaxes four times more slowly than the mechanical stress. This shows that stress relaxation involves processes on length scales smaller than a polymer coil. Our work demonstrates how optical spectroscopy and microscopy can be used to study the coupling of individual fluorescent probe molecules to their embedding polymeric matrix and to an external mechanical stimulus on the single-molecule level.