next up previous contents index
Next: Appendix A Up: Integral equations for polymer Previous: Molecular liquids

Polymer RISM

We now apply the RISM  formalism to polymeric liquids . For simplicity we restrict the presentation to polyethylene , each CH2 will be viewed as one ''atom''. It is clear that if boundary effects may be neglected, all monomers along the chain are equivalent, i.e. that $G_{\alpha \beta }(r)$ and all other functions do not depend on $\alpha $ or $\beta $. Remember that $\alpha $and $\beta $ each run from 0 to N, where (N+1) is the number of monomers per chain. The monomer-monomer radial distribution function then is

 \begin{displaymath}G(r)=
\frac{1}{\rho _{m}}\sum_{\beta }G_{\alpha \beta }(r)\rho _{\beta }=G_{\alpha
\beta }(r)
\end{displaymath} (2.35)

where $\rho _{m}=(N+1)\rho _{\beta }$ is the monomer density. Introducing Eq. (2.34) into Eq. (2.35) we get
 
G(r) = $\displaystyle \frac{1}{\rho _{m}}\{\omega (r)-\delta (r)\}+g(r)$ (2.36)
$\displaystyle \omega (r)$ = $\displaystyle \sum_{\beta }\Omega _{\alpha \beta }(r)$ (2.37)

where $g(r)=g_{\alpha \beta }(r)$. Using $h(r)=h_{\alpha \beta }(r)$ and $
c(r)=c_{\alpha \beta }(r)$ in Eq. (2.31) we derive

\begin{displaymath}\hat{h}(k)=\hat{\omega}(k)\hat{c}(k)\hat{\omega}(k)+\rho _{m}\hat{\omega}(k)
\hat{c}(k)\hat{h}(k)
\end{displaymath} (2.38)

$\hat{\omega}(k)$ may be calculated using the methods of Chapter 1. Some details are given in Appendix A.

We now apply the method to polyethylene, which will be modelled by a RIS model with characteristics given in section 2.6. The interaction between two monomers on different chains is modelled by a hard sphere interaction with diameter $\sigma = 3.9 \mbox{\AA}$. Moreover N=6416 and $\rho_m = 33.58
\mathrm{{nm^{-3}}}$. Calculations were done with PY closure.

In Fig. (2.3) is given the intermolecular distribution function. We clearly see what is called the ''correlation hole''  by de Gennes ; the intermolecular distribution gradually rises to its limit value 1, with only little structure. The scattering function is given in Fig. (2.4), and is in perfect agreement with X-ray experiments, except at very small values of k.

We give two more results, from MD simulations  this time. In Fig. (2.5 ) is given the intramolecular oxygen-oxygen correlation in PEO . The simulation was done with GROMOS united atom potentials. The box consisted of two chains of 800 monomers each; $\rho
_{m}=14.45\mathrm{{nm^{-3}}}$ and T=400K. The corresponding intermolecular correlation function is given in Fig. (2.6).


  
Figure 2.3: Intermolecular distribution function.
\scalebox{0.5}{\includegraphics[14,85][580,385]{fig2_3.eps}}


  
Figure 2.4: Scattering function.
\scalebox{0.5}{\includegraphics[14,85][580,385]{fig2_4.eps}}


  
Figure 2.5: Intramolecular oxygen-oxygen correlation in polyethyleneoxide.
\scalebox{0.5}{\includegraphics[14,450][580,760]{fig2_5.eps}}


  
Figure 2.6: Intermolecular correlation function in polyethyleneoxide.
\scalebox{0.5}{\includegraphics[14,450][580,760]{fig2_6.eps}}


next up previous contents index
Next: Appendix A Up: Integral equations for polymer Previous: Molecular liquids
W.J. Briels