Wen-Sheng Xu

Influence of Cohesive Energy and Chain Stiffness on Polymer Glass Formation

The generalized entropy theory is applied to assess the joint influence of the microscopic cohesive energy and chain stiffness on glass formation in polymer melts using a minimal model containing a single bending energy and a single (monomer averaged) nearest neighbor van der Waals energy. The analysis focuses on the combined impact of the microscopic cohesive energy and chain stiffness on the magnitudes of the isobaric fragility parameter mP and the glass transition temperature Tg. The computations imply that polymers with rigid structures and weak nearest neighbor interactions are the most fragile, while Tg becomes larger when the chains are stiffer and/or nearest neighbor interactions are stronger. Two simple fitting formulas summarize the computations describing the dependence of mP and Tg on the microscopic cohesive and bending energies. The consideration of the combined influence of the microscopic cohesive and bending energies leads to the identification of some important design concepts, such as i...

Lattice cluster theory for polymer melts with specific interactions

Despite the long-recognized fact that chemical structure and specific interactions greatly influence the thermodynamic properties of polymer systems, a predictive molecular theory that enables systematically addressing the role of chemical structure and specific interactions has been slow to develop even for polymer melts. While the lattice cluster theory (LCT) provides a powerful vehicle for understanding the influence of various molecular factors, such as monomer structure, on the thermodynamic properties of polymer melts and blends, the application of the LCT has heretofore been limited to the use of the simplest polymer model in which all united atom groups within the monomers of a species interact with a common monomer averaged van der Waals energy. Thus, the description of a compressible polymer melt involves a single van der Waals energy. As a first step towards developing more realistic descriptions to aid in the analysis of experimental data and the design of new materials, the LCT is extended here to treat models of polymer melts in which the backbone and side groups have different interaction strengths, so three energy parameters are present, namely, backbone-backbone, side group-side group, and backbone-side group interaction energies. Because of the great algebraic complexity of this extension, we retain maximal simplicity within this class of models by further specializing this initial study to models of polymer melts comprising chains with poly(n-α-olefin) structures where only the end segments on the side chains may have different, specific van der Waals interaction energies with the other united atom groups. An analytical expression for the LCT Helmholtz free energy is derived for the new model. Illustrative calculations are presented to demonstrate the degree to which the thermodynamic properties of polymer melts can be controlled by specific interactions.

Generalized Entropy Theory of Glass Formation in Polymer Melts with Specific Interactions

Chemical structure has been long recognized to greatly influence polymer glass formation, but a general molecular theory that predicts how chemical structure determines the properties of glass-forming polymers has been slow to develop. While the generalized entropy theory (GET) explains the influence of various molecular details on polymer glass formation, the application of the GET has heretofore been limited to the use of the simplest polymer model in which all united atom groups within the monomers of a species interact with a common monomer averaged van der Waals energy. However, energetic heterogeneities are ubiquitous within the monomers of real polymers, and their implications for polymer glass formation remain to be investigated theoretically. This paper uses an extension of the GET to explore the influence of energetic heterogeneities within monomers upon the nature of polymer glass formation. This extension of the GET is achieved by combining the Adam–Gibbs theory relating the structural relaxat...

Thermodynamic scaling of dynamics in polymer melts: predictions from the generalized entropy theory.

Many glass-forming fluids exhibit a remarkable thermodynamic scaling in which dynamic properties, such as the viscosity, the relaxation time, and the diffusion constant, can be described under different thermodynamic conditions in terms of a unique scaling function of the ratio ρ(γ)∕T, where ρ is the density, T is the temperature, and γ is a material dependent constant. Interest in the scaling is also heightened because the exponent γ enters prominently into considerations of the relative contributions to the dynamics from pressure effects (e.g., activation barriers) vs. volume effects (e.g., free volume). Although this scaling is clearly of great practical use, a molecular understanding of the scaling remains elusive. Providing this molecular understanding would greatly enhance the utility of the empirically observed scaling in assisting the rational design of materials by describing how controllable molecular factors, such as monomer structures, interactions, flexibility, etc., influence the scaling exponent γ and, hence, the dynamics. Given the successes of the generalized entropy theory in elucidating the influence of molecular details on the universal properties of glass-forming polymers, this theory is extended here to investigate the thermodynamic scaling in polymer melts. The predictions of theory are in accord with the appearance of thermodynamic scaling for pressures not in excess of ~50 MPa. (The failure at higher pressures arises due to inherent limitations of a lattice model.) In line with arguments relating the magnitude of γ to the steepness of the repulsive part of the intermolecular potential, the abrupt, square-well nature of the lattice model interactions lead, as expected, to much larger values of the scaling exponent. Nevertheless, the theory is employed to study how individual molecular parameters affect the scaling exponent in order to extract a molecular understanding of the information content contained in the exponent. The chain rigidity, cohesive energy, chain length, and the side group length are all found to significantly affect the magnitude of the scaling exponent, and the computed trends agree well with available experiments. The variations of γ with these molecular parameters are explained by establishing a correlation between the computed molecular dependence of the scaling exponent and the fragility. Thus, the efficiency of packing the polymers is established as the universal physical mechanism determining both the fragility and the scaling exponent γ.

