Orientation of Phenyl Rings and Methylene Bisectors at the Free Surface of Atactic Polystyrene

Predicting Polymer Brush Behavior in Solvents Using the Steepest-Entropy-Ascent Quantum Thermodynamic Framework

The steepest-entropy-ascent quantum thermodynamic (SEAQT) framework is utilized to study the effects of temperature on polymer brushes. The brushes are represented by a discrete energy spectrum, and energy degeneracies obtained through the replica-exchange Wang-Landau algorithm. The SEAQT equation of motion is applied to the density of states to establish a unique kinetic path from an initial thermodynamic state to a stable equilibrium state. The kinetic path describes the brush's evolution in state space, as it interacts with a thermal reservoir. The predicted occupation probabilities along the kinetic path are used to determine the expected thermodynamic and structural properties. The polymer density profile of a polystyrene brush in cyclohexane solvent is predicted using the equation of motion, and it agrees qualitatively with the experimental density profiles. The Flory-Huggins parameter chosen to describe brush-solvent interactions affects the solvent distribution in the brush but has a minimal impact on the polymer density profile. Three types of nonequilibrium kinetic paths with differing amounts of entropy production are considered: a heating path, a cooling path, and a heating-cooling path. Properties such as tortuosity, radius of gyration, brush density, solvent density, and brush chain conformations are calculated for each path.

Modeling of polystyrenic nanoparticles driven β-trans-crystalline efficiency in isotactic polypropylene

Addition of four styrenic soft nano-particles induced relative β-polymorphism as much as 20, 27, 34 and 10% during 10 minutes annealing of iPP at 116 °C. The kinetic and thermodynamic aspects of the phenomenon rationalized by a bi-exponential function

Molecular structure of poly(methyl methacrylate) surface. I. Combination of interface-sensitive infrared-visible sum frequency generation, molecular dynamics simulations, and ab initio calculations

The chemical composition and molecular structure of polymeric surfaces are important in understanding wetting, adhesion, and friction. Here, we combine interface-sensitive sum frequency generation spectroscopy (SFG), all-atom molecular dynamics (MD) simulations, and ab initio calculations to understand the composition and the orientation of chemical groups on poly(methyl methacrylate) (PMMA) surface as a function of tacticity and temperature. The SFG spectral features for isotactic and syndiotactic PMMA surfaces are similar, and the dominant peak in the spectra corresponds to the ester-methyl groups. The SFG spectra for solid and melt states are very similar for both syndiotactic and isotactic PMMA. In comparison, the MD simulation results show that both the ester-methyl and the α-methyl groups of syndiotactic-PMMA are ordered and tilted toward the surface normal. For the isotactic-PMMA, the α-methyl groups are less ordered compared to their ester-methyl groups. The backbone methylene groups have a broad angular distribution and are disordered, independent of tacticity and temperature. We have compared the SFG results with theoretical spectra calculated using MD simulations and ab initio calculations. Our analysis shows that the weaker intensity of α-methyl groups in SFG spectra is due to a combination of smaller molecular hyperpolarizability, lower ordering, and lower surface number density. This work highlights the importance of combining SFG spectroscopy with MD simulations and ab initio calculations in understanding polymer surfaces.

Interfacial properties of oxidized polystyrene and its interaction with water

