Debye function


Fit a Debye function to the small angle neutron scattering data of ideal polymer chains.

  • Function
    \[I(q) = \frac {(b_H-b_D)^2}{V} N \phi(1-\phi) g_D(R_g,q) + \mathrm{Bckgrnd}\]
    • \(N=M_w/M_\mathrm{mono}\) is the degree of polymerization of the chain (\(M_w\) is a parameter of the experimental data)

    • \(\phi\) is the volume fraction of deuterated chains (read from the file)

    • \(g_D(R_g,q)\) is the Debye function, given by

    \[g_D(R_g,q) = \frac{2}{(q^2R_g^2)^2}\left( q^2R_g^2 + exp(-q^2R_g^2) -1 \right)\]
  • Parameters
    • Contrast: This sets the magnitude of the scattering and is equal to \((b_H-b_D)^2/V\) where \(b_{H/D}\) is the scattering cross-section of the hydrogenous/deuterated monomer and \(V\) is the monomer volume.

    • \(C_{gyr}\): This sets the scale of the radius of gyration of the chain. For a given molecular weight, the radius of gyration is \(R_g^2=C_{gyr}M_w\). For many polymers, this value is available in the literature, but small adjustments may still be necessary to optimize the agreement with the experimental data.

    • \(M_\mathrm{mono}\): The molecular weight of a single monomer (should be known from the chain chemistry).

    • Bckgrnd: This sets the level of the background scattering. It can, in principle, be computed from known incoherent scattering cross sections but, in practice, there are often many other unknown contributions and therefore fitting is necessary.

    • \(\lambda\): It applies a simple strain measure by shifting the radius of gyration by a constant factor for all \(q\) values, \(R_g\to \lambda R_g\) (the Stretched button must be checked). This can be used to compare the microscopic deformation with the effect of a fully affine bulk deformation or to fit to the low \(q\) data to produce an effective radius of gyration under flow. Compression perpendicular to the flow direction can be modelled by setting \(\lambda<1\).

    • \(\chi\): Parameter to model the effect of a weak interaction between the hydrogenous and the deuterated monomers on the scattering, modelled within the random phase approximation [3] (the Non-Ideal Mix button must be checked). The scattered intensity is calculated according to the function below, in which \(\chi\) is independent of \(M_w\) and \(\phi\) but is expected to change with temperature. Typically, the effect of \(\chi\) is small, but this depends upon the temperature, degree of polymerization and deuterated fraction. For deuterated/hydrogenated polystyrene \(\chi\approx 1.7\cdot 10^{-4}\) at 160 degrees C, and it is expected to be smaller with increasing temperature.

      \[I(q) = \frac {(b_H-b_D)^2}{V} \left( \frac{1}{N \phi(1-\phi) g_D(R_g,q)}-2\chi \right)^{-1} + \mathrm{Bckgrnd} \]


This theory accounts for the small angle neutron scattering from melts of ideal Gaussian chains. It is intended to simultaneously fit quiescent scattering data for monodisperse linear polymers, for several molecular weights. The theory uses a set of molecular parameters which can be obtained from the literature, or can be fitted to data if unavailable elsewhere. For full details of the background theory see [1]; for an example of its application see [2]; for details about the effect of non-ideal mixing, see [3].