Interaction between plates in a polymer melt

Structure and dynamics of a polymer melt at an attractive surface

We study the structural and dynamic properties of a polymer melt in the vicinity of an adhesive solid substrate by means of Molecular Dynamics simulation at various degrees of surface adhesion. The properties of the individual polymer chains are examined as a function of the distance to the interface and found to agree favorably with theoretical predictions. Thus, the adsorbed amount at the adhesive surface is found to scale with the macromolecule length as Γ is proportional to √N, regardless of the adsorption strength. For chains within the range of adsorption we analyze in detail the probability size distributions of the various building blocks: loops, tails and trains, and find that loops and tails sizes follow power laws while train lengths decay exponentially thus confirming some recent theoretical results. The chain dynamics as well as the monomer mobility are also investigated and found to depend significantly on the proximity of a given layer to the solid adhesive surface with onset of vitrification for sufficiently strong adsorption.

Crystal nucleation enhanced at the diffuse interface of immiscible polymer blends

We report dynamic Monte Carlo simulations of immiscible binary polymer blends, which exhibit weakly enhanced crystal nucleation near interfaces between two phase-separated polymers. We found that this enhancement is not accompanied by any preferred crystal orientation, implying its origin is mainly of enthalpic rather than entropic nature. Mean-field theory of polymer blends predicts that for immiscible polymers the melting point of the crystallizable component increases upon dilution in the other component, while it normally decreases for miscible blends. A local dilution is forced to occur at the diffuse interface of immiscible polymers; therefore the melting point of crystallizable polymers rises, which, in turn, enhances the thermodynamic driving force for crystal nucleation near the interface.

Understanding crystal orientation in quasi-one-dimensional polymer systems

We performed molecular simulations of polymer crystallization confined in rigid nanotubes. We found that the neutrally repulsive walls of nanotubes enhance the crystal nucleation. This enhancement gives uniform crystal orientation parallel to the tube axis at high temperatures. However, the nanotube walls that are wetted by the polymer, i.e. 'sticky' walls, switch the dominant crystal orientation from parallel to perpendicular. Since only perpendicularly oriented crystals can grow along the tube to give high crystallinities, the sticky walls provide a promising route towards the mass production of uniformly oriented nanocrystals in quasi-one-dimensional polymer systems.

Forces between thiolate-modified gold surfaces in a melt of end-functionalized polymers

To better understand surface forces across polymer melts, we measured the force between two chemically well-defined solid surfaces in a melt of polymer chains with a functional end group. As for surfaces, we used self-assembed monolayers (SAMs) of alkyl thiols with different end groups (methyl, amino, and hydroxyl) on gold. The polymer was a hydroxyl-terminated polyisoprene. To measure the force, an atomic force microscope was used. Between methyl-terminated SAMs, a weak and short-range repulsion was detected. Between hydroxyl or amino-terminated SAMs, a strong and long-range repulsion was observed up to distances of 16 nm. This indicates that the hydroxyl group of the polymer binds to the hydroxyl or amino groups of the SAMs. It forms a brush-like structure, which leads to steric repulsion. On amino-terminated SAMs, force-versus-distance curves on approach and retraction were monotonically repulsive and reversible. With hydroxyl-terminated SAMs, a jump was observed on approach when the load exceeded a certain threshold. On retraction, an adhesion had to be overcome. We interpret the jump as a rupture of the polymer layer. It indicates that the kinetics of bond and brush formation is faster on OH-SAMs than on NH2-SAMs.

Monte Carlo simulations of conformations of chain molecules in a cylindrical pore

Off-lattice Monte Carlo simulations are employed to study the behavior of chain molecules confined in a long cylindrical pore under the condition of hard-body interaction. The emphasis is placed on the chain and bead distributions as well as the location dependence of the chain conformation and anisotropy. The simulation results show that the chains very near the pore surface tend to wrap around the surface in various configurations. This behavior is qualitatively similar to that of the chains near but outside a cylindrical rod. Moreover, the bead concentration near the pore surface increases with increasing surface curvature.

Atomistic simulations to compute surface properties of poly(N-vinyl-2-pyrrolidone) (PVP) and blends of PVP/chitosan

Atomistic simulations were performed on poly(N-vinyl-2-pyrrolidone) (PVP) and its blends with chitosan (CS) in different ratios using molecular mechanics (MM) and molecular dynamics (MD) simulations in three-dimensionally periodic and effective two-dimensionally periodic condensed phases. Four independent microstructures were generated to analyze their surface properties. The calculated surface-energy values for PVP compared quite well with the experimental data reported in the literature. The density profile was analyzed, and the structure of the films showed an interior region of the bulk density. Various components of the energetic interactions (torsional, van der Waals, etc.) were examined to gain deeper insight into the nature of regular and anomalous interactions between the bulk and the surface films. Surface energies of PVP/CS blends were computed by MD simulations using the bulk pressure-volume-temperature (PVT) parameters. Bulk properties such as the cohesive energy density (CED) and solubility parameter (delta) were calculated using MM and MD simulations in the NVT ensemble under periodic boundary conditions. The Flory equation of state was used to compute the thermal expansion coefficient as well as PVT parameters. These surface-energy values agreed well with the surface-energy data calculated using the Zisman equation, which were also in accordance with the experimental observations. The results from this study suggest that computer simulations would provide valuable information on polymers and polymer-blend surfaces.