All-atom molecular dynamics simulations have been carried out to study the wetting of atactic polystyrene (aPS) thin films by water droplets. The effect of oxidation of the aPS surface on the contact angle has been studied as a function of oxygen concentration. Oxidation of aPS has been achieved by randomly replacing with oxygen the ortho and/or meta hydrogens on the aromatic rings within 1 nm of the aPS surface until the desired concentration of oxygen is reached. The simulated contact angle is found to decrease monotonically with increasing degree of oxidation, consistent with recent experimental results. The number of hydrogen bonds between water molecules and polystyrene at the interface is found to monotonically increase with oxygen concentration. By use of a modified Good-Girafalco-Fowkes-Young equation, the contribution of nondispersion interactions, γsl(P), to the interfacial energy at the aPS/water interface has been determined as a function of the degree of oxidation. The values of γsl(P) extracted appear to follow a quadratic dependence on oxygen concentration of the aPS surface. The roughness of the polystyrene surface appears to be independent of oxygen concentration when the polystyrene is exposed to vacuum, and it appears to increase slightly when it is in contact with water. The orientational ordering of the phenyl rings at the polystyrene surface exhibits no dependence on oxygen concentration for polystyrene in vacuum. However, the ordering appears to decrease slightly with increasing oxygen concentration when the polystyrene is in contact with water.

Surface orientation of magainin 2: molecular dynamics simulation and sum frequency generation vibrational spectroscopic studies

We combined molecular dynamics based free energy calculations with sum frequency generation (SFG) spectroscopy to study the orientational distribution of solvated peptides near hydrophobic surfaces. Using a simplified atomistic model of the polystyrene (PS) surface, molecular dynamics simulations have been applied to compute the orientational probability of an α-helical peptide, magainin 2, with respect to the PS/water interface. Free energy calculations revealed that the preferred (horizontal) peptide orientation was driven by the favorable interactions between the hydrophobic PS surface and the hydrophobic residues on the helix, and additional simulations examined the importance of small aggregate formation. Concentration-dependent measurements obtained via SFG vibrational spectroscopy suggest that, at very low peptide concentrations, magainin molecules tend to lie down at the PS/solution interface, which correlates well with the simulation results. When the concentration is increased, peptides exhibit behavior not captured by MD simulations using single helical peptides. A combination of simulations and experiments was shown to yield more reliable results with molecular-level insights into interaction between peptides and polymer surfaces.

Adsorption of phenylphosphonic acid on GaAs (100) surfaces

The adsorption of phenylphosphonic acid (PPA) on GaAs (100) surfaces from solutions in acetonitrile/water mixtures was studied using Fourier transform infrared spectroscopy in attenuated total reflection in multiple internal reflections (ATR/MIR), X-ray photoelectron spectroscopy (XPS), high-resolution electron energy loss spectroscopy (HREELS), and atomic force microscopy (AFM). ATR/MIR in situ showed that the accumulation of PPA molecules near the GaAs surface increased with the water concentration in the solution. For water contents lower than 4%, ATR/MIR and XPS results are consistent with the formation of a low-density monolayer. A mechanism is proposed for H2O percentages lower than 4% involving the creation of interfacial bonds through a Brønsted acid-base reaction, which involves the surface hydroxyl groups most probably bound to Ga. It was found that the morphology of the final layer depended strongly on the water concentration in the adsorbing solution. For water concentrations equal to or higher than 5%, the amount of adsorbed molecules drastically increased and was accompanied by modifications in the infrared spectral region corresponding to P-O and P=O. This sudden change indicates a deprotonation of the acid. XPS studies revealed the presence of extra oxygen atoms as well as gallium species in the layer, leading to the conclusion that phosphonate and hydrogenophosphonate ions are present in the PPA layer intercalated with H3O+ and Ga3+ ions. This mechanism enables the formation of layers approximately 10 times thicker than those obtained with lower H2O percentages. HREELS indicated that the surface is composed of regions covered by PPA layers and uncovered regions, but the uncovered regions disappeared for water contents equal to or higher than 5%. XPS results are interpreted using a model consisting of a monolayer partially covering the surface and a thick layer. This model is consistent with AFM images revealing roughness on the order of 7 nm for the thick layer and 0.2-0.5 nm for the thin layer. Sonication proves to be an effective method for reducing layer thickness.

Mapping Atomistic Simulations to Mesoscopic Models: A Systematic Coarse-Graining Procedure for Vinyl Polymer Chains