Generalized entropy theory of glass-formation in fully flexible polymer melts

The generalized entropy theory (GET) offers many insights into how molecular parameters influence polymer glass-formation. Given the fact that chain rigidity often plays a critical role in understanding the glass-formation of polymer materials, the GET was originally developed based on models of semiflexible chains. Consequently, all previous calculations within the GET considered polymers with some degree of chain rigidity. Motivated by unexpected results from computer simulations of fully flexible polymer melts concerning the dependence of thermodynamic and dynamic properties on the cohesive interaction strength (ϵ), the present paper employs the GET to explore the influence of ϵ on glass-formation in models of polymer melts with a vanishing bending rigidity, i.e., fully flexible polymer melts. In accord with simulations, the GET for fully flexible polymer melts predicts that basic dimensionless thermodynamic properties (such as the reduced thermal expansion coefficient and isothermal compressibility) are universal functions of the temperature scaled by ϵ in the regime of low pressures. Similar scaling behavior is also found for the configurational entropy density in the GET for fully flexible polymer melts. Moreover, we find that the characteristic temperatures of glass-formation increase linearly with ϵ and that the fragility is independent of ϵ in fully flexible polymer melts, predictions that are again consistent with simulations of glass-forming polymer melts composed of fully flexible chains. Beyond an explanation of these general trends observed in simulations, the GET for fully flexible polymer melts predicts the presence of a positive residual configurational entropy at low temperatures, indicating a return to Arrhenius relaxation in the low temperature glassy state.

Spatial Rearrangement and Mobility Heterogeneity of an Anionic Lipid Monolayer Induced by the Anchoring of Cationic Semiflexible Polymer Chains

We use Monte Carlo simulations to investigate the interactions between cationic semiflexible polymer chains and a model fluid lipid monolayer composed of charge-neutral phosphatidyl-choline (PC), tetravalent anionic phosphatidylinositol 4,5-bisphosphate (PIP2), and univalent anionic phosphatidylserine (PS) lipids. In particular, we explore how chain rigidity and polymer concentration influence the spatial rearrangement and mobility heterogeneity of the monolayer under the conditions where the cationic polymers anchor on the monolayer. We find that the anchored cationic polymers only sequester the tetravalent PIP2 lipids at low polymer concentrations, where the interaction strength between the polymers and the monolayer exhibits a non-monotonic dependence on the degree of chain rigidity. Specifically, maximal anchoring occurs at low polymer concentrations, when the polymer chains have an intermediate degree of rigidity, for which the PIP2 clustering becomes most enhanced and the mobility of the polymer/PIP2 complexes becomes most reduced. On the other hand, at sufficiently high polymer concentrations, the anchoring strength decreases monotonically as the chains stiffen—a result that arises from the pronounced competitions among polymer chains. In this case, the flexible polymers can confine all PIP2 lipids and further sequester the univalent PS lipids, whereas the stiffer polymers tend to partially dissociate from the monolayer and only sequester smaller PIP2 clusters with greater mobilities. We further illustrate that the mobility gradient of the single PIP2 lipids in the sequestered clusters is sensitively modulated by the cooperative effects between anchored segments of the polymers with different rigidities. Our work thus demonstrates that the rigidity and concentration of anchored polymers are both important parameters for tuning the regulation of anionic lipids.

Lattice model of linear telechelic polymer melts. I. Inclusion of chain semiflexibility in the lattice cluster theory

The lattice cluster theory (LCT) for the thermodynamics of polymer systems has recently been reformulated to treat strongly interacting self-assembling polymers composed of fully flexible linear telechelic chains [J. Dudowicz and K. F. Freed, J. Chem. Phys. 136, 064902 (2012)]. Here, we further extend the LCT for linear telechelic polymer melts to include a description of chain semiflexibility, which is treated by introducing a bending energy penalty whenever a pair of consecutive bonds from a single chain lies along orthogonal directions. An analytical expression for the Helmholtz free energy is derived for the model of semiflexible linear telechelic polymer melts. The extension provides a theoretical tool for investigating the influence of chain stiffness on the thermodynamics of self-assembling telechelic polymers, and for further exploring the influence of self-assembly on glass formation in such systems.

Effects of Concentration and Ionization Degree of Anchoring Cationic Polymers on the Lateral Heterogeneity of Anionic Lipid Monolayers

We employed coarse-grained Monte Carlo simulations to investigate a system composed of cationic polymers and a phosphatidyl-choline membrane monolayer, doped with univalent anionic phosphatidylserine (PS) and tetravalent anionic phosphatidylinositol 4,5-bisphosphate (PIP2) lipid molecules. For this system, we consider the conditions under which multiple cationic polymers can anchor onto the monolayer and explore how the concentration and ionization degree of the polymers affect the lateral rearrangement and fluidity of the negatively charged lipids. Our work shows that the anchoring cationic polymers predominantly bind the tetravalent anionic PIP2 lipids and drag the PIP2 clusters to migrate on the monolayer. The polymer/PIP2 binding is found to be drastically enhanced by increasing the polymer ionization fraction, which causes the PIP2 lipids to form into larger clusters and reduces the mobility of the polymer/PIP2 complexes. As expected, stronger competition effects between anchoring polymers occur at higher polymer concentrations, for which each anchoring polymer partially dissociates from the monolayer and hence sequesters a smaller PIP2 cluster. The desorbed segments of the anchored polymers exhibit a faster mobility on the membrane, whereas the PIP2 clusters are closely restrained by the limited adhering cationic segments of anchoring polymers. We further demonstrate that the PIP2 molecules display a hierarchical mobility in the PIP2 clusters, which is regulated by the synergistic effect between the cationic segments of the polymers. The PS lipids sequester in the vicinity of the polymer/PIP2 complexes if the tetravalent PIP2 lipids cannot sufficiently neutralize the cationic polymers. Finally, we illustrate that the increase in the ionic concentration of the solution weakens the lateral clustering and the mobility heterogeneity of the charged lipids. Our work thus provides a better understanding of the fundamental biophysical mechanism of the concentration gradients and the hierarchical mobility of the anionic lipids in the membrane caused by the cationic polymer anchoring on length and time scales that are generally inaccessible by atomistic models. It also offers insight into the development and design of novel biological applications on the basis of the modulation of signaling lipids.

Glass formation in a mixture of hard disks and hard ellipses

We present an event-driven molecular dynamics study of glass formation in two-dimensional binary mixtures composed of hard disks and hard ellipses, where both types of particles have the same area. We demonstrate that characteristic glass-formation behavior appears upon compression under appropriate conditions in such systems. In particular, while a rotational glass transition occurs only for the ellipses, both types of particles undergo a kinetic arrest in the translational degrees of freedom at a single density. The translational dynamics for the ellipses is found to be faster than that for the disks within the same system, indicating that shape anisotropy promotes the translational motion of particles. We further examine the influence of mixtures composition and aspect ratio on the glass formation. For the mixtures with an ellipse aspect ratio of k = 2, both translational and rotational glass transition densities decrease with increasing the disk concentration at a similar rate, and hence, the two glass transitions remain close to each other at all concentrations investigated. By elevating k, however, the rotational glass transition density diminishes at a faster rate than the translational one, leading to the formation of an orientational glass for the ellipses between the two transitions. Our simulations imply that mixtures of particles with different shapes emerge as a promising model for probing the role of particle shape in determining the properties of glass-forming liquids. Furthermore, our work illustrates the potential of using knowledge concerning the dependence of glass-formation properties on mixtures composition and particle shape to assist in the rational design of amorphous materials.

Lattice model of linear telechelic polymer melts. II. Influence of chain stiffness on basic thermodynamic properties

The lattice cluster theory (LCT) for semiflexible linear telechelic melts, developed in Paper I, is applied to examine the influence of chain stiffness on the average degree of self-assembly and the basic thermodynamic properties of linear telechelic polymer melts. Our calculations imply that chain stiffness promotes self-assembly of linear telechelic polymer melts that assemble on cooling when either polymer volume fraction ϕ or temperature T is high, but opposes self-assembly when both ϕ and T are sufficiently low. This allows us to identify a boundary line in the ϕ-T plane that separates two regions of qualitatively different influence of chain stiffness on self-assembly. The enthalpy and entropy of self-assembly are usually treated as adjustable parameters in classical Flory-Huggins type theories for the equilibrium self-assembly of polymers, but they are demonstrated here to strongly depend on chain stiffness. Moreover, illustrative calculations for the dependence of the entropy density of linear telechelic polymer melts on chain stiffness demonstrate the importance of including semiflexibility within the LCT when exploring the nature of glass formation in models of linear telechelic polymer melts.