Surface forces in a confined polymer melt: self-consistent field analysis of full and restricted equilibrium cases

In full equilibrium the self-consistent field theory for a homopolymer melt confined between two surfaces predicts pronounced oscillatory interaction forces on the monomer length scale. However, when not all the polymer molecules can reversibly equilibrate with the bulk, the trapped molecules may be squeezed, adding a repulsive contribution to the interaction energy. The classical constrained or restricted equilibrium approach by Scheutjens and Fleer two decades ago to deal with this for polymers adsorbed from dilute solutions, breaks down in semidilute and concentrated polymer solutions. We present a generalized restricted equilibrium ansatz applicable also for concentrated polymer solutions. The key idea is that only the adsorbed polymer molecules, i.e., molecules that touch the surface at least once, are forced to remain inside the gap, whereas the nonadsorbing chains are free to move out of the gap when the surfaces approach each other. As in dilute solutions, the forces found in confined melt with trapped adsorbed chains become repulsive. We analyse the dependence of the interaction forces both in full as well as in restricted equilibrium cases as a function of the chain length and the interactions with the surface for a compressible polymer melt.

Confinement free energy and chain conformations of homopolymers confined between two repulsive walls

Lattice Monte Carlo simulations of polymer solutions confined between two parallel plates were performed. The confinement free energy Deltamicro(conf) per chain and the radius of gyrations of the chains parallel and perpendicular to the plates were obtained. When the concentration of the confined solution is above the overlap concentration, Deltamicro(conf) is found to scale with Na/D in a power law, betaDeltamicro(conf) approximately (Na/D)(m), with an exponent m=1.10+/-0.02 for athermal walls where N is the number of monomers in a chain, D is the slit width, and a is the lattice spacing. The presence of a weak attractive polymer/wall interaction epsilon(w) does not change the scaling variable, but the exponent m increases slightly. Extrapolating the results to melt would suggest that the predictions made by de Gennes [C. R. Acad. Sci. Paris II 305, 1181 (1987)] about the confinement free energy cost per chain in polymer melt is correct as far as the scaling variable is concerned, but is incorrect about the exponent m observed. The implication of this result on the predicted force between plates immersed in polymer melt is discussed. The parallel dimensions of the confined chain is expanded when the slit width D is narrow, however, the expansion is reduced at high concentration. It is conceivable that in melt the chain is not expanded when confined in a repulsive slit.

Thermodynamics and partitioning of homopolymers into a slit-A grand canonical Monte Carlo simulation study

Grand canonical ensemble Monte Carlo simulation (GCMC) combined with the histogram reweighting technique was used to study the thermodynamic equilibrium of a homopolymer solution between a bulk and a slit pore. GCMC gives the partition coefficients that agree with those from canonical ensemble Monte Carlo simulations in a twin box, and it also gives results that are not accessible through the regular canonical ensemble simulation such as the osmotic pressure of the solution. In a bulk polymer solution, the calculated osmotic pressure agrees very well with the scaling theory predictions both for the athermal polymer solution and the theta solution. However, one cannot obtain the osmotic pressure of the confined solution in the same way since the osmotic pressure of the confined solution is anisotropic. The chemical potentials in GCMC simulations were found to differ by a translational term from the chemical potentials obtained from canonical ensemble Monte Carlo simulations with the chain insertion method. This confirms the equilibrium condition of a polymer solution partition between the bulk and a slit pore: the chemical potentials of the polymer chain including the translational term are equal at equilibrium. The histogram reweighting method enables us to obtain the partition coefficients in the whole range of concentrations based on a limited set of simulations. Those predicted bulk-pore partition coefficient data enable us to perform further theoretical analysis. Scaling predictions of the partition coefficient at different regimes were given and were confirmed by the simulation data.

Properties of branched confined polymers

A model of star-branched polymer chains confined in a slit formed by two parallel surfaces was studied. The chains were embedded to a simple cubic lattice and consisted of f=3 branches of equal length. The macromolecules had the excluded volume and the confining surfaces were impenetrable for polymer segments. No attractive interactions between polymer segments and then between polymer segments and the surfaces were assumed and therefore the system was a thermal. Monte Carlo simulations were carried out employing the sampling algorithm based on chain's local changes of conformation. Lateral diffusion of star-branched chains was studied. Dynamic properties of star-branched chains between the walls with impenetrable rod-like obstacles were also studied and compared to the previous case. The density profiles of polymer segments on the slit were determined. The analysis of contacts between the polymer chain and the surfaces was also carried out.

Structure in Thin and Ultrathin Spin-Cast Polymer Films

The molecular organization in ultrathin polymer films (thicknesses less than 1000 angstroms) and thin polymer films (thicknesses between 1000 and 10,000 angstroms) may differ substantially from that of bulk polymers, which can lead to important differences in resulting thermophysical properties. Such constrained geometry films have been fabricated from amorphous poly(3-methyl-4-hydroxy styrene) (PMHS) and semicrystalline poly(di-n-hexyl silane) (PD6S) by means of spin-casting. The residual solvent content is substantially greater in ultrathin PMHS films, which suggests a higher glass transition temperature that results from a stronger hydrogen-bonded network as compared with that in thicker films. Crystallization of PD6S is substantially hindered in ultrathin films, in which a critical thickness of 150 angstroms is needed for crystalline morphology to exist and in which the rate of crystallization is initially slow but increases rapidly as the film approaches 500 angstroms in thickness.