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 (
) of polymer chain to its length (
) via
. A prerequisite for application of these theories is that
at θ-state are known, which serves as an essential reference state. Although the scaling exponent at θ-state has been determined as
, 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
in a symmetric external field due to the fixed segment, as depicted in
Figure 1.
The free molecules are represented by
, composed of
segments; the chemically identical short chain molecules are represent by
, composed of
segments; while the tethered fragments are represented by
and
, composed of
and
segments, respectively. At equilibrium, the density distributions of free chains, short chains, and the tethered fragments satisfy the variational relations
where
stands for the grand free-energy,
(
) is the density of the
segment on chain
at position
. For a homogeneous polymer melt, the contribution of the
component is automatically deleted.
The segment distributions of the free () and short () chains around the fixed segment are related to the intermolecular site–site correlation function , where is the density profile of segments on molecule or around the fixed segment , and is the bulk density of segment. Whereas the distribution of segments from the fragments and are related to the intramolecular correlation function , where is the density profile of segment on the tethered chain ( or ) around the fixed segment . There is only one tethered polymer chain, therefore the normalization condition is . 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 , and the average intramolecular correlation function with .
The system considered above is equivalent to a mixture of four polymeric components
in a symmetric external field due to the fixed segment. The grand potential is related to the Helmholtz free-energy
through the Legendre transform
where
is the chemical potential of that segment
,
is the external field acting on that segment
. The first summation is over all chains
in the mixture
, and the second summation is over all segments of chain
. 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
, given by
; here segments
and
are nearest neighbors from the same molecule, and
is the Dirac-delta function. The bond length (
), 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
. 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]
where
is the free-energy density, which stems from the modified fundamental theory [
36], including both the scalar and vector contributions
where
with
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
where
is the interaction potential between any two species
and
, and can be represented by a cut-and-shifted LJ potential with a Weeks–Chandler–Anderson separation [
37,
38].
The functional
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
with the approximation
. Here
and
denote the direct correlation functions of semiflexible and flexible polymer chains. They are calculated from the PRISM integral equation [
39]
where
is the total correlation function and
is the intramolecular correlation function.
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
bonded monomers, enforced by giving each segment a label and allowing segments to exclusively bond to their specific matching segments. It is indicated as
in which the first summation is over all the segments
, and the second is over all the association sites on segment
as
, representing the set of all associating sites on segment
.
is the fraction of segment
that are not bonded at their association site A.
where
denotes the neighboring segment that will bond to segment
, and
. Here
is a geometric constant that accounts for the volume available for bonding between segments, and
represents the association Mayer
f-function.
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
. As a result, the Euler–Lagrange Equation is written as
where
is the set of all segments bonded to segment
. This equation can be rewritten to give the density profile
with
In Equation (11),
is the bulk chemical potential of chain
. The multiple integrals
and
for the free and short chains
are solved in a recursive fashion and are given by
where
is the Heaviside step function. When applied to the tethered chains
,
and
can be expressed as
The chemical potentials of the tethered fragments
and
can be determined using the normalization conditions
and
, where
. For the first tethered segment (segment “1” of
),
which yields
. Using other segments yields equivalent results for
. Substituting this into Equation (11), we obtain
which matches the known condition for the tethered segment. Solving for the other segments
gives
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, (if is even) or (if is odd) calculations are required for predicting the detailed local structures of polymer chains consisting of 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 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 . 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
Å 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
can be calculated from the intramolecular correlation function [
40,
43]
to quantitatively analyze polymer conformation.
Figure 3a shows
for PS chains at
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
. 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]
where
is the distance between segments
and
, and
denotes the magnitude of the scattering wave vector
. The expression for the form factor of a free chain is the well-known Debye equation
Figure 4a shows the form factors for PS melt with molecular weight 1 × 10
4 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 (
). 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
, which is defined as
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 (
),
is very sensitive to the degree of polymerization, both in magnitude and functional form. Increase of polymerization degree leads to enhanced
. The qualitative behavior of
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 (
) for
(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 (
), all the different PS melts have the same structure,
i.e., the same
.
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
) immersed in a monodisperse melt of shorter, chemically identical chains (polymerization degree
). 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
is much larger than
, 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
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 (
). In the early theoretical investigations [
9,
15], it was predicted that polymer chains (with segment number
) can be expanded in its homologue (with segment number
) in the range of
. 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
has been enlarged for the stretching of PS chain when compared to the ideal chain blends.