In [18]:

## SAXS: Bimodal silica particles¶

In this notebook we demonstrate modelling of small angle scattering data from bimodal silica particles. More details as well as experimental data you can find in the original paper from Applied Crystallography Journal

http://journals.iucr.org/j/issues/2015/05/00/vg5026/

This notebook is a demonstration of the core functionality of ESCAPE package.

The investigated population of silica nanoparticles consists of particles of two sizes: small particles of type "A" and large particles of type "B", each with its own Gaussian distribution of sizes with characteristic mean radii ${R}_{A}$, ${R}_{B}$ and their distribution widths ${\sigma }_{A}$ and ${\sigma }_{B}$.

The resulted intensity is calculated as a sum of the scattering intensity contributed by each scatterer, i.e. each assembly of particles as following

${I}_{mod}\left(q\right)=\sum _{c}^{N}\left[{\int }_{a}^{b}P\left(q,x\right)f\left(x\right)dx\right]S\left(q\right)$

Each contribution $c$ consists of a form-factor $P\left(q,x\right)$ determining the shape of a scatterer, $f\left(x\right)$ - distribution of certain parameter of the form-factor in the case of disperse aspects. $S\left(q\right)$ is a structure factor for each contribution, which is not relevant for our case and is equal unity.

We start implementation of our model with definition of model parameters. In the comments below we give a short description of each parameter. For the details please have a look on the original publication.

In [19]:

Now we can implement form-factors for our model. These are standard form-factors for homegeneous particles with spherical shape.

In [20]:

${P}_{A}$ and ${P}_{B}$ are the form-factors for the corresponding scatterers without taking into account the size distribution. Since ESCAPE operates mostly with simple functor-type objects, the user has the ability to look at and analyze the behavior of each object, before the final simulation. For example, to visually estimate the contribution of any part of the model or locate a calculation error.

In [21]:
8.522
1.017
37.65
1.017

After implementing of form-factors and checking their curves we are to finalize our model and to implement integration over size distribution. Here we had to reformulate the expressions, because in ESCAPE there is no Gaussian number-weighted size distribution. For that purposes we introduced expressions for ${C}_{A}$ and ${C}_{B}$ normalization constants which we will use as multipliers for the form-factor expressions. over size distribution.

In [22]:
8.522
2.007
1.030
1.017
37.65
8.296
6.511
0.000