Self-Consistent Field Theory and Its Applications by M. W. Matsen
Self-Consistent Field Theory and Its Applications by M. W. Matsen
Self-Consistent Field Theory and Its Applications by M. W. Matsen
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
74 1 <strong>Self</strong>-consistent field theory <strong>and</strong> its applications<br />
where φ A (r) is the concentration from n A noninteracting A-type polymers subjected to the<br />
complex field, W A (r) =W + (r)+W − (r). In this case, φ A (r) is a complex quantity with no<br />
physical significant. The same procedure can be adapted equally well to the partition functions<br />
of other systems such as the diblock copolymer melts considered in Section 1.8.1.<br />
Since the integr<strong>and</strong> of Z in Eq. (1.316) involves a Boltzmann factor, exp(−F/k B T ),<br />
st<strong>and</strong>ard simulation techniques of statistical mechanics can be applied, where the coordinates<br />
are now the fluctuating fields, W − (r) <strong>and</strong> W + (r). Ganesan <strong>and</strong> Fredrickson (2001) have used<br />
Langevin dynamics while Düchs et al. (2003) have formulated a Monte Carlo algorithm for<br />
generating a sequence 〉 of configurations with the appropriate Boltzmann weights. The ensemble<br />
average,<br />
〈ˆφA (r) , can then be approximated <strong>by</strong> averaging φ A (r) over the configurations<br />
generated <strong>by</strong> the simulation. Even though φ A (r) is generally complex, the imaginary part will<br />
average to zero. The need to consider two separate fluctuating fields is still computationally<br />
dem<strong>and</strong>ing, but in practice the functional integral over the field, W + (r), can be performed<br />
with the saddle-point approximation, which implies that as W − (r) fluctuates, W + (r) is continuously<br />
adjusted so that φ A (r) +φ B (r) =1(Düchs et al. 2003). Even with the problem<br />
reduced to a single fluctuating field, it is still an immense challenge to perform simulations<br />
in full three-dimensional space. However, with the rapid improvement in computational performance<br />
<strong>and</strong> the development of improved algorithms, it is just a matter of time before these<br />
field-theoretic simulations become viable.<br />
1.11 Appendix: The Calculus of Functionals<br />
SCFT is a continuum field theory that relies heavily on the calculus of functionals. For those<br />
not particularly familiar with this more obscure variant of calculus, this section reviews the<br />
essential aspects required for a comprehensive underst<strong>and</strong>ing of SCFT. In short, the term<br />
functional refers to a function of a function. A specific example of a functional is<br />
F[f] ≡<br />
∫ 1<br />
0<br />
dx exp(f(x)) (1.319)<br />
This definition provides a well-defined rule for producing a scalar number from the input<br />
function, f(x). For example, F[f] returns the value, e, for the input, f(x) =1, <strong>and</strong> it returns<br />
(e − 1) for f(x) =x. Note that we follow a convention, where the arguments of a functional<br />
are enclosed in square brackets, rather than the round brackets generally used for ordinary<br />
functions.<br />
The calculus of functionals is, in fact, very much like multi-variable calculus, where the<br />
value of f(x) for each x acts as a separate variable. Take the case where 0 ≤ x ≤ 1,<br />
<strong>and</strong> imagine that the x coordinate is divided into a discrete array of equally spaced points,<br />
x m ≡ m/M for m =0, 1, 2, ..., M. Provided that M is large, the function, f(x), is well<br />
represented <strong>by</strong> the set of values, {f 0 ,f 1 , ..., f M }, where f m ≡ f(x m ). In turn, the functional<br />
in Eq. (1.319) is well approximated <strong>by</strong> the multi-variable function,<br />
F({f m }) ≡ exp(f 0)<br />
2M<br />
+ 1 M<br />
M−1<br />
∑<br />
m=1<br />
exp(f m )+ exp(f M)<br />
2M<br />
(1.320)