# Renormalized One-loop Theory of Correlations in Polymer Blends

###### Abstract

The renormalized one-loop theory is a coarse-grained theory of corrections to the self-consistent field theory (SCFT) of polymer liquids, and to the random phase approximation (RPA) theory of composition fluctuations. We present predictions of corrections to the RPA for the structure function and to the random walk model of single-chain statics in binary homopolymer blends. We consider an apparent interaction parameter that is defined by applying the RPA to the small limit of . The predicted deviation of from its long chain limit is proportional to , where is chain length. This deviation is positive (i.e., destabilizing) for weakly non-ideal mixtures, with , but negative (stabilizing) near the critical point. The positive correction to for low values of is a result of the fact that monomers in mixtures of shorter chains are slightly less strongly shielded from intermolecular contacts. The depression in near the critical point is a result of long-wavelength composition fluctuations. The one-loop theory predicts a shift in the critical temperature of , which is much greater than the predicted width of the Ginzburg region. Chain dimensions deviate slightly from those of a random walk even in a one-component melt, and contract slightly with increasing . Predictions for and single-chain properties are compared to published lattice Monte Carlo simulations.

## I Introduction

Equilibrium properties of polymer mixtures that involve correlations over distances greater than a few Angstroms are rather well described by a family of closely related self-consistent field theories. Among these are Flory-Huggins (FH) theory for the free energy of a homogeneous mixture and the random phase approximation (RPA) theory of composition fluctuations. Both of these are special cases of a more general self-consistent field theory (SCFT) of inhomogenous polymer liquids.

Theoretical treatments of polymer mixtures that go beyond SCFT may be
roughly classified as either microscopic liquid state theories
or coarse-grained theories of correlation effects.
Liquid state theories such as the polymer reference interaction site
model (PRISM)
^{1}
and
Born-Green-Yvon theories
^{2}
are designed to treat models with strong repulsive interactions, and
attempt to predict all details of liquid structure from first principles.
Coarse-grained theories are based on more idealized models of chain
conformations and interactions, e.g., as continuous Gaussian chains
with point-like interactions, similar to those used in coarse-grained
SCFT.

Over the past 20 years, a
series of coarse-grained theories of correlation effects in
polymer mixtures have made use of a functional integral
representation of the partition function that was originally
introduced by Edwards.
^{3, 4, 5, 6, 7, 8, 9}
All of these ”field theoretic” studies of mixtures have thus
far used an approximation in which the distribution of field
fluctuations is approximated by a Gaussian functional. This
approximation is referred to in the jargon of field theory
(which we have adopted) as a “one-loop” approximation. An
alternative way of deriving this theory, which we find more
informative, is discussed in Sec. (VI) of
this paper.

### i.1 Coarse-Graining and Renormalization

The reason for studying correlations in coarse-grained
models is the hope that, by doing so, one might describe
phenomena that are not captured by SCFT, but that are, in
some sense, independent of short-wavelength details. The
first several studies of coarse-grained models of blends
^{3, 4, 5, 6, 7, 8} were, however,
plagued by a sensitivity of all predictions to the treatment
of very short wavelength correlations, at length scales which
the coarse-grained Hamiltonian was never intended to describe.
The way around this problem, which was first applied in this
context by Wang ^{9}, is a renormalization procedure
in which the sensitivity to monomer scale structure is
absorbed into the values of a few phenomenological parameters.
Wang ^{9} considered the one-loop prediction for the
free energy of a homogeneous polymer blend, and showed that the
sensitivity of this quantity to local fluid structure could be
absorbed into an appropriate renormalization of the Flory-Huggins
interaction parameter. More recently, we ^{10}
showed how a generalization of this procedure can be used to
obtain unambiguous results for one-loop corrections to both
collective and single-chain correlation functions at nonzero
wavenumbers, using the same formalism for polymer blends and
diblock copolymer melts.

Development of this renormalization procedure required two insights that were absent from the interpretation of earlier studies of the one-loop theory:

The first was that the SCF approximation should be identified with the infinite chain limit of the one-loop theory, rather than a simple random-mixing approximation. Prior to Wang’s work, the SCF approximation had often been conceptually identified with a simple mean-field approximation that is obtained from a saddle-point approximation for the underlying functional integral. This mean-field theory ignores all correlations except those arising from intramolecular connectivity, which is never a quantitatively useful approximation. Defining SCF parameters by a process of extrapolation to the infinite chain limit instead takes advantage of the fact that local correlations in liquids of long polymers depend only weakly on chain length, even in strongly correlated liquids. The difference between the free energy of a system of finite chains and that of a hypothetical reference system of infinite chains is largely the result of slight differences in long-wavelength correlations. These differences may be adequately described by a coarse-grained model. By treating the interaction free energy of an infinite-chain reference system as an input to the theory (which must be inferred from experimental or simulation data), one can obtain a renormalized theory that makes unambiguous predictions of corrections to SCF phenomenology.

A second necessary insight was that the one-loop prediction for
the structure factor at nonzero wavenumber is affected
by a renormalization of statistical segment length, as well as
of the interaction free energy. The value of the statistical
segment length of a flexible polymer in a dense mixture is
generally different than the “bare” value required as an input
to the theory. The one-loop predictions for the magnitude of
this difference ^{8}, are sensitive to short-wavelength
correlations. This renormalization of has no effect upon the
predicted value of structure factor at exactly , and
so did not need to be taken into account in Wang’s interpretation
of one-loop predictions for . It becomes an unavoidable
complication, however, in calculations of single-chain correlations,
and in applications to diblock copolymer melts, for which the
dominant fluctuations occur at nonzero .

In our previous paper ^{10}, we focused primarily
on clarifying the physical reasoning underlying our renormalization
procedure, and on demonstrating that it yields a physically consistent
interpretation. In this paper, we present predictions for binary
homopolymer blends. Predictions for diblock copolymer melts will
be presented elsewhere.

### i.2 Related Prior Work

The analysis given here has several particularly close antecedents.

Our analysis of the free-energy density of a homogeneous blend and
the long wavelength limit of , which is given in Secs.
II and III, is closely related to
that of Wang.^{9} The main differences arise from
mathematical approximations that Wang introduced in order to
obtain a completely analytic approximation for behavior near the
spinodal. We have avoided all approximations other than the
one-loop approximation itself, at the cost of some additional
numerical analysis. The additional rigor has allowed us to calculate
some additional quantities (e.g., the absolute shift in the
critical point) and to identify physical phenomena that are
important far from the spinodal.

Our analysis of corrections to Gaussian chain statistics in blends,
given in Sec. IV, is closely related to the recent
theoretical analysis of chain statistics in one-component melts by
Beckrich, Johner, Semenov, and coworkers ^{11, 12}.
Like us, but prior to our work, these authors introduced a renormalized
statistical segment length, defined by extrapolation to ,
as part of a systematic calculation of small corrections to random-walk
statistics. We initially failed to notice the connection to our own
work, ^{10} or to cite their work. These authors
presented a beautiful analysis of universal corrections to Gaussian
statistics in an incompressible melt at intermediate wavelengths,
less than the coil size but greater than a monomer size. They also
obtained analytic results for one-loop corrections to the intramolecular
correlation function at all wavenumbers in melts of equilibrium
polymers, with an exponential distribution of chain lengths. By
introducing a numerical procedure for extracting UV convergent parts
of divergent integrals, we are also able to numerically evaluate
one-loop predictions for melts and blends of monodisperse homopolymers
and (elsewhere) for block copolymer melts. Here, we consider
statistics of monodisperse polymers in binary blends.

The most important antecedent of both our work and the recent work
of the Strasbourg group is the much larger, and (generally) more
sophisticated theoretical literature on correlations in semi-dilute
polymer solutions.^{13} At a technical level, the
closest analog of our approach within this literature seems to
be the analysis of mono-disperse semi-dilute solutions by Ohta
and Nakanishi.^{14, 15}

### i.3 Outline

The paper is organized as follows: Sec. II provides a brief review of the renormalized one-loop theory for the free energy density, and of our notation. Sec. III presents predictions for corrections to the RPA theory of composition fluctuations. Sec. IV presents predictions for corrections to random walk statistics in polymer blends. Sec. V presents comparisons of theoretical predictions to results of published lattice Monte Carlo simulations. Sec. VI discusses the physical relationship between the one-loop theory for the free energy and RPA theory of correlations. Sec. VII clarifies the physical origin of some features of the one-loop theory by relating them to corresponding features of the underlying RPA theory.

## Ii Theoretical Formulation

### ii.1 Model and Notation

We consider a coarse-grained model of a binary homopolymer blend. Let
be the average volume per monomer, for either species. Let the volume
fraction, degree of polymerization, statistical segment length and packing
length ^{16} for species be denoted by , ,
, and , respectively, for species
and 2.

Let the potential energy be the sum of an intramolecular potential and a monomer-monomer pair interaction

(1) |

Here, is the concentration of monomers of type , in which is the position of monomer on molecule of species . The strength of the pair interaction is taken to be of the form

(2) |

Here, is a function with a range and an integral , so that approaches a Dirac -function as . Here, is a dimensionless compression energy and is a dimensionless “bare” interaction parameter. The incompressible limit may obtained by considering the behavior when is very large.

Let

(3) |

where , and let denote the Fourier transform. Let be the intramolecular contribution to from chains of type . A continuous random walk model yields

(4) |

where is the radius of gyration for species and is the Debye function.

Composition fluctuations in an effectively incompressible blend may be described by a scalar correlation function

(5) |

Here, we follow Schweizer and Curro and our previous work, by
expressing this scalar correlation function in terms of what we
(and Wang ^{9}) refer to as an “apparent” interaction
parameter , by defining

(6) |

This is, of course, a generalization of the RPA expression for , in which it is assumed that chains are random walks and that is independent of , , and .