The paper introduces a systematic procedure to coarse-grain atomistic models of the largest family of synthetic polymers into a mesoscopic model that is able to keep detailed information about chain stereosequences. The mesoscopic model consists of sequences of superatoms centered on methylene carbons of two different types according to the kind of diad (m or r) they belong to. The corresponding force-field contains three different bonds, six angle and three nonbonded terms. Recently developed analytical potentials, based on sums of Gaussians for bond and angle terms of the mesoscale force field have been used. For the nonbonded part, numerical potentials optimized by pressure-corrected iterative Boltzmann inversion have been used. As test case we coarse-grained an atomistic all-atom model of atactic polystyrene. The proposed mesoscale model has been successfully tested against structural and dynamical properties for different chain lengths and opens the possibility of relaxing melts of high molecular weight vinyl polymers.

Interior segment regrowth configurational-bias algorithm for the efficient sampling and fast relaxation of coarse-grained polyethylene and polyoxyethylene melts on a high coordination lattice

We demonstrate the application of a modified form of the configurational-bias algorithm for the simulation of chain molecules on the second-nearest-neighbor-diamond lattice. Using polyethylene and poly(ethylene-oxide) as model systems we show that the present configurational-bias algorithm can increase the speed of the equilibration by at least a factor of 2-3 or more as compared to the previous method of using a combination of single-bead and pivot moves along with the Metropolis sampling scheme [N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller, J. Chem. Phys. 21, 1087 (1953)]. The increase in the speed of the equilibration is found to be dependent on the interactions (i.e., the polymer being simulated) and the molecular weight of the chains. In addition, other factors not considered, such as the density, would also have a significant effect. The algorithm is an extension of the conventional configurational-bias method adapted to the regrowth of interior segments of chain molecules. Appropriate biasing probabilities for the trial moves as outlined by Jain and de Pablo for the configurational-bias scheme of chain ends, suitably modified for the interior segments, are utilized [T. S. Jain and J. J. de Pablo, in Simulation Methods for Polymers, edited by M. Kotelyanskii and D. N. Theodorou (Marcel Dekker, New York, 2004), pp. 223-255]. The biasing scheme satisfies the condition of detailed balance and produces efficient sampling with the correct equilibrium probability distribution of states. The method of interior regrowth overcomes the limitations of the original configurational-bias scheme and allows for the simulation of polymers of higher molecular weight linear chains and ring polymers which lack chain ends.

Effect of rubbing-induced polymer chain alignment on adhesion and friction of glassy polystyrene surfaces

The friction and adhesion properties of polystyrene surfaces are studied below the glass transition temperature by means of atomic force microscopy in argon. Even at a temperature far below the glass transition, the repeated sliding of a polystyrene bead tip on the non-cross-linked polystyrene surface causes significant reduction of friction and adhesion forces. There is no measurable wear of the polystyrene surface due to repeated sliding. These decreases are associated with the alignment of the outermost polymer segments induced by repeated rubbing. There are only little changes in friction and adhesion on the cross-linked polystyrene surface in which the covalent cross-linking prevents chain realignment.

Calculation of pressure using the virtual-volume-variation method and the virial method from chain conformations obtained by Monte Carlo simulations on the second nearest neighbor diamond lattice

For a model system of polyethylene of chain lengths 40 and 100 carbon atoms, we calculated the pressure at different densities and compared them with the experimental values. The simulation was conducted on the second nearest neighbor diamond lattice, and the pressure was calculated using the virtual-volume-variation method after the system was reverse mapped to its fully atomistic form in continuous space and energy minimized. In addition, the pressure was also calculated from the virial route by conducting a short molecular dynamics simulation starting from the energy minimized structure. We show that the pressure obtained from our simulations is quite reasonable in the length of simulation time (in Monte Carlo steps) normally employed in our group. These results provide additional evidence for the equilibration of our model systems, and methodology to calculate the pressure in our lattice models.