Next Article in Journal
Modification of Spherical Polyelectrolyte Brushes by Layer-by-Layer Self-Assembly as Observed by Small Angle X-ray Scattering
Next Article in Special Issue
Modelling and Validation of Synthesis of Poly Lactic Acid Using an Alternative Energy Source through a Continuous Reactive Extrusion Process
Previous Article in Journal
Analysis and Testing of Bisphenol A—Free Bio-Based Tannin Epoxy-Acrylic Adhesives
Previous Article in Special Issue
Decorating Nanoparticle Surface for Targeted Drug Delivery: Opportunities and Challenges
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Density Functional Theory of Polymer Structure and Conformations

1
State Key Laboratory of Organic-Inorganic Composites, Beijing University of Chemical Technology, Beijing 100029, China
2
Key Laboratory of Beijing City on Preparation and Processing of Novel Polymer Materials, Beijing University of Chemical Technology, Beijing 100029, China
*
Authors to whom correspondence should be addressed.
Polymers 2016, 8(4), 121; https://doi.org/10.3390/polym8040121
Submission received: 8 March 2016 / Revised: 26 March 2016 / Accepted: 30 March 2016 / Published: 15 April 2016
(This article belongs to the Special Issue Computational Modeling and Simulation in Polymer)

Abstract

:
We present a density functional approach to quantitatively evaluate the microscopic conformations of polymer chains with consideration of the effects of chain stiffness, polymer concentration, and short chain molecules. For polystyrene (PS), poly(ethylene oxide) (PEO), and poly(methyl methacrylate) (PMMA) melts with low-polymerization degree, as chain length increases, they display different stretching ratios and show non-universal scaling exponents due to their different chain stiffnesses. In good solvent, increase of PS concentration induces the decline of gyration radius. For PS blends containing short ( m 1 = 1 100 ) and long ( m = 100 ) chains, the expansion of long chains becomes unobvious once m 1 is larger than 40, which is also different to the scaling properties of ideal chain blends.

Graphical Abstract

1. Introduction

The description of macromolecular conformations in various environments is an outstanding problem in polymer physics [1,2,3,4,5,6,7,8]. The scaling properties may be altered in different media and polymer ratios. Theoretical concepts of polymer physics have frequently been used to establish the connection between polymer microstructure and dimension [9,10,11]. On the other hand, the advances in experimental techniques, such as neutron scattering methods, allow a direct measurement of the intramolecular and intermolecular structure, thus providing, in principle, the possibility of testing these theories [11,12,13,14].
The stretching of polymers has been treated exhaustively by a number of Flory-type mean field theories based on general principles that relate the dimensions and the length of a chain to its free-energy [15,16,17,18]. These theories provide reasonable scaling exponent that correlate the gyration radius ( R g ) of polymer chain to its length ( N ) via R g N v . A prerequisite for application of these theories is that R g at θ-state are known, which serves as an essential reference state. Although the scaling exponent at θ-state has been determined as ν = 1 / 2 , the pre-coefficient of the relational expression is still unavailable. Moreover, since the free-energy expression is independent of microscopic inter- and intramolecular structures, the quantitative application of these theories and their comparison to experimental results in different conditions are still far from satisfactory.
Polymer reference interaction site model (PRISM) integral equation provides a framework for describing the microscopic structures of uniform polymeric fluids [19,20]. The intermolecular correlation functions are calculated for a given set of intramolecular correlation functions after the equation is formulated and coupled with a closure relation. It has been successfully applied to describe the intermolecular structure of a broad range of polymeric systems [21,22,23]. Since the intramolecular correlation function is typically unknown, the equation suffers from a serious self-consistency problem. Although this problem can be partially solved by involving a semiflexible chain model or inserting a single-chain molecular simulation to integrate the intra- and intermolecular correlation functions, this strategy is restricted to improve the accuracy of intermolecular correlation function, whereas the intramolecular correlation function cannot be derived from the equation. Besides, the integral equation approach suffers inconsistency among different routes for thermodynamic properties. As such, the structure description provided by the PRISM is generally qualitative. In order to present a quantitative evaluation to the intramolecular correlation function, an explicit free-energy analysis of the system is undoubtedly necessary.
Classical density functional theory (DFT) provides another route to represent both structural and thermodynamic properties of polymeric fluids [24,25,26]. The theory for fluids is based on the minimization of grand free-energy functional, and gives prediction for the equilibrium free-energy and microscopic structure of the components. In particular, the theory yields self-contained structural and thermodynamic properties and employs no molecular simulations as input. In recent years, an accurate density functional approach has been developed for polymer systems by combining the modified interfacial statistical association fluid theory [27,28,29,30]. The excess Helmholtz free-energy due to polymerization is related to the association equilibrium through the multibody cavity correlation functions. In order to characterize the intramolecular correlation function, a test-particle method has been integrated into a relatively simple DFT to calculate the local inhomogeneous density profile of the polymeric fluid in the external field of one segment fixed at the origin [31,32]. Such DFT is based on the first order thermodynamic perturbation theory for polyatomic molecules. It avoids molecular simulations as input and shows the advantages of self-consistency among inter- and intramolecular correlation functions, and is very accurate in comparison with simulation data for freely jointed hard-sphere chains. The drawback of the theory lies in its ideal chain model. In fact, different polymers in experiments have different degrees of stiffness to their backbones.
There is considerable interest in developing a theory capable of accurately predicting the microscopic conformation of complicated polymer fluids to decipher the stretching properties. In this paper, we present a density functional approach by combing the PRISM equation, the test-particle method, and the modified interfacial statistical association fluid theory [32] to deal with actual polymer conformations, particularly the intramolecular correlation functions. The radii of gyration of polystyrene (PS), poly(ethylene oxide) (PEO), and poly(methyl methacrylate) (PMMA) are calculated to analyze their expansion or contraction under different conditions. We utilize coarse-graining chain for modeling the conformational behavior of the polymers with chain stiffness, reflecting the dependence of bending energy on the angle between two nearby bond or tangent vectors. Since all parameters are taken from the general force field [33,34,35], all results given by the theory are strict predictions.
The rest of the paper is organized as follows: Section 2 presents the whole molecular model and the DFT calculation method. Section 3 shows the calculation results of gyration radius and structure factor. The expansion effects are evaluated for polymer melts, concentration, and blends containing long and short chains. Section 4 summarizes the final conclusion.

