Lattice theory for binding of linear polymers to a solid substrate from polymer melts: I. Influence of chain connectivity on molecular binding and adsorption.

Most theories of the binding of molecules to surfaces or for the association between molecules treat the binding species as structureless entities and neglect their rigidity and the changes in their stiffness induced by the binding process. The binding species are also taken to be "ideal," meaning that the existence of van der Waals interactions and changes in these interactions upon molecular binding are also neglected. An understanding of the thermodynamics of these multifunctional molecular binding processes has recently come into focus in the context of the molecular binding of complex molecules, such as dendrimers and DNA grafted nanoparticles, to surfaces where the degree of binding cooperativity and selectivity, as well as the location of the binding transition, are found to be sensitive to the number of binding units constrained to a larger scale polymeric scaffold. We address the fundamental problem of molecular binding by extending classical Langmuir theory to describe the particular example of the reversible binding of semiflexible polymer chains to a solid substrate under melt conditions. The polymer chains are assumed to have a variable number N of binding units (segments) and to exhibit variable bending energies and van der Waals interactions in the bulk and on the surface, in addition to strong directional interactions with the surface. The resulting generalized Langmuir theory is applied to the examination of the influence of the chain connectivity of ideal polymers on the surface coverage Θ, transition binding temperature T1/2 at which Θ = 1/2, and on the derivative |dΘ/dT|T=T1/2 and the constant volume specific heat of binding, Cv bind, measures of the cooperativity and "sharpness" of the binding transition, respectively. Paper II is devoted to the impact of the van der Waals attractive interactions and chain stiffness on the reversible binding of nonideal polymer chains to a solid surface, including the enthalpy-entropy compensation phenomenon observed experimentally in many molecular and particle binding processes.

Lattice theory for binding of linear polymers to a solid substrate from polymer melts. II. Influence of van der Waals interactions and chain semiflexibility on molecular binding and adsorption.

The polymeric Langmuir theory, developed in Paper I [J. Dudowicz et al., J. Chem. Phys. 151, 124706 (2019)], is employed to investigate the influence of van der Waals interactions and chain rigidity on the thermodynamics of the reversible molecular binding to interfaces in one component polymer fluids (polymer melts). Both van der Waals interactions and chain stiffness are found to influence the temperature variation of the surface coverage Θ, the binding transition itself, and the cooperativity of molecular binding. Re-entrancy of the surface coverage Θ(T) is found to arise when the intermolecular interactions are sufficiently attractive to cause a liquid-vapor like phase separation in the interfacial region, a phenomenon that can occur in the binding of both small molecules and polymer chains to surfaces.

Lattice theory of competitive binding: Influence of van der Waals interactions on molecular binding and adsorption to a solid substrate from binary liquid mixtures.

The reversible binding of molecules to surfaces is one of the most fundamental processes in condensed fluids, with obvious applications in the molecular separation of materials, chromatographic characterization, and material processing. Motivated in particular by the ubiquitous occurrence of binding processes in molecular biology and self-assembly, we have developed a lattice type theory of competitive molecular binding to solid substrates from binary mixtures of two small molecule liquids that interact between themselves by van der Waals forces in addition to exhibiting binding interactions with the solid surface. The derived theory, in contrast to previously existing theoretical frameworks, enables us to investigate the influence of van der Waals interactions on interfacial binding and selective molecular adsorption. For reference, the classic Langmuir theory of adsorption is recovered when all van der Waals interaction energies between the molecules in the bulk liquid phase and those on the surface are formally set to zero. Illustrative calculations are performed for the binding of molecules to a solid surface from pure liquids and from their binary mixtures. The properties analyzed include the surface coverage θ, the binding transition temperature Tbind, the individual surface coverages, θA and θC, and the relative surface coverages, σAC≡θA/θC or σCA≡θC/θA. The latter two quantities coincide with the degrees of adsorption directly determined from experimental adsorption measurements. The Langmuir theory is shown to apply formally under a wide range of conditions where the original enthalpies (Δh or ΔhA and ΔhC) and entropies (Δs or ΔsA and ΔsC) of the binding reactions are simply replaced by their respective "effective" counterparts (Δheff or ΔhAeff and ΔhCeff and Δseff or ΔsAeff and ΔsCeff), whose values depend on the strength of der Waals interactions and of the "bare" free energy parameters (Δh or ΔhA and ΔhC, and Δs or ΔsA and ΔsC). Numerous instances of entropy-enthalpy compensation between these effective free energy parameters follow from our calculations, confirming previous reports on this phenomenon obtained from experimental studies of molecular binding processes in solution.

Mixtures of two self- and mutually-associating liquids: Phase behavior, second virial coefficients, and entropy-enthalpy compensation in the free energy of mixing.

The theoretical description of the phase behavior of polymers dissolved in binary mixtures of water and other miscible solvents is greatly complicated by the self- and mutual-association of the solvent molecules. As a first step in treating these complex and widely encountered solutions, we have developed an extension of Flory-Huggins theory to describe mixtures of two self- and mutually-associating fluids comprised of small molecules. Analytic expressions are derived here for basic thermodynamic properties of these fluid mixtures, including the spinodal phase boundaries, the second osmotic virial coefficients, and the enthalpy and entropy of mixing these associating solvents. Mixtures of this kind are found to exhibit characteristic closed loop phase boundaries and entropy-enthalpy compensation for the free energy of mixing in the low temperature regime where the liquid components are miscible. As discussed by Widom et al. [Phys. Chem. Chem. Phys. 5, 3085 (2003)], these basic miscibility trends, quite distinct from those observed in non-associating solvents, are defining phenomenological characteristics of the "hydrophobic effect." We find that our theory of mixtures of associating fluids captures at least some of the thermodynamic features of real aqueous mixtures.

Relation Between Solvent Quality and Phase Behavior of Ternary Mixtures of Polymers and Two Solvents that Exhibit Cononsolvency.

The phase boundaries of polymer solutions in mixed solvents can be extremely complex due to the many competing van der Waals and associative interactions that can arise in these ubiquitous and technologically important complex fluids. The present paper focuses specific attention on ternary solutions of polymers (B) dissolved in a mixture of two solvents (A, C) that competitively associate with the polymer. We are particularly concerned with explaining the origin of the peculiar phenomenon of cononsolvency in mixed solvents, where a mixture of two individually good solvents behaves effectively as a poor solvent. Our computations are based on a recently developed generalization of Flory-Huggins theory that incorporates the competitive solvation of a polymer by two associating solvents. On the basis of this framework, we evaluate the limit of polymer phase stability (spinodal curves) and the second osmotic virial coefficient [Formula: see text] as a function of temperature and the composition of the pure solvent mixture that is maintained in osmotic equilibrium with the ternary solution. The calculations provide new insights into the miscibility patterns of ternary A/B/C polymer solutions exhibiting cononsolvency.

Phase behavior and second osmotic virial coefficient for competitive polymer solvation in mixed solvent solutions.

We apply our recently developed generalized Flory-Huggins (FH) type theory for the competitive solvation of polymers by two mixed solvents to explain general trends in the variation of phase boundaries and solvent quality (quantified by the second osmotic virial coefficient B2) with solvent composition. The complexity of the theoretically predicted miscibility patterns for these ternary mixtures arises from the competitive association between the polymer and the solvents and from the interplay of these associative interactions with the weak van der Waals interactions between all components of the mixture. The main focus here lies in determining the influence of the free energy parameters for polymer-solvent association (solvation) and the effective FH interaction parameters {χαß} (driving phase separation) on the phase boundaries (specifically the spinodals), the second osmotic virial coefficient B2, and the relation between the positions of the spinodal curves and the theta temperatures at which B2 vanishes. Our classification of the predicted miscibility patterns is relevant to numerous applications of ternary polymer solutions in industrial formulations and the use of mixed solvent systems for polymer characterization, such as chromatographic separation where mixed solvents are commonly employed. A favorable comparison of B2 with experimental data for poly(methyl methacrylate)/acetonitrile/methanol (or 1-propanol) solutions only partially supports the validity of our theoretical predictions due to the lack of enough experimental data and the neglect of the self and mutual association of the solvents.

Communication: Cosolvency and cononsolvency explained in terms of a Flory-Huggins type theory.

Standard Flory-Huggins (FH) theory is utilized to describe the enigmatic cosolvency and cononsolvency phenomena for systems of polymers dissolved in mixed solvents. In particular, phase boundaries (specifically upper critical solution temperature spinodals) are calculated for solutions of homopolymers B in pure solvents and in binary mixtures of small molecule liquids A and C. The miscibility (or immiscibility) patterns for the ternary systems are classified in terms of the FH binary interaction parameters {χαß} and the ratio r = ϕ A /ϕ C of the concentrations ϕ A and ϕ C of the two solvents. The trends in miscibility are compared to those observed for blends of random copolymers (AxC1-x) with homopolymers (B) and to those deduced for A/B/C solutions of polymers B in liquid mixtures of small molecules A and C that associate into polymeric clusters {ApCq}i, (i = 1, 2, …, ∞). Although the classic FH theory is able to explain cosolvency and cononsolvency phenomena, the theory does not include a consideration of the mutual association of the solvent molecules and the competitive association between the solvent molecules and the polymer. These interactions can be incorporated in refinements of the FH theory, and the present paper provides a foundation for such extensions for modeling the rich thermodynamics of polymers in mixed solvents.

Theory of competitive solvation of polymers by two solvents and entropy-enthalpy compensation in the solvation free energy upon dilution with the second solvent.

We develop a statistical mechanical lattice theory for polymer solvation by a pair of relatively low molar mass solvents that compete for binding to the polymer backbone. A theory for the equilibrium mixture of solvated polymer clusters {AiBCj} and free unassociated molecules A, B, and C is formulated in the spirit of Flory-Huggins mean-field approximation. This theoretical framework enables us to derive expressions for the boundaries for phase stability (spinodals) and other basic properties of these polymer solutions: the internal energy U, entropy S, specific heat CV, extent of solvation Φsolv, average degree of solvation 〈Nsolv〉, and second osmotic virial coefficient B2 as functions of temperature and the composition of the mixture. Our theory predicts many new phenomena, but the current paper applies the theory to describe the entropy-enthalpy compensation in the free energy of polymer solvation, a phenomenon observed for many years without theoretical explanation and with significant relevance to liquid chromatography and other polymer separation methods.

The meaning of the "universal" WLF parameters of glass-forming polymer liquids.

Although the Williams-Landell-Ferry (WLF) equation for the segmental relaxation time τ(T) of glass-forming materials is one of the most commonly encountered relations in polymer physics, its molecular basis is not well understood. The WLF equation is often claimed to be equivalent to the Vogel-Fulcher-Tammann (VFT) equation, even though the WLF expression for τ(T) contains no explicit dependence on the fragility parameter D of the VFT equation, while the VFT equation lacks any explicit reference to the glass transition temperature Tg, the traditionally chosen reference temperature in the WLF equation. The observed approximate universality of the WLF parameters C1((g)) and C2((g)) implies that τ(T) depends only on T-Tg, a conclusion that seems difficult to reconcile with the VFT equation where the fragility parameter D largely governs the magnitude of τ(T). The current paper addresses these apparent inconsistencies by first evaluating the macroscopic WLF parameters C1((g)) and C2((g)) from the generalized entropy theory of glass-formation and then by determining the dependence of C1((g)) and C2((g)) on the microscopic molecular parameters (including the strength of the cohesive molecular interactions and the degree of chain stiffness) and on the molar mass of the polymer. Attention in these calculations is restricted to the temperature range (Tg < T < Tg + 100 K), where both the WLF and VFT equations apply.

Advances in the generalized entropy theory of glass-formation in polymer melts.

The generalized entropy theory (GET) of polymeric glass-forming liquids is reformulated into a computationally simpler and more natural formalism than the original version of this theory. The new theoretical framework greatly facilitates establishing essential trends in the dependence of the segmental relaxation time τ, fragility, characteristic temperatures of glass-formation, etc., on the combined influences of monomer molecular structure, chain rigidity, and cohesive interaction strength. Special attention is placed on the estimating the parameters of the phenomenological Vogel-Fulcher-Tammann relations for describing segmental relaxation in diverse liquids in the low temperature range of glass-formation, Tg > T > Tc (or Tg < T < Tg + 100 K), where Tg and Tc are, respectively, the glass transition temperature and the crossover temperature separating the high and low temperature regimes of glass-formation. Finally, we discuss how the molecular energetic interaction parameters of the GET can be estimated from experimental data. Illustrative calculations are performed for the stiffness factor σ and the cohesive energy density u as a first step in this direction.

Two glass transitions in miscible polymer blends?

In contrast to mixtures of two small molecule fluids, miscible binary polymer blends often exhibit two structural relaxation times and two glass transition temperatures. Qualitative explanations postulate phenomenological models of local concentration enhancements due to chain connectivity in ideal, fully miscible systems. We develop a quantitative theory that explains qualitative trends in the dynamics of real miscible polymer blends which are never ideal mixtures. The theory is a synthesis of the lattice cluster theory of blend thermodynamics, the generalized entropy theory for glass-formation in polymer materials, and the Kirkwood-Buff theory for concentration fluctuations in binary mixtures.

Concentration fluctuations in miscible polymer blends: influence of temperature and chain rigidity.

In contrast to binary mixtures of small molecule fluids, homogeneous polymer blends exhibit relatively large concentration fluctuations that can strongly affect the transport properties of these complex fluids over wide ranges of temperatures and compositions. The spatial scale and intensity of these compositional fluctuations are studied by applying Kirkwood-Buff theory to model blends of linear semiflexible polymer chains with upper critical solution temperatures. The requisite quantities for determining the Kirkwood-Buff integrals are generated from the lattice cluster theory for the thermodynamics of the blend and from the generalization of the random phase approximation to compressible polymer mixtures. We explore how the scale and intensity of composition fluctuations in binary blends vary with the reduced temperature τ ≡ (T - T(c))/T (where T(c) is the critical temperature) and with the asymmetry in the rigidities of the components. Knowledge of these variations is crucial for understanding the dynamics of materials fabricated from polymer blends, and evidence supporting these expectations is briefly discussed.

Solvation of polymers as mutual association. I. General theory.

A Flory-Huggins (FH) type lattice theory of self-assembly is generalized to describe the equilibrium solvation of long polymer chains B by small solvent molecules A. Solvation is modeled as a thermally reversible mutual association between the polymer and a relatively low molar mass solvent. The FH Helmholtz free energy F is derived for a mixture composed of the A and B species and the various possible mutual association complexes AiB, and F is then used to generate expressions for basic thermodynamic properties of solvated polymer solutions, including the size distribution of the solvated clusters, the fraction of solvent molecules contained in solvated states (an order parameter for solvation), the specific heat (which exhibits a maximum at the solvation transition), the second and the third osmotic virial coefficients, and the boundaries for phase stability of the mixture. Special attention is devoted to the analysis of the "entropic" contribution χ(s) to the FH interaction parameter χ of polymer solutions, both with and without associative interactions. The entropic χ(s) parameter arises from correlations associated with polymer chain connectivity and disparities in molecular structure between the components of the mixture. Our analysis provides the first explanation of the longstanding enigma of why χ(s) for polymer solutions significantly exceeds χ(s) for binary polymer blends. Our calculations also reveal that χ(s) becomes temperature dependent when interactions are strong, in sharp contrast to models currently being used for fitting thermodynamic data of associating polymer-solvent mixtures, where χ(s) is simply assumed to be an adjustable constant based on experience with solutions of homopolymers in nonassociating solvents.

Solvation of polymers as mutual association. II. Basic thermodynamic properties.

The theory of equilibrium solvation of polymers B by a relatively low molar mass solvent A, developed in the simplest form in Paper I, is used to explore some essential trends in basic thermodynamic properties of solvated polymer solutions, such as the equilibrium concentrations of solvated polymers AiB and free solvent molecules A, the mass distribution φ(AiB)(i) of solvated clusters, the extent of solvation of the polymer Φ(solv), the solvation transition lines T(solv)(φB(o)), the specific heat C(V), the osmotic second virial coefficient B2, phase stability boundaries, and the critical temperatures associated with closed loop phase diagrams. We discuss the differences between the basic thermodynamic properties of solvated polymers and those derived previously for hierarchical mutual association processes involving the association of two different species A and B into AB complexes and the subsequent polymerization of these AB complexes into linear polymeric structures. The properties of solvated polymer solutions are also compared to those for solutions of polymers in a self-associating solvent. Closed loop phase diagrams for solvated polymer solutions arise in the theory from the competition between the associative and van der Waals interactions, a behavior also typical for dispersed molecular and nanoparticle species that strongly associate with the host fluid. Our analysis of the temperature dependence of the second osmotic virial coefficient reveals that the theory must be generalized to describe the association of multiple solvent molecules with each chain monomer, and this complex extension of the present model will be developed in subsequent papers aimed at a quantitative rather than qualitative treatment of solvated polymer solutions.

Lattice cluster theory of associating telechelic polymers. III. Order parameter and average degree of self-assembly, transition temperature, and specific heat.

The lattice cluster theory of strongly interacting, structured polymer fluids is applied to determine the thermodynamic properties of solutions of telechelic polymers that may associate through bifunctional end groups. Hence, this model represents a significant albeit natural extension of a diverse array of prior popular equilibrium polymerization models in which structureless "bead" monomers associate into chain-like clusters under equilibrium conditions. In particular, the thermodynamic description of the self-assembly of linear telechelic chains in small molecule solvents (initiated in Paper II) is systematically extended through calculations of the order parameter Φ and average degree of self-assembly, the self-assembly transition temperature T(p), and the specific heat C(V) of solutions of telechelic molecules. Special focus is placed on examining how molecular and thermodynamic parameters, such as the solution composition φ, temperature T, microscopic interaction energies (ε(s) and ε), and length M of individual telechelic chains, influence the computed thermodynamic quantities that are commonly used to characterize self-assembling systems.

Lattice cluster theory of associating polymers. IV. Phase behavior of telechelic polymer solutions.

The newly developed lattice cluster theory (in Paper I) for the thermodynamics of solutions of telechelic polymers is used to examine the phase behavior of these complex fluids when effective polymer-solvent interactions are unfavorable. The telechelics are modeled as linear, fully flexible, polymer chains with mono-functional stickers at the two chain ends, and these chains are assumed to self-assemble upon cooling. Phase separation is generated through the interplay of self-assembly and polymer/solvent interactions that leads to an upper critical solution temperature phase separation. The variations of the boundaries for phase stability and the critical temperature and composition are analyzed in detail as functions of the number M of united atom groups in a telechelic chain and the microscopic nearest neighbor interaction energy ε(s) driving the self-assembly. The coupling between self-assembly and unfavorable polymer/solvent interactions produces a wide variety of nontrivial patterns of phase behavior, including an enhancement of miscibility accompanying the increase of the molar mass of the telechelics under certain circumstances. Special attention is devoted to understanding this unusual trend in miscibility.

Lattice cluster theory of associating polymers. I. Solutions of linear telechelic polymer chains.

The lattice cluster theory (LCT) for the thermodynamics of a wide array of polymer systems has been developed by using an analogy to Mayer's virial expansions for non-ideal gases. However, the high-temperature expansion inherent to the LCT has heretofore precluded its application to systems exhibiting strong, specific "sticky" interactions. The present paper describes a reformulation of the LCT necessary to treat systems with both weak and strong, "sticky" interactions. This initial study concerns solutions of linear telechelic chains (with stickers at the chain ends) as the self-assembling system. The main idea behind this extension of the LCT lies in the extraction of terms associated with the strong interactions from the cluster expansion. The generalized LCT for sticky systems reduces to the quasi-chemical theory of hydrogen bonding of Panyioutou and Sanchez when correlation corrections are neglected in the LCT. A diagrammatic representation is employed to facilitate the evaluation of the corrections to the zeroth-order approximation from short range correlations.

Lattice cluster theory of associating polymers. II. Enthalpy and entropy of self-assembly and Flory-Huggins interaction parameter χ for solutions of telechelic molecules.

The lattice cluster theory for solutions of telechelic polymer chains, developed in paper I, is applied to determine the enthalpy Δh(p) and entropy Δs(p) of self-assembly of linear telechelics and to evaluate the Flory-Huggins (FH) interaction parameter χ governing the phase behavior of these systems. Particular focus is placed on examining how these interaction variables depend on the composition of the solution, temperature, van der Waals and local "sticky" interaction energies, and the length of the individual telechelic chains. The FH interaction parameter χ is found to exhibit an entropy-enthalpy compensation effect between the "entropic" and "enthalpic" portions as either the composition or mass of the telechelic species is varied, providing unique theoretical insights into this commonly reported, yet, enigmatic phenomenon.

Can the Miscibility of Telechelic Polymer Solutions Increase with Polymer Chain Length?

Increasing the molar mass of the polymers in blends and in solutions tends to decrease miscibility, but application of the lattice cluster theory for strongly interactiong polymer systems to telechelic polymer solutions explains why this usual trend can be inverted, a situation actually observed in some telechelic polymer solutions and blends.

General approach to polymer chains confined by interacting boundaries.

Polymer chains, confined to cavities or polymer layers with dimensions less than the chain radius of gyration, appear in many phenomena, such as gel chromatography, rubber elasticity, viscolelasticity of high molar mass polymer melts, the translocation of polymers through nanopores and nanotubes, polymer adsorption, etc. Thus, the description of how the constraints alter polymer thermodynamic properties is a recurrent theoretical problem. A realistic treatment requires the incorporation of impenetrable interacting (attractive or repulsive) boundaries, a process that introduces significant mathematical complications. The standard approach involves developing the generalized diffusion equation description of the interaction of flexible polymers with impenetrable confining surfaces into a discrete eigenfunction expansion, where the solutions are normally truncated at the first mode (the "ground state dominance" approximation). This approximation is mathematically well justified under conditions of strong confinement, i.e., a confinement length scale much smaller than the chain radius of gyration, but becomes unreliable when the polymers are confined to dimensions comparable to their typically nanoscale size. We extend a general approach to describe polymers under conditions of weak to moderate confinement and apply this semianalytic method specifically to determine the thermodynamics and static structure factor for a flexible polymer confined between impenetrable interacting parallel plate boundaries. The method is first illustrated by analyzing chain partitioning between a pore and a large external reservoir, a model system with application to chromatography. Improved agreement is found for the partition coefficients of a polymer chain in the pore geometry. An expression is derived for the structure factor S(k) in a slit geometry to assist in more accurately estimating chain dimensions from scattering measurements for thin polymer films.