### ii.2 One-loop Free Energy

The Helmholtz free energy per volume for this coarse-grained model may be expressed, in the incompressible limit, as a sum

(7) |

where

(8) |

is the ideal Flory-Huggins entropy of mixing, and

(9) | |||||

is a simple mean-field approximation for the interaction energy density. The remainder, denoted , is a correlation free energy. Here, and throughout the remainder of this paper except Sec. V, where we compare to simulation results, we set .

A one-loop approximation for has been obtained in previous work by considering a Gaussian approximation for the distribution of fluctuations of a complex chemical potential field within the context of the Edwards auxiliary field formalism. In Sec. VI, we discuss an alternative derivation in which this approximation is instead obtained by thermodynamic integration of the RPA correlation function. However it is obtained, the resulting approximation for is given, in the general case of a compressible liquid, by a Fourier integral

(10) |

Here, the bold faced symbols and indicate matrices with elements and , respectively, denotes the identity, and is a determinant. Here and hereafter, we use the shorthand for Fourier integrals.

In the nearly-incompressible limit of interest here, Eq. 10 for may be approximately reduced to

(11) |

To obtain the integrand of Eq. (11), in a model with a pair interaction of the form given in Eq. (2), we have both taken the limit in Eq. (10), and also replaced the Fourier transform of by unity. This yields an expression for the integrand that is valid for . In what follows, we regularize the one-loop approximation for by treating as a cutoff wavenumber, and simply restricting the integral in Eq. (11) to . This is equivalent to the use of a pair potential in which the Fourier transform of is taken to be unity for , and zero for . The purpose of the renormalization procedure discussed below is to extract predictions that are independent of such details of the treatment of short-wavelength correlations.

### ii.3 Divergence and Renormalization

It has long been known that, in a model of continuous Gaussian
chains, the idealized limit of
point-like interactions yields a one-loop expression for
that diverges as the cube of .
^{3, 4, 5, 9},
In previous work, both Wang ^{9} and we ^{10}
have analyzed the structure of the high- divergence of this free
energy density, and its physical significance.
We found that Eq. (11) for can be
expressed as a sum

(12) |

of an “ultra-violet” (UV) divergent contribution
, which increases with increasing ,
and a contribution that is independent of . The UV
divergent contribution may be expressed as
a sum
^{10}

(13) |

where

(14) |

and where . The contribution
, which is proportional to , was identified
^{10} as an excess free energy for chain ends. The
quantity , which is independent of , is
instead a correction to the “bulk” interaction free energy density
of a hypothetical reference system of infinitely long chains.

To remove the explicit cutoff dependence, we re-interpret the one-loop free energy (excluding end effects) as a sum

(15) |

in which

(16) |

is identified as the interaction free energy of the relevant form of SCFT. The contribution is a non-universal chain end contribution that could (and generally should) be added to SCFT, to reflect the fact that the fluid structure is perturbed in the immediate vicinity of any chain end. The remaining cutoff-independent contribution , is identified as a universal correction to the SCF phenomenology.

Eq. (16) for simplifies considerably in the case of structurally symmetric polymers, with . In this case is given, to within a composition-independent constant, by a function of the form

(17) |

where

(18) |

is an effective Flory-Huggins interaction parameter. This UV
divergent expression for was first obtained by Olvera de la
Cruz et al. ^{3}. The physical origin of the
UV divergence of this quantity is discussed in Sec. VII.

We focus in the remainder of this paper on the UV convergent correction . It is convenient to introduce a notation for and related quantities, in which we use to denote a renormalized Fourier integral. The renormalized Fourier integral of an integrand is defined by a difference

(19) |

between the unrenormalized integral and the UV
divergent part of the integral, denoted by .
In a theory of continuous chains, the “UV divergent” part of Fourier
integral is obtained by constructing a high- asymptotic expansion
of the integrand, in decreasing powers of , and retaining only
those terms in the expansion that lead to a UV divergent integral.
^{10}

With this notation, we may write

(20) |

The -dependence of may be isolated by non-dimensionalizing the renormalized integral in terms of a reference length , which we take to be the root-mean-squared end-to-end of species 1. It may be shown that all of the quantities in the integrand can then be expressed as functions of and of a set of dimensionless parameters , , , . By this method, we find that may be expressed as a function of the form

(21) |

in which is a dimensionless function given by the value of the non-dimensionalized renormalized integral. The required set of dimensionless parameters is the same as that required by the standard non-dimensionalized SCFT for an incompressible blend. Note that, in Eqs. (20) and (21), we have written as a product of a parameter , but have not specified whether by we mean , , or (perhaps) the low- limit of . In fact, different choices lead to different variants of the theory, as discussed in more detail below.

Results for and related quantities that are reported here have all been obtained by numerically evaluating the relevant renormalized Fourier integrals. To do this, for each quantity of interest, we first numerically evaluate the unrenormalized integral with a large but finite cutoff , and then subtract an analytic result for the UV divergent part of the regularized integral. The Fourier integrals required by the one-loop theory may all be expressed as either one dimensional integrals with respect to a wavenumber (for free energy density and its derivatives) or two-dimensional integrals with respect to and a polar angle (for correlation functions at nonzero wavenumber). Subtraction of the UV divergent part from an unrenormalized integral yields an estimate of the renormalized integral that depends only weakly on the numerical cutoff . Our final result for the value of a renormalized integral is obtained by repeating this procedure for several values of , and numerically extrapolating to .

### ii.4 Structure Function

The long-wavelength limit of is related to the composition dependence of the free energy density by a statistical mechanical theorem:

(22) |

The quantity is also related to the long-wavelength limit of by

(23) |

Here, and hereafter, we use to denote the limit of the quantity defined in Eq. (6). The quantity is the spinodal value of , given by

(24) |

Using decomposition (15) for , and hereafter neglecting the explicit end effects, the one-loop approximation for may be expressed as a sum

(25) |

in which

(26) |

is the interaction parameter of phenomenological SCFT, and

(27) |

is a correction to SCFT arising from long-wavelength fluctuations.

An explicit Fourier integral expression for may be obtained
by differentiating the integrand of Eq. (20) for
twice with respect to . The resulting approximation for
may be expressed in terms of either Edwards’ screened
interaction ^{10} or in terms of the correlation
function

(28) |

When expressed in terms of , it is found that

(29) |

where

(30) |

The quantities and are proportional to the third and fourth derivatives of with respect to , respectively. Since , and are independent of , the proximity to the spinodal in this expression is controlled by alone.

By non-dimensionalizing the integral in Eq. (29), we may show that it yields a correction to of the form

(31) |

in which is a non-dimensional function that is defined by the non-dimensionalized integral. Here,

(32) |

is an invariant degree of polymerization, in which (by convention) and .

In Eqs. (20), (21), and (31), we have intentionally expressed and as functions of an interaction parameter , without specifying whether this input parameter should be taken to be the bare parameter or some type of a renormalized value. In fact, different variants of the renormalized one-loop theory can be obtained by different choices for the interaction parameter used within these Fourier integrals. If we define by simply subtracting the UV divergent part of the one-loop expression for , we obtain a theory that is no longer UV divergent, but that retains an explicit dependence upon the unmeasurable “bare” parameter . Following the logic used to construct renormalized expansions of field theories, however, we may also replace by a renormalized parameter within each integrand to obtain a theory from which all reference to the bare parameter have been removed. This substitution may be formally justified, within the context of a one-loop approximation, by observing that the difference between approximations obtained by using a bare or renormalized parameter in the one-loop theory differ only at second order in a systematic loop expansion.

In the remainder of this paper, we follow Wang ^{9}
by considering a self-consistent one-loop approximation. In this
approximation, also known as a “Hartree” approximation,
is replaced within the Fourier integral expression for
by . In the resulting theory,
is given by a function of the form

(33) |

Here, is related to the SCF (i.e., ) interaction free energy by Eq. (26). For purposes of numerical evaluation, it is convenient to rewrite Eq. (33) as an explicit expression for as an explicit function of .

## Iii Collective Fluctuations

### iii.1 Symmetric Blends, Critical Composition

Figs. 1 displays predictions of the renormalized Hartree theory for in blends of structurally symmetric polymers ( and ) at their critical composition (). It shows predictions for as a function of blends with three different chain lengths, corresponding to = 64, 128, and 512. The difference between from is largest for the shortest chains, and vanishes in the limit . The SCF predicts a critical point for a symmetric . The actual critical point for such a blend is reached when .

Fig. 2 shows the underlying dimensionless function as a function of . Note that both and vanish in the limit for any structurally symmetric mixture, with and , because this is an exact result for such a ideal mixture, in which the two species are physically indistinguishable.

Notably, the dependence of upon is not monotonic: For small values of , the deviation is positive and increases with increasing . Closer to the critical point, begins to decrease with increasing , changes sign at , and is negative at the critical point. Note that all the curves in Fig. (1) intersect at , where . For any structurally symmetric blend, with and , we find that the initial slope is

(34) |

for any composition.

At the critical point of such a blend ( and ), the one-loop prediction for approaches a negative limiting value

(35) |

with a slope that diverges at the critical point. In the Hartree approximation, this yields a predicted critical value of

(36) |

If is a decreasing function of temperature, this corresponds to a depression in the critical temperature by an amount

The predicted decrease in near the critical point is not surprising: Near the critical point, long wavelength correlations decrease the level of local mixing between and monomers, thereby decreasing the apparent interaction parameter and stabilizing the homogeneous phase.

The origin of the positive values obtained for at small values of is less obvious. In Sec. VII of this paper we show that this prediction is the result of a subtle dependence of local correlations in a polymer liquid upon overall chain length, which is discussed in much greater detail in an accompanying paper. In a dense liquid of long polymers, the immediate vicinity of each monomer is crowded with monomers that belong to the same chain. This intramolecular “self-concentration” reduces the space available for monomers from other chains, thus creating a so-called correlation hole in the intermolecular correlation function. In a mixture, the correlation hole reduces the magnitude of the potential energy of interaction between monomers of types and . We show in the accompanying paper that the depth of the correlation hole, decreases slightly with decreasing chain length. For from the spinodal, this simple packing effect causes to increase with decreasing chain length, reflecting the fact that shorter chains are less strongly shielded from intermoecular contacts. Because we define to be the limiting value of in the limit of infinitely long chains, this effect tends makes a positive contribution to the difference for any finite . In the accompanying paper, we analyze this effect without reference to the one-loop theory, and derive the prefactor in Eq.34 by very different reasoning than that used here.

A shift in the critical temperature by an amount was predicted previously by Holyst
and Vilgis ^{6}. The one-loop theory of
Holyst and Vilgis did not use a systematic renormalization
procedure, but instead introduced an ad hoc cutoff
to avoid the UV divergence. This approach
captured the correct scaling, but does not allow a meaningful
calculation of numerical values for one-loop corrections,
because it yields corrections whose values are sensitive
to the exact numerical value chosen for an arbitrary cutoff
length.

### iii.2 Vicinity of the Spinodal

We now discuss the behavior of the Hartree theory near the spinodal.
This part of our analysis is similar to that given by Wang
^{9}, who focused on behavior near the spinodal, but who
made several approximations that we avoid.

A rigorous asymptotic expansion of the behavior of near the spinodal may be obtained by expanding , and in Eq. (29) around their values at and . Let

(37) |

The required asymptotic expansion of is

(38) |

where

(39) |

By also expanding and around , we obtain an asymptotic expansion of the form

(40) |

where , , and are dimensionless coefficients that are independent of , but that that depend upon , , and . Here, is an arbitrary choice of a reference degree of polymerization, which could be taken to be or or an appropriate average of the two.

The coefficients and in this expansion can be calculated analytically, by considering the singular behavor of the integral near as . The coefficient is given

(41) | |||||

Here, is the third derivative with respect to
of the ideal mixture free energy per monomer .
This quantity vanishes at the Flory-Huggins critical composition,
at which , and so (as noted
by Wang ^{9}) also vanishes at the critical
composition . At the critical composition, we find

(42) |

These values for and agree with corresponding results of
Wang^{9}.

In the Hartree approximation, the value of (the reduced “temperature”) at the critical point is determined by the value of the constant at the critical composition: The critical point occurs when at , or when

(43) |

If is nonzero, an expansion of this form thus yields a shift in the critical value of of order .

The coefficient cannot be obtained from an asymptotic analysis of the low- behavior of the integral in the Eq.(29). To understand why, recall that the problem actually involves three length scales: the monomer size (or cutoff length ), the coil size , and correlation length . Very near the spinodal, . Our renormalization procedure removes all dependence on , but leaves a dependence on both and . Away from the critical composition, the integral in Eq. (29) develops an infra-red (IR) divergence of the form as . This divergence is the result of fluctuations with wavenumbers , and so its prefactor can be obtained by considering the behavior of the integrand near . At the critical composition, the coefficient can be isolated by considering the IR divergence of the derivative of Eq. (29) with respect to , which then diverges as . The actual value of , however, convergences even at the critical point, and is dominated by fluctuations with . As a result, the one-loop prediction for the shift in the critical point can only be calculated by numerically evaluating the renormalized integral, as done here, without introducing approximations for and that are valid only for .

At the critical composition, the asymptotic expansion of the Hartree theory may be expressed more compactly as a sum

(44) |

where , with . The theory exhibits strongly non-classical critical behavior for values of less than a crossover value . Approximating by the value at which the fluctuation correction is equal to yields a reduced crossover temperature

(45) |

as first noted by de Gennes. ^{17}
As noted by Holyst and Vilgis ^{6},
the predicted with of the critical region
is thus much less than the predicted
magnitude of the shift of the critical temperature from
its SCF value.

Very close to the critical point, , the Hartree theory yields , or . The predicted critical exponent of is much larger than the known Ising critical exponent of . This is a well known defect of the Hartree theory.

### iii.3 Off-Critical Blends

At any composition except the critical composition, the Hartree theory
considered here does not have a spinodal: At off-critical compositions, because
diverges like with decreasing , there is no value of
for which . As noted by Wang ^{9}, in the Hartree
approximation, this causes to become a non-monotonic function of
or (equivalently) to become a multi-valued function
. An example of this behavior is shown in Fig. 3,
which shows predictions of renormalized Hartree theory for vs.
in a series of off-critical blends of various compositions with , , and . The turning point at which
reaches its maximum value, which is also where is maximum, is what
Wang referred to as a “pseudo-spinodal”. We do not attach any physical
significance to this turning point, except as a particularly obvious failure of
the Hartree approximation.

i

## Iv Single Chain Properties

The one-loop theory ^{18, 10} also predicts
slight deviations from random-walk statistics for individual chains.
The one-loop approximation for the intramolecular correlations function
is a form of renormalized perturbation theory in which each pair of
monomers on a chain interact via Edwards’ screened interaction,
and in which the effect of this interaction is taken into account
to first order in powers of the strength of the screened
interaction^{19}.

We consider the single-chain correlation function . This is given for a Gaussian homopolymer by

(46) |

where is the radius of gyration
of species and is the Debye function.
The one-loop correction to in a homopolymer blend
is given by ^{18, 10}

(47) |

Here, is the screened interaction between any two monomers, which is given in an incompressible blend at wavelengths by

(48) |

The function is given by

where is the 4-point intramolecular correlation
function for an ideal Gaussian chain^{10}.
It is straightforward to verify that
vanishes in the limit , as required to retain consistency
with the requirement that .

If the high- behavior of the above integral is regularized introducing
a sharp cutoff at , and using Eq. (48)
for all , the resulting integral is found to increase
linearly with as is increased. We
^{10} and Beckrich et al. ^{11}
have shown, however, that this UV divergence can be absorbed into
renormalization of the statistical segment length. After this UV
divergent part is subtracted, the remaining UV convergent correction,
which we will denote by may be used to
characterize the deviation from random-walk statistics. More
precisely, is the deviation of the
intramolecular correlation function from that
of a Gaussian chain with a renormalized statistical segment length
characteristic of a dense one-component liquid of
infinitely long chains. The correction
is found to be smaller than by a factor proportional
to , and exhibits a nontrivial dependence on .

The fractional deviation of the polymer radius gyration is given by (appendix A)

(49) |

where is the prediction for a Gaussian chain with statistical segment length . This fractional deviation is proportional to .

### iv.1 Monodisperse Melt

We first consider the case of a monodisperse one-component melt, which
we obtain by setting . The one-loop approximation used
here reduces in the case of one-component liquids to that investigated
by Semenov et al.
^{20, 12, 21, 11}.
The only difference in our treatment of this case is our use of a
numerical integration and renormalization procedure to obtain
accurate results for monodisperse polymers.

Fig.4 shows the predicted correction to the single chain correlation function in a monodisperse melt. From the low behavior (see appendix A) we find that the predicted fractional change of the radius of gyration is

(50) |

We thus predict a contraction of finite chains, relative to a hypothetical
Gaussian chain with a statistical segment length , in
agreement with the conclusions of Ref. ^{12} [see Eq. (20)].

### iv.2 Bidisperse Melt

We next consider the behavior of a trace concentration of a species of chains of length (component 1) in a matrix of chemically homologous shorter chains of length (component 2). This limit is obtained by calculating the limit of . To represent chemically similar species, we set and . This case reduces to that of a mono-disperse melt when , and to a dilute solution of polymer in oligomeric solvent when .

In the previous subsection, we noted that the size of a single chain in the pure melt decreases as compared to the unperturbed value. It is also well known that a single chain swells in the dilute solutions with chemically similar small molecule solvents, which create a good solvent. In Fig.5, we show one-loop predictions for the fractional change of the radius of gyration of a test chain of length in a matrix of chains of length . The change is indeed negative for close to 1 and becomes positive for close to 0, passing through zero at .

The limit corresponds to a polymer solution in a good
small-molecule solvent. The curve in Fig.5 diverges as
near (appendix A). In this
limit, where the matrix chains are point-like, Eq.
(48) for the screened interaction between a
pair of monomers on a longer chain reduces to a constant
. Recall that the definition of a “monomer”
in an incompressible coarse-grained theory is simply the length of a
chain that occupies an arbitrary reference volume . If, to simplify
notation, the reference volume is taken to be equal to the excluded
volume of a solvent molecule (i.e., a matrix chain), by
setting , we obtain an interaction
characterized by an excluded volume parameter equal to the solvent
volume. Predictions of the theory in this limit are thus identical
to those obtained in the first order perturbation theory for a single
chain with a point-like interactions .^{22}
It is straightforward to show that this perturbative result is
consistent with the prediction ^{23} that the longer
chain will undergo substantial expansion if ,
signalling a breakdown of first order perturbation theory.

### iv.3 Binary Blends

We now consider binary blends with . For simplicity, we focus on symmetric systems, with and . Fig. 6 shows the one-loop predictions for the fractional change of the radius of gyration squared, rescaled by , as function of for various compositions. Increasing always causes both species in a blend to contract. The dependence of upon is strongest when species is a tracer in a matrix of (). There is no dependence of upon in the opposite limit of essentially pure homopolymer 1. In the dilute limit, , the one-loop theory predicts a deviation that is a strictly linear function of . Furthermore, it can be shown that the slope of the limiting line is -0.62, which is -2 times that of the coefficient to the divergence for the bidisperse melt in the limit of small-molecule good solvent discussed in the previous subsection (appendix A).

It is interesting to note that at the critical composition, , the deviation induced by increasing always remains less than the chain-length dependent deviation from ideal chain statistics that is already present in the one-component melt (). Also note that the radius of gyration is a smooth, nearly linear function of even near the critical point. This reflects the fact that the long-wavelength composition fluctuations that appear near the critical point have characteristic wavelengths much larger than the radius of gyration, and so couple very weakly to the conformations of individual chains.

## V Comparison to Simulation

Here, we compare the predictions of the one-loop theory to the results
of lattice Monte-Carlo simulation. We compare to bond-fluctuation model
(BFM) simulations of symmetric blends at Binder, Deutsch,
and Müller, and to BFM simulations of single-chain properties by
Beckrich et al. ^{11}.
The bond fluctuation model is a lattice model on a simple cubic lattice
in which each monomer occupies a set of 8 sites at the vertices of a
cube. In all of the simulations discussed here, the monomer concentration
is 1/16 monomer per lattice site, which is half the maximum possible
concentration in this model. The statistical segment length extrapolated
for the infinitely long chain is 3.244, which has been reported in ^{12}.
In simulations of polymer blends, monomers
of types and on neighboring sites interact with an energy ,
with and . Binder and coworkers have
carried out semi-grand canonical ensemble simulations using two
different variants of the rules for what constitute a “neighbor”:
One in which each monomer can interact only with monomers located at
the 6 closest allowed sites, and another in which each monomer can
interact monomers with monomers at any of the 54 closest allowed
sites.

### v.1 Determining

To compare predictions of a coarse-grained theory to simulation data,
one must somehow establish a relationship between the interaction
parameters used in the simulations and the phenomenological
interaction parameter used in both Flory-Huggins theory the one-loop
theory. In an accompanying paper, we propose a way to do this, which
we use here. The analysis given there refines and justifies a form of
“modified Flory-Huggins theory” that was originally proposed by
Müller and Binder.^{24}
There, we consider a perturbation theory for a lattice model of the
type described above, in which free energy of mixing is expanded in
powers of a perturbation parameter .
We identify the SCF interaction free energy per monomer
by considering the behavior of the free energy of mixing in the limit
of infinitely long chains. We show there that is given
to first order in a power series in by a function

(51) |

in which

(52) |

Here, is an “effective coordination number” that is obtained by evaluating a value of the average number of inter-molecular nearest neighbors per monomer in a one-component reference liquid, with , and extrapolating the value to the limit of infinitely long chains. The values of used here were reported by Müller and Binder: for the model with 54 interacting sites, and for the model with 6 interacting sites. A more complete discussion of the reasoning underlying this prescription is given elsewhere.

Our use of a first order perturbation theory to estimate in the one-phase region of the phase diagram is justified by the fact that the critical value of decreases as with increasing chain length. The fractional error in our estimate of arising from the use of a first order expansion is thus expected to be of order near the critical point. For , this error is thus expected to be much smaller than the corrections to predicted by the one loop theory, which are of order . The resulting errors in our estimate of the “bulk” interaction parameter are expected to be comparable in importance to the corrections that arise from end effects.

### v.2 Composition Fluctuations in Blends

Fig. 7 shows the results of Müller for
plotted vs. for bond-fluctuation model
simulations of chains of length ^{25}.
These simulation used a variant of the BFM in which interactions
extend over 54 neighboring sites, and in which .
Values of were extracted from the semi-grand equation
of state. We have calculated the SCF effective parameter
for each simulation from Eq. (52) using the
value of reported by Mueller and Binder.
^{24}

^{25}. The calculation of values of for simulation results is discussed in the text.

In the upper panel of Fig. 7 the solid line is the RPA prediction and the dashed line is the one-loop prediction obtained using the appropriate value of (The value was reported as 240, which were computed from the statistical segment length of chains). The lower shows a comparison of simulation results for the difference multiplied by to the universal curve shown in Fig. 2. We see that the self-consistent one-loop theory gives an extremely accurate description of deviations from SCFT far from the critical point, for , but appears to overestimate the stabilizing effect of critical fluctuations. For this system, the suppression of by long-wavelength fluctuations near the critical point is almost exactly cancelled by the enhancement of at smaller values of , giving a critical value of that is closer to the SCFT value than to the value predicted by the one-loop theory.

Fig. 8 compares the predictions of the
one-loop theory for the ration for
different values of to the results of simulations by
Deutsch and Binder ^{26}. Here,
is the
critical temperature that is predicted by SCF if we use
Eq. (52) for . The critical temperatures
reported by Deutsch and Binder were identified for each chain
length by a finite-size scaling analysis of simulations with
different simulation cell sizes. Simulation results are shown
for variants of the BFM model in which the neighbor interactions
extend over 6 neighbors (circles) and over 54 neighbors (triangles),
For both versions of the models,
extrapolates to unity as . This confirms
that first order perturbation theory prediction for is
indeed asymptotically exact in the limit .

For both versions of the BFM, however, the measured deviations of are significantly smaller than those predicted by the one-loop theory. This quantitative discrepancy must be, in part, a result of the inadequacy of the one-loop theory as a theory of critical phenomenon, and the need for a renormalization group approach near the critical point. This data suggests that, if it is indeed possible to express universal corrections to that arise from critical fluctuations as a power series in , the coefficient of any correction must be much smaller than predicted by the one-loop theory.

It is important to also note that the slight deviations of from unity in Fig. 8 are not the same for the two variants of the BFM, and appear to of opposite sign. This model-dependence cannot be so easily explained as a failure of the one-loop approximation near the critical point, since we expect that a more sophisticated treatment of critical phenomena in a simple coarse-grained model would also predict a universal dependence of on the dimensionless parameter . The observed non-universality could, however, be the result of corrections arising from contribution to that are second order in , which are neglected in our estimate of , and/or the result of end effects. Both of these effects are neglected in our analysis, and both are expected to yield model dependent corrections to that are proportional to . It appears to us that the dominant contribution to the deviation for this data may well be a nonuniversal correction proportional to (which appears to provide a better description of the dependence, as originally suggested by Deutsch and Binder) rather than a universal correction proportional to .

### v.3 Single-Chain Correlations

In Fig.9 we compare one-loop predictions and simulation
results of Beckrich et al. ^{11} for the deviation
of the correlation function for a monodisperse melt
from the Debye function. Simulation results for 5 different chain
lengths were taken from Ref. ^{11}.
As noted by Beckrich et al., the collapse of data for from
different chain lengths when is
plotted vs. already shows that the deviation is proportional
to , as predicted by the one-loop theory.
Simulations and predictions agree quite well, though there do appear
to be some small systematic discrepancies in the high regime.
The agreement is somewhat better than that obtained by the Strasbourg
group using a more approximate treatment of a monodisperse melt.
A similar level of agreement was obtained previously by that group
in comparisons of the one-loop theory to simulations of polydisperse
equilibrium polymers, for which they obtained analytic predictions
for .
^{20, 12, 21, 11, 27}.

Müller has examined the dimension of homopolymer chains of various
lengths in binary blends with varying compositions and interactions
^{28}, using Monte Carlo simulations of the BFM with the
same parameters as the one used to study composition fluctuations
^{25}. In Fig.10 we re-plotted
times his results for the fractional deviation of the
mean-squared end-to-end vector for component over a range of
values of for a fixed value ,
corresponding to , for 6 different chain lengths,
adapted from Fig.5b in ^{28}. Here represents
the deviation of in a mixture from the ideal Gaussian chain value
calculated using the statistical segment length of a sufficiently long
chain. The original results were presented as the difference between
the value in the blend and that in a melt of chains of the same length.
We have switched the normalization by assuming that the radius of gyration
measured in the melt is smaller compared to the ideal value by the amount
.
The fact that results for different chain lengths
collapse confirms by itself that these deviations are proportional to
. The method of obtaining the theoretical curve in Fig.10
is sketched in appendix B. The qualitative trend of the
simulation results were relatively well described with the theory.

## Vi Thermodynamic Integration of the RPA

The one-loop approximation has most often been presented within
the context of the Edwards formalism, as a Gaussian approximation
for the functional integral representation for the partition
function. Here, we discuss an alternative derivation that makes
no reference to the Edwards formalism, in which the one-loop free
energy is obtained by thermodynamic integration with respect to
the strength of either the pair potential or the bare interaction
parameter, while using the RPA approximation for .
The derivation is closely analogous to the usual derivation
of the Debye-Hückel approximation for excess free energy
of an electrolyte solutions, which is obtained by combining
an expression for the free energy as an integral with respect
to a ”charging” parameter with the Debye-Hückel linear
self-consistent field approximation for the correlation function
(see, e.g., McQuarrie ^{29}, pgs. 328-339).

### vi.1 Variation of the Pair Potential

Our analysis starts with the following exact expression for the dependence of the free energy upon the pair potential in any model with a pair potential: The perturbation in the free energy induced by an infinitesimal perturbation in the pair potential is given by an integral

(53) |

The corresponding perturbation of the correlation free energy density may expressed as an integral

(54) |

or as a Fourier integral

(55) |

where ,

By explicitly evaluating the corresponding variation of Eq. (10) for the one-loop correlation free energy for a compressible liquid, we find that, in this approximation,

(56) |

where

(57) |

is the RPA (or mean-field) approximation for in a compressible liquid. We thereby find that the one-loop approximation for the free energy density is thermodynamically consistent with an RPA approximation for .

It follows that the one-loop free energy may be obtained by a thermodynamic integration similar to that used to obtain the free energy in the Debye-Hückel theory of electrolytes: Consider a sequence of systems with a reduced pair potential , with a charging parameter in the range . The one-loop correlation energy (which vanishes in the ideal gas state) may be expressed as an integral