2. Model and Approach

To simplify the complicated polymer conformation, we use a coarse-grained model to describe polymer chains. Each polymer chain is represented by a series of bonded segments and obeys the semiflexible chain conformation. The segment diameter is given by repeating unit. Suppose that one segment from an arbitrarily selected chain is fixed at the origin. The system considered is equivalent to a mixture of four polymeric components ( F ,   S ,   C   and   D ) in a symmetric external field due to the fixed segment, as depicted in Figure 1.
The free molecules are represented by F , composed of m F segments; the chemically identical short chain molecules are represent by S , composed of m S segments; while the tethered fragments are represented by C and D , composed of m C and m D segments, respectively. At equilibrium, the density distributions of free chains, short chains, and the tethered fragments satisfy the variational relations
δ Ω δ ρ i ( F ) ( r ) = δ Ω δ ρ i ( S ) ( r ) = δ Ω δ ρ i ( C ) ( r ) = δ Ω δ ρ i ( D ) ( r ) = 0
where Ω stands for the grand free-energy, ρ i ( l ) ( r ) ( l = F ,   S ,   C ,   and   D ) is the density of the i th segment on chain l at position r . For a homogeneous polymer melt, the contribution of the S component is automatically deleted.
The segment distributions of the free ( F ) and short ( S ) chains around the fixed segment are related to the intermolecular site–site correlation function g i j ( r ) = ρ i , j ( r ) / ρ b , where ρ i , j ( r ) is the density profile of segments i on molecule F or S around the fixed segment j , and ρ b is the bulk density of segment. Whereas the distribution of segments from the fragments C and D are related to the intramolecular correlation function ω i j ( r ) = ρ i , j ( r ) , where ρ i , j ( r ) is the density profile of segment i on the tethered chain ( C or D ) around the fixed segment j . There is only one tethered polymer chain, therefore the normalization condition is ω i , j ( r ) d r = 1 . As a consequence, the site–site intra- and intermolecular correlation functions specify all detail microscopic structures of a polymeric fluid. In particular, the gyration radius is a version of the intramolecular correlation function. From the site–site correlation functions, we can calculate the average intermolecular correlation function with g ( r ) = 1 m F 2 i = 1 m F j = 1 m F g i j ( r ) , and the average intramolecular correlation function with ω ( r ) = 1 m F i = 1 m F j = 1 m F ω i j ( r ) .
The system considered above is equivalent to a mixture of four polymeric components ( F + S + C + D ) in a symmetric external field due to the fixed segment. The grand potential is related to the Helmholtz free-energy A [ ρ ( r ) ] through the Legendre transform
Ω [ ρ i ( F ) ( r ) , ρ i ( S ) ( r ) , ρ i ( C ) ( r ) , ρ i ( D ) ( r ) ] = A [ ρ i ( F ) ( r ) , ρ i ( S ) ( r ) , ρ i ( C ) ( r ) , ρ i ( D ) ( r ) ] l = F , S , C , D i = 1 m l d r ( μ i ( l ) V ext i , l ( r ) ) ρ i ( l ) ( r )
where μ i ( l ) is the chemical potential of that segment i , V ext i , l is the external field acting on that segment i . The first summation is over all chains l in the mixture ( F , S , C , D ) , and the second summation is over all segments of chain l . For any segment that is not immediately bonded with the fixed segment, the external potential is identical to the Lennard–Jones (LJ) interaction potential. While for two segments that are directly connected to the fixed segment, the external potential includes the bonding potential v bond , given by exp [ β v bond ( r i , r j ) ] = δ ( | r i r j | σ ) / 4 π σ 2 ; here segments i and j are nearest neighbors from the same molecule, and δ is the Dirac-delta function. The bond length ( L ), bond angle ( θ ), and nonbonded parameters ( σ and ε ) for PS, PEO, and PMMA are summarized in Table 1.
The total Helmholtz free-energy functional can be decomposed into ideal and excess contributions [27]. The ideal contribution is generally represented by A id = d r ρ i ( l ) ( r ) [ ln [ ρ i ( l ) ( r ) ] 1 ] . The excess contribution consists of hard-sphere repulsion, long-range attraction, chain stiffness, and chain connectivity over the ideal gas state of the atomic mixture.
The free-energy due to hard-sphere repulsion is given by the fundamental measure theory [24]
A hs [ ρ i ( l ) ( r ) ] = d r Φ hs [ n γ ( r ) ]
where Φ hs [ n γ ( r ) ] is the free-energy density, which stems from the modified fundamental theory [36], including both the scalar and vector contributions
Φ hs [ n γ ( r ) ] = [ n 0 ln ( 1 n 3 ) + n 1 n 2 n V 1 n V 2 1 n 3 + 1 36 π ( n 3 ln ( 1 n 3 ) + n 3 2 ( 1 n 3 ) 2 ) n 2 3 3 n 2 n V 2 n V 2 n 3 3 ]
where n γ ( r ) with γ = 0 ,     1 ,     2 ,     3 ,     V 1 ,     V 2 are the weighted densities. The details have been given elsewhere [36].
The long-range attractive contribution to the free-energy functional can be simplified as
A att [ ρ i ( l ) ( r ) ] = 1 2 i j d r d r ρ i ( l ) ( r ) ρ j ( l ) ( r ) u i j att ( | r - r | )
where u i j att ( r ) is the interaction potential between any two species i and j , and can be represented by a cut-and-shifted LJ potential with a Weeks–Chandler–Anderson separation [37,38].
The functional A stiff is constructed to account for the contribution of conformational entropy given by the polymer stiffness, which is implemented to improve the accuracy for description of real polymer chains
A stiff [ ρ i ( l ) ( r ) ] = 1 2 i j d r d r ρ i ( l ) ( r ) ρ j ( l ) ( r ) c i j stiff ( | r r | )  
with the approximation c stiff ( r ) = c semiflexible ( r ) c flexible ( r ) . Here c semiflexible ( r ) and c flexible ( r ) denote the direct correlation functions of semiflexible and flexible polymer chains. They are calculated from the PRISM integral equation [39]
h ( r ) = d r d r ω ( | r r | ) c ( | r r | ) [ ω ( r ) + ρ h ( r ) ]
where h ( r ) is the total correlation function and ω ( r ) is the intramolecular correlation function. ω ( r ) is cursorily represented by the Koyama model [40] for flexible or semiflexible chains. To solve the equation, we adopt the Kovalenko–Hirata approximation [41].
To compute the free-energy contribution due to chain connectivity, we use the Tripathi–Chapman functional [27], in which the chains are treated as a sequence of m bonded monomers, enforced by giving each segment a label and allowing segments to exclusively bond to their specific matching segments. It is indicated as
A chain [ ρ i ( l ) ( r ) ] = d r i = 1 m l ρ i ( l ) ( r ) × A Γ ( i ) ( ln X A i , l ( r ) X A i , l ( r ) 2 + 1 2 )
in which the first summation is over all the segments i , and the second is over all the association sites on segment i as Γ ( i ) , representing the set of all associating sites on segment i . X A i , l is the fraction of segment i that are not bonded at their association site A.
X A i . l ( r ) = 1 1 + d r X B j , l ( r ) Δ i j ( l ) ( r , r ) ρ j ( l ) ( r )
where j denotes the neighboring segment that will bond to segment i , and Δ i j ( l ) ( r , r ) = K F i j ( l ) ( r , r ) y i j ( l ) ( r , r ) . Here K is a geometric constant that accounts for the volume available for bonding between segments, and F i j ( l ) ( r , r ) = exp ( β ε 0 β v bond i j ( r , r ) ) 1 represents the association Mayer f-function. y i j ( l ) ( r , r ) is the cavity correlation function. The details have been given elsewhere [27].
The functional derivatives of the free-energies are required to obtain the equilibrium density profiles, which are given as Ω / ρ i ( l ) ( r ) = 0 . As a result, the Euler–Lagrange Equation is written as
ln ρ i ( l ) ( r ) + A Γ ( i ) ln X A i , l ( r ) 1 2 q = F , S , C , D γ = 1 m q γ { γ } ρ γ ( q ) ( r ) δ ln y γ γ ( r ) δ ρ i ( l ) ( r ) d r + δ β A hs δ ρ i ( l ) ( r ) + δ β A att δ ρ i ( l ) ( r ) + δ β A stiff δ ρ i ( l ) ( r ) = β ( μ i ( l ) V ext i , l ( r ) )
where { γ } is the set of all segments bonded to segment γ . This equation can be rewritten to give the density profile
ρ i ( l ) ( r ) = exp ( β μ M l ) exp [ D i ( l ) ( r ) β V ext i , l ( r ) ] I 1 , i ( l ) ( r ) I 2 , i ( l ) ( r )
with
D i ( l ) ( r ) = 1 2 q = F , S , C , D γ = 1 m q γ { γ } ρ γ ( q ) ( r ) δ ln y γ γ ( r ) δ ρ i ( l ) ( r ) d r δ β A hs δ ρ i ( l ) ( r ) δ β A att δ ρ i ( l ) ( r ) δ β A stiff δ ρ i ( l ) ( r )
In Equation (11), μ M l is the bulk chemical potential of chain l . The multiple integrals I 1 , i ( l ) ( r )   and   I 2 , i ( l ) ( r ) for the free and short chains ( l = F , S ) are solved in a recursive fashion and are given by
{ I 1 , 1 ( l ) ( r ) = 1 I 1 , i ( l ) ( r ) = I 1 , i 1 ( l ) ( r ) exp [ D i 1 ( l ) ( r ) β V ext i 1 , l ( r ) ] Δ i 1 , i ( l ) ( r , r ) ( r Θ ( σ | r r | ) r ) d r
{ I 2 , i ( l ) ( r ) = I 2 , i + 1 ( l ) ( r ) exp [ D i + 1 ( l ) ( r ) β V ext i + 1 , l ( r ) ] Δ i , i + 1 ( l ) ( r , r ) ( r Θ ( σ | r r | ) r ) d r I 2 , m F ( l ) ( r ) = 1
where Θ ( r ) is the Heaviside step function. When applied to the tethered chains ( l = C , D ) , I 1 , i ( l ) ( r )   and I 2 , i ( l ) ( r ) can be expressed as
{ I 1 , 1 ( l ) ( σ ) = 1 I 1 , 2 ( l ) ( r ) = exp [ D 1 ( l ) ( σ ) β V ext 1 , l ( σ ) ] Δ 1 , 2 ( l ) ( σ , r ) ( σ Θ ( σ | σ r | ) r ) I 1 , i ( l ) ( r ) = I 1 , i 1 ( l ) ( r ) exp [ D i 1 ( l ) ( r ) β V ext i 1 , l ( r ) ] Δ i 1 , i ( l ) ( r , r ) ( r Θ ( σ | r r | ) r ) d r
{ I 2 , 1 ( l ) ( σ ) = I 2 , 2 ( l ) ( r ) exp [ D 2 ( l ) ( r ) β V ext 2 , l ( r ) ] Δ 1 , 2 ( l ) ( σ , r ) ( r Θ ( σ | r σ | ) σ ) d r I 2 , i ( l ) ( r ) = I 2 , i + 1 ( l ) ( r ) exp [ D i + 1 ( l ) ( r ) β V ext i + 1 , l ( r ) ] Δ i , i + 1 ( l ) ( r , r ) ( r Θ ( σ | r r | ) r ) d r I 2 , m l ( l ) ( r ) = 1
The chemical potentials of the tethered fragments C and D can be determined using the normalization conditions 4 π r 2 ρ i ( C ) ( r ) d r = 1 and 4 π r 2 ρ i ( D ) ( r ) d r = 1 , where i = 1 , 2 , 3 , ... m l ,   for   l = C   or   D . For the first tethered segment (segment “1” of l = C   or   D ),
4 π r 2 { exp ( β μ M l ) exp [ D 1 ( l ) ( r ) β V ext 1 , l ( r ) ] I 1 , 1 ( l ) ( r ) I 2 , 1 ( l ) ( r ) } d r = 1
which yields exp ( β μ M l ) = 1 / [ exp ( D 1 ( l ) ( σ ) ) I 1 , 1 ( l ) ( σ ) I 2 , 1 ( l ) ( σ ) ] . Using other segments yields equivalent results for μ M l . Substituting this into Equation (11), we obtain
ρ 1 ( C ) ( r ) = ρ 1 ( D ) ( r ) = δ ( r σ ) 4 π σ 2
which matches the known condition for the tethered segment. Solving for the other segments ( i = 2 , 3 , ... m l ) gives
ρ i ( l ) ( r ) = 1 exp [ D 1 ( l ) ( σ ) ] I 1 , 1 ( l ) ( σ ) I 2 , 1 ( l ) ( σ ) exp [ D i ( l ) ( r ) β V ext i , l ( r ) ] I 1 , i ( l ) ( r ) I 2 , i ( l ) ( r )
In calculating inter- and intramolecular correlation functions, we fix the segments of a polymer chain one by one and the density distributions around the fixed segment are calculated with Equations (11)–(19). Because of symmetry, m F / 2 (if m F is even) or ( m F + 1 ) / 2 (if m F is odd) calculations are required for predicting the detailed local structures of polymer chains consisting of m F identical segments. The density profiles are solved using the Picard-type iterative method. In the theoretical calculations, the computation domain is divided into equally-spaced grid points along the dimension normal to the surface. The grid spacing of 0.02 σ is used in our calculation. At every iteration step, a new estimate to the density profiles is calculated with Equation (11), and then is mixed with the old one. The result acts as the new guess for the next iteration. This procedure is iterated until numerical self-consistency is achieved, in the sense that the difference of density profile between the preceding step and the present step is less than 5.0 × 10 4 . All the convolution results of the equations are directly evaluated in the real space. In order to save computational time, we focus on those low-polymerized systems, where the monomer number is not more than 120. Unless noted otherwise, all calculations are performed at a total packing fraction of 0.40, which corresponds to a typical dense melt value with a realistic dimensionless isothermal compressibility.

3. Results and Discussion

Figure 2 presents the inter- and intramolecular correlation functions for the PS chain given by the DFT, along with the simulation results that are attainable from the literature [42]. It is shown that the overall shape of theoretical curves and the positions of the main peaks can match the simulation results, indicating that the current theoretical model is suitable to quantitatively evaluate the microscopic structure of polymer melts. The deviations are probably due to the approximation of the model, in which the tacticity or torsional angle effects are overlooked. As shown in Figure 2a, on short length scales, the local structure of the PS chain is seen to have “liquid-like” tendency, displaying shells of first-nearest neighbors. Over longer distances, the curves exhibit a spatially slowly varying “correlation hole”, corresponding to the relative absence of neighboring sites due to intramolecular screening. Figure 2b presents the corresponding “nonbonded” intramolecular correlation functions. The discontinuity at r = 7 Å is due to the direct interaction between next nearest neighbors along the polymer chain.
As the essential criterion to evaluate polymer structure and dimension, the gyration radius R g can be calculated from the intramolecular correlation function [40,43]
R g = 2 π m 0 r 4 ω ( r ) d r
to quantitatively analyze polymer conformation. Figure 3a shows R g for PS chains at T = 500   K as function of chain length in comparison with the available experimental data [44]. It is clear that both the PRISM equation and the Flory theory cannot provide reasonable intramolecular structure description. In contrast, a good agreement has been achieved between DFT calculations and experimental values in the whole range of chain length. In particular, the scaling exponents derived from the Flory theory and the DFT are fitted in Figure 3b. The values are 0.50 and 0.63, respectively. It shows clearly that, at a low degree of polymerization, the exponent from DFT is obviously higher than the universal value (0.50). In other words, actual polymers at a low degree of polymerization display non-universal scaling behavior and have a relatively larger scaling exponent, which can be attributed to their rigid characteristics. Moreover, as the chain becomes stiffer, the excluded volume increases, leading to increasing R g . A similar result was also derived from molecular simulations [45,46].
Another consequence of the gradual ramp-up of polymer stiffness manifests itself through the form factor, which is particularly useful for studying the scaling behavior. Generally, the form factor is derived from light scattering experiments. In the theoretical model, it can be directly calculated from the Fourier transform of the intramolecular correlation function [47]
P ( q ) = 1 m 2 i = 1 m j = 1 m sin ( q r i j ) q r i j = 1 m 2 i = 1 m j = 1 m ω i j ( q ) = 1 m ω ( q )
where r i j is the distance between segments i and j , and q denotes the magnitude of the scattering wave vector q . The expression for the form factor of a free chain is the well-known Debye equation
P ( q ) = 2 ( q 2 R g 2 ) 2 ( exp ( q 2 R g 2 ) 1 + q 2 R g 2 )
Figure 4a shows the form factors for PS melt with molecular weight 1 × 104 in θ solvent. It seems that a good agreement can be achieved between DFT calculations and the experimental values [48], whereas the Debye expression is not always valid in the fractal regime ( q 1.0 ). This can be seen by the gradual transition from a slope of −2 for the Debye equation to a slope of −1 for PS chains. Figure 4b shows the form factor of several different lengths of PS melt. In particular, we are interested in the so-called fractal regime. This regime provides information related to the chain statistics inside the coil, which can reveal details about stiffness and self-avoiding behavior. In the regime, the results given by two equations are closer as the polymerization degree increases. One can expect that on a scale which is large enough, they will again appear as flexible coils.
A standard way to examine the structure of a polymeric chain at all distances is through the static structure factor S ( q ) , which is defined as
S ( q ) = ω ( q ) + ρ h ( q ) = ω ( q ) + ρ [ g ( q ) 1 ]
and the results for PS chains at different degrees of polymerization are presented in Figure 5. As can be seen, the general shape of the plots could be divided into three distinct regions. In the small wave vector regime ( q 0.2 ), S ( q ) is very sensitive to the degree of polymerization, both in magnitude and functional form. Increase of polymerization degree leads to enhanced S ( q ) . The qualitative behavior of S ( q ) at small wave vectors is controlled by a competition between the intramolecular fluctuation and the intermolecular correlation. As the wave vector increases, there is a plateau-like regime ( 0.2 q 0.4 ) for S ( q ) (self-similar structure of the chains). This regime increases as the degree of polymerization increases. The plateau regime of high wave vectors is compound in nature, reflecting both the fixed bond length constraint and local intermolecular structure in the melt. In the high wave vector regime ( q 0.4 ), all the different PS melts have the same structure, i.e., the same S ( q ) .
Figure 6 shows the predicted gyration radii of PEO and PMMA chains in comparison respectively with the simulation [34] and experimental ones [49]. In Figure 6a, the theoretical predictions are generally in good agreement with the reported values. The results indicate that the current model is quantitatively reliable for describing the scaling properties in a broad range of polymeric systems. The scaling exponents for the two polymers with low polymerization are fitted in Figure 6b. The values are 0.59 and 0.56, which are smaller than the exponent for low-polymerized PS chains. The scaling exponents give an indication of stretching ratios. Among these three polymer melts, PS displays the maximum stretching ratio, whereas PEO has the minimum one. In this regard, the stretching ratio increases with enhancing chain stiffness.
Figure 7 gives the gyration radii of PS chains with various packing fractions ( η ). Calculations are carried out for η in the range of 0.04 to 0.40. The highest packing mimics neat melt condition, while the lowest is more representative of a dilute solution. At low concentrations, the chains are extended but not rodlike. As the concentration increases, the size of PS chains decreases because of the greater screening of the intramolecular interactions. Another major result is that the change in averaged chain dimension is on the order of ~10% due to an incomplete cancellation of the long-range intra-chain excluded volume and condensed-phase-induced interactions. These modest nonideal effects are of the same qualitative size as discovered in diblock and triblock copolymer melts [50,51].
Finally, we consider the statistical properties of relatively long chains (polymerization degree m = 100 ) immersed in a monodisperse melt of shorter, chemically identical chains (polymerization degree m 1 = 1 ~ 100 ). Such systems can act as a fundamental test of the standard model for polymer conformation. Figure 8 describes the gyration radius of a long chain as a function of chain length of a short-chain molecule. If m is much larger than m 1 , one sees that the short chain molecules acts as a good solvent and the long chains are swollen, owing to the excluded volume effect. As m 1 increases, the swelling extent of long chains decreases, since the excluded volume effects are screened by the surrounding polymer chains. The expansion of long chains becomes insignificant as the segment number of shorter chains is larger than a certain value ( m 1 = 40 ). In the early theoretical investigations [9,15], it was predicted that polymer chains (with segment number m ) can be expanded in its homologue (with segment number m 1 ) in the range of m 1 < m 1 / 2 . However, these investigations are based on an ideal chains model. Actual polymer chains have larger excluded volume. From the present theoretical model, we can find that the length range of m 1 has been enlarged for the stretching of PS chain when compared to the ideal chain blends.

4. Conclusions

We present a density functional approach to study the microscopic conformation of polymer chains. The form factors, structure factors, and gyration radii have been calculated based on the cited force field parameters of the coarse-grained model. Firstly, the effects of chain length to the scaling properties of PS, PMMA, and PEO melts have been compared to clarify the contribution of chain stiffness to the stretching ratios. The scaling exponents are 0.63, 0.59, and 0.56, respectively. The results reveal that polymer chains display non-universal behaviors due to their different stiffnesses. On the other hand, the effects of concentration on the stretching of PS chains have been calculated to evaluate the influence of medium. As the packing fraction of PS chain increases, the gyration radius of low polymerized PS declines. Finally, the swelling of PS chains blended with shorter PS chains has been analyzed, and an enlarged expansion regime has been observed, which is different from the ideal chain blends.

Acknowledgments

We would like to express our sincere thanks to the National Natural Science Foundation of China (Grant No. 51525301) for financial support. This work was also supported by the National Natural Science Foundation of China (Nos. 21476007 and 51221002), the Doctoral Science Research Foundation of the Education Ministry of China under Grant (No. 20130010110005), and by Chemcloudcomputing of Beijing University of Chemical Technology.

Author Contributions

Ming Tian and Jianguo Mi proposed the concept of the paper and developed the model. Zhaoyang Wei and Nanying Ning performed the theoretical computations. Liqun Zhang derived the constitutive model. Zhaoyang Wei and Jianguo Mi analyzed the results and wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hara, Y.; Mayama, H.; Yamaguchi, Y.; Takenaka, Y.; Fukuda, R. Direct observation of periodic swelling and collapse of polymer chain induced by the Belousov–Zhabotinsky reaction. J. Phys. Chem. B 2013, 117, 14351–14357. [Google Scholar] [CrossRef] [PubMed]
  2. Doi, M.; Edwards, S.F. The Theory of Polymer Dynamics; Oxford University Press: Oxford, UK, 1988. [Google Scholar]
  3. Doi, M. Introduction to Polymer Physics; Oxford University Press: Oxford, UK, 1996. [Google Scholar]
  4. Chen, Q.; Kooij, E.S.; Sui, X.; Padberg, C.J.; Hempenius, M.A.; Schön, P.M.; Vancso, G.J. Collapse from the top: Brushes of poly(N-isopropylacrylamide) in co-nonsolvent mixtures. Soft Matter 2014, 10, 3134–3142. [Google Scholar] [CrossRef] [PubMed]
  5. Willott, J.D.; Humphreys, B.A.; Murdoch, T.J.; Edmondson, S.; Webber, G.B.; Wanless, E.J. Hydrophobic effects within the dynamic pH-response of polybasic tertiary amine methacrylate brushes. Phys. Chem. Chem. Phys. 2015, 17, 3880–3890. [Google Scholar] [CrossRef] [PubMed]
  6. Sui, X.; Chen, Q.; Hempenius, M.A.; Vancso, G.J. Probing the collapse dynamics of poly(N-isopropylacrylamide) brushes by AFM: Effects of co-nonsolvency and grafting densities. Small 2011, 7, 1440–1447. [Google Scholar] [CrossRef] [PubMed]
  7. Polotsky, A.A.; Plamper, F.A.; Borisov, O.V. Collapse-to-swelling transitions in pH-and thermoresponsive microgels in aqueous dispersions: The thermodynamic theory. Macromolecules 2013, 46, 8702–8709. [Google Scholar] [CrossRef]
  8. Noriega, R.; Rivnay, J.; Vandewal, K.; Koch, F.P.; Stingelin, N.; Smith, P.; Toney, M.F.; Salleo, A. A general relationship between disorder, aggregation and charge transport in conjugated polymers. Nat. Mater. 2013, 12, 1038–1044. [Google Scholar] [CrossRef] [PubMed]
  9. Lang, M.; Rubinstein, M.; Sommer, J.-U. Conformations of a long polymer in a melt of shorter chains: Generalizations of the flory theorem. ACS Macro Lett. 2015, 4, 177–181. [Google Scholar] [CrossRef] [PubMed]
  10. Heyda, J.; Muzdalo, A.; Dzubiella, J. Rationalizing polymer swelling and collapse under attractive cosolvent conditions. Macromolecules 2013, 46, 1231–1238. [Google Scholar] [CrossRef]
  11. Crawford, M.; Smalley, R.; Cohen, G.; Hogan, B.; Wood, B.; Kumar, S.; Melnichenko, Y.B.; He, L.; Guise, W.; Hammouda, B. Chain conformation in polymer nanocomposites with uniformly dispersed nanoparticles. Phys. Rev. Lett. 2013, 110, 196001. [Google Scholar] [CrossRef] [PubMed]
  12. McEwan, M.E.; Egorov, S.A.; Ilavsky, J.; Green, D.L.; Yang, Y. Mechanical reinforcement of polymer nanocomposites: Theory and ultra-small angle X-ray scattering (USAXS) studies. Soft Matter 2011, 7, 2725–2733. [Google Scholar] [CrossRef]
  13. Abbott, L.J.; McDermott, A.G.; Del Regno, A.; Taylor, R.G.; Bezzu, C.G.; Msayib, K.J.; McKeown, N.B.; Siperstein, F.R.; Runt, J.; Colina, C.M. Characterizing the structure of organic molecules of intrinsic microporosity by molecular simulations and X-ray scattering. J. Phys. Chem. B 2012, 117, 355–364. [Google Scholar] [CrossRef] [PubMed]
  14. Jouault, N.; Dalmas, F.; Said, S.; Di Cola, E.; Schweins, R.; Jestin, J.; Boué, F. Direct measurement of polymer chain conformation in well-controlled model nanocomposites by combining SANS and SAXS. Macromolecules 2010, 43, 9881–9891. [Google Scholar] [CrossRef]
  15. Flory, P.J. The configuration of real polymer chains. J. Chem. Phys. 1949, 17, 303–310. [Google Scholar] [CrossRef]
  16. Lifshitz, I.; Grosberg, A.Y.; Khokhlov, A. Some problems of the statistical physics of polymer chains with volume interaction. Rev. Mod. Phys. 1978, 50, 683–713. [Google Scholar] [CrossRef]
  17. Paturej, J.; Milchev, A.; Egorov, S.A.; Binder, K. Star polymers confined in a nanoslit: A simulation test of scaling and self-consistent field theories. Soft Matter 2013, 9, 10522–10531. [Google Scholar] [CrossRef]
  18. Grosberg, A.Y. Annealed lattice animal model and Flory theory for the melt of non-concatenated rings: Towards the physics of crumpling. Soft Matter 2014, 10, 560–565. [Google Scholar] [CrossRef] [PubMed]
  19. Curro, J.G.; Schweizer, K.S. Theory of polymer melts: An integral equation approach. Macromolecules 1987, 20, 1928–1934. [Google Scholar] [CrossRef]
  20. Honnell, K.G.; McCoy, J.D.; Curro, J.G.; Schweizer, K.S.; Narten, A.H.; Habenschuss, A. Local structure of polyethylene melts. J. Chem. Phys. 1991, 94, 4659–4662. [Google Scholar] [CrossRef]
  21. Schweizer, K.S.; Curro, J.G. Integral equation theories of the structure, thermodynamics, and phase transitions of polymer fluids. Adv. Chem. Phys. 1997, 98, 1–142. [Google Scholar]
  22. Pütz, M.; Curro, J.G.; Grest, G.S. Self-consistent integral equation theory for polyolefins: Comparison to molecular dynamics simulations and X-ray scattering. J. Chem. Phys. 2001, 114, 2847–2860. [Google Scholar] [CrossRef]
  23. Frischknecht, A.L.; McGarrity, E.S.; Mackay, M.E. Expanded chain dimensions in polymer melts with nanoparticle fillers. J. Chem. Phys. 2010, 132, 204901. [Google Scholar] [CrossRef] [PubMed]
  24. Tripathi, S.; Chapman, W.G. Microstructure and thermodynamics of inhomogeneous polymer blends and solutions. Phys. Rev. Lett. 2005, 94, 087801. [Google Scholar] [CrossRef] [PubMed]
  25. Wu, J. Density functional theory for chemical engineering: From capillarity to soft materials. AIChE J. 2006, 52, 1169–1193. [Google Scholar] [CrossRef]
  26. Wu, J.; Li, Z. Density-functional theory for complex fluids. Annu. Rev. Phys. Chem. 2007, 58, 85–112. [Google Scholar] [CrossRef] [PubMed]
  27. Jain, S.; Dominik, A.; Chapman, W.G. Modified interfacial statistical associating fluid theory: A perturbation density functional theory for inhomogeneous complex fluids. J. Chem. Phys. 2007, 127, 244904. [Google Scholar] [CrossRef] [PubMed]
  28. Jain, S.; Jog, P.; Weinhold, J.; Srivastava, R.; Chapman, W.G. Modified interfacial statistical associating fluid theory: Application to tethered polymer chains. J. Chem. Phys. 2008, 128, 154910. [Google Scholar] [CrossRef] [PubMed]
  29. Bymaster, A.; Chapman, W.G. An i SAFT density functional theory for associating polyatomic molecules. J. Phys. Chem. B 2010, 114, 12298–12307. [Google Scholar] [CrossRef] [PubMed]
  30. Emborsky, C.P.; Feng, Z.; Cox, K.R.; Chapman, W.G. Recent advances in classical density functional theory for associating and polyatomic molecules. Fluid Phase Equilibr. 2011, 306, 15–30. [Google Scholar] [CrossRef]
  31. Yu, Y.-X.; Wu, J. Extended test-particle method for predicting the inter- and intramolecular correlation functions of polymeric fluids. J. Chem. Phys. 2003, 118, 3835–3842. [Google Scholar] [CrossRef]
  32. Bymaster, A. Molecular Modeling the Microstructure and Phase Behavior of Bulk and Inhomogeneous Complex Fluids. Ph.D. Thesis, Rice University, Houston, TX, USA, 2009. [Google Scholar]
  33. Milano, G.; Müller-Plathe, F. Mapping atomistic simulations to mesoscopic models: A systematic coarse-graining procedure for vinyl polymer chains. J. Phys. Chem. B 2005, 109, 18609–18619. [Google Scholar] [CrossRef] [PubMed]
  34. Lee, H.; de Vries, A.H.; Marrink, S.-J.; Pastor, R.W. A coarse-grained model for polyethylene oxide and polyethylene glycol: Conformation and hydrodynamics. J. Phys. Chem. B 2009, 113, 13186–13194. [Google Scholar] [CrossRef] [PubMed]
  35. Xia, W.; Mishra, S.; Keten, S. Substrate vs. free surface: Competing effects on the glass transition of polymer thin films. Polymer 2013, 54, 5942–5951. [Google Scholar] [CrossRef]
  36. Yu, Y.-X.; Wu, J. Structures of hard-sphere fluids from a modified fundamental-measure theory. J. Chem. Phys. 2002, 117, 10156–10164. [Google Scholar] [CrossRef]
  37. Chandler, D.; Weeks, J.D. Equilibrium structure of simple liquids. Phys. Rev. Lett. 1970, 25, 149–152. [Google Scholar] [CrossRef]
  38. Weeks, J.D.; Chandler, D.; Andersen, H.C. Role of repulsive forces in determining the equilibrium structure of simple liquids. J. Chem. Phys. 1971, 54, 5237–5247. [Google Scholar] [CrossRef]
  39. Schweizer, K.S.; Curro, J.G. Integral-equation theory of the structure of polymer melts. Phys. Rev. Lett. 1987, 58, 246–249. [Google Scholar] [CrossRef] [PubMed]
  40. Honnell, K.G.; Curro, J.G.; Schweizer, K.S. Local structure of semiflexible polymer melts. Macromolecules 1990, 23, 3496–3505. [Google Scholar] [CrossRef]
  41. Kovalenko, A.; Hirata, F. Self-consistent description of a metal–water interface by the Kohn–Sham density functional theory and the three-dimensional reference interaction site model. J. Chem. Phys. 1999, 110, 10095–10112. [Google Scholar] [CrossRef]
  42. Zhao, L.; Li, Y.-G.; Mi, J.; Zhong, C. Integral equation theory for atactic polystyrene melt with a coarse-grained model. J. Chem. Phys. 2005, 123, 124905. [Google Scholar] [CrossRef] [PubMed]
  43. Yethiraj, A.; Hall, C.K. Monte Carlo simulations and integral equation theory for microscopic correlations in polymeric fluids. J. Chem. Phys. 1992, 96, 797–807. [Google Scholar] [CrossRef]
  44. Ding, Y.; Kisliuk, A.; Sokolov, A. When does a molecule become a polymer? Macromolecules 2004, 37, 161–166. [Google Scholar] [CrossRef]
  45. Cifra, P.; Bleha, T. Free energy of deformation of the radius of gyration in semiflexible chains. Macromol. Theory Simul. 2007, 16, 501–512. [Google Scholar] [CrossRef]
  46. Cifra, P.; Benková, Z.; Bleha, T. Persistence lengths and structure factors of wormlike polymers under confinement. J. Phys. Chem. B 2008, 112, 1367–1375. [Google Scholar] [CrossRef] [PubMed]
  47. Tree, D.R.; Muralidhar, A.; Doyle, P.S.; Dorfman, K.D. Is DNA a good model polymer? Macromolecules 2013, 46, 8369–8382. [Google Scholar] [CrossRef] [PubMed]
  48. Koyama, H.; Yoshizaki, T.; Einaga, Y.; Hayashi, H.; Yamakawa, H. Scattering function of oligo- and polystyrenes in dilute solutions. Macromolecules 1991, 24, 932–937. [Google Scholar] [CrossRef]
  49. Tamai, Y.; Konishi, T.; Einaga, Y.; Fujii, M.; Yamakawa, H. Mean-square radius of gyration of oligo- and poly(methyl methacrylate)s in dilute solutions. Macromolecules 1990, 23, 4067–4075. [Google Scholar] [CrossRef]
  50. Almdal, K.; Rosedale, J.H.; Bates, F.S.; Wignall, G.D.; Fredrickson, G.H. Gaussian-to stretched-coil transition in block copolymer melts. Phys. Rev. Lett. 1990, 65, 1112. [Google Scholar] [CrossRef] [PubMed]
  51. Fried, H.; Binder, K. The microphase separation transition in symmetric diblock copolymer melts: A Monte Carlo study. J. Chem. Phys. 1991, 94, 8349–8366. [Google Scholar] [CrossRef]
Figure 1. Schematic of the polymer model used in this work. Here, a middle segment from a polymer chain (filled black sphere) is fixed at the origin. The density distributions of segments from the tethered fragments ( C and D ), free ( F ), and short ( S ) chain molecules are related to the intra- and intermolecular segment–segment correlation functions. In a homogeneous polymer melt, the short chain no longer exists.
Figure 1. Schematic of the polymer model used in this work. Here, a middle segment from a polymer chain (filled black sphere) is fixed at the origin. The density distributions of segments from the tethered fragments ( C and D ), free ( F ), and short ( S ) chain molecules are related to the intra- and intermolecular segment–segment correlation functions. In a homogeneous polymer melt, the short chain no longer exists.
Polymers 08 00121 g001
Figure 2. Comparison of the average (a) inter- and (b) intramolecular correlation functions of the PS chain obtained from theory and molecular simulations [42] at 413.2 K.
Figure 2. Comparison of the average (a) inter- and (b) intramolecular correlation functions of the PS chain obtained from theory and molecular simulations [42] at 413.2 K.
Polymers 08 00121 g002
Figure 3. (a) Gyration radius as a function of chain length for PS chains at 500 K; (b) Determination of the scaling exponents. The circles are calculated results and the solid line is the best linear regression of the circles.
Figure 3. (a) Gyration radius as a function of chain length for PS chains at 500 K; (b) Determination of the scaling exponents. The circles are calculated results and the solid line is the best linear regression of the circles.
Polymers 08 00121 g003
Figure 4. (a) Comparison between experimental data of the form factors and calculated ones for PS melt with the molecular weight 1 × 104 in θ solvent; (b) Form factor for different lengths of PS melt at 500 K. Solid lines correspond to DFT calculations in Equation (21). Symbols show the curve obtained from the Debye expression.
Figure 4. (a) Comparison between experimental data of the form factors and calculated ones for PS melt with the molecular weight 1 × 104 in θ solvent; (b) Form factor for different lengths of PS melt at 500 K. Solid lines correspond to DFT calculations in Equation (21). Symbols show the curve obtained from the Debye expression.
Polymers 08 00121 g004
Figure 5. Structure factor of PS chain with three values of polymerization degree at 500 K.
Figure 5. Structure factor of PS chain with three values of polymerization degree at 500 K.
Polymers 08 00121 g005
Figure 6. (a) Gyration radius as a function of chain length for PEO and PMMA melts at 300 K; (b) Determination of the scaling exponents for the two polymers. The squares and circles are calculated results and the solid line is the best linear regression.
Figure 6. (a) Gyration radius as a function of chain length for PEO and PMMA melts at 300 K; (b) Determination of the scaling exponents for the two polymers. The squares and circles are calculated results and the solid line is the best linear regression.
Polymers 08 00121 g006
Figure 7. Gyration radius of PS chain as a function of packing fraction at 500 K. The chain length is fixed at m = 100 .
Figure 7. Gyration radius of PS chain as a function of packing fraction at 500 K. The chain length is fixed at m = 100 .
Polymers 08 00121 g007
Figure 8. Gyration radius of long chain PS ( m = 100 ) as a function of chain length of short chain PS ( m 1 = 1 100 ) at 500 K.
Figure 8. Gyration radius of long chain PS ( m = 100 ) as a function of chain length of short chain PS ( m 1 = 1 100 ) at 500 K.
Polymers 08 00121 g008
Table 1. Bond and nonbond coarse-grained force field parameters.
Table 1. Bond and nonbond coarse-grained force field parameters.
Speciesσ (Å)ε(k)θ(deg)L (Å)
PS a5.0862.51402.46
PEO b4.30405.81303.30
PMMA c6.50150.91222.80
a Reported in Reference [33]; b Reported in Reference [34]; c Reported in Reference [35].

Share and Cite

MDPI and ACS Style

Wei, Z.; Ning, N.; Zhang, L.; Tian, M.; Mi, J. Density Functional Theory of Polymer Structure and Conformations. Polymers 2016, 8, 121. https://doi.org/10.3390/polym8040121

AMA Style

Wei Z, Ning N, Zhang L, Tian M, Mi J. Density Functional Theory of Polymer Structure and Conformations. Polymers. 2016; 8(4):121. https://doi.org/10.3390/polym8040121

Chicago/Turabian Style

Wei, Zhaoyang, Nanying Ning, Liqun Zhang, Ming Tian, and Jianguo Mi. 2016. "Density Functional Theory of Polymer Structure and Conformations" Polymers 8, no. 4: 121. https://doi.org/10.3390/polym8040121

APA Style

Wei, Z., Ning, N., Zhang, L., Tian, M., & Mi, J. (2016). Density Functional Theory of Polymer Structure and Conformations. Polymers, 8(4), 121. https://doi.org/10.3390/polym8040121

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop