Opened 4 years ago

Closed 4 years ago

Last modified 3 years ago

#1105 closed defect (fixed)

mass_surface_fractal model documentation, unit tests, defaults are wrong

Reported by: richardh Owned by: richardh
Priority: critical Milestone: SasView 4.2.0
Component: sasmodels Keywords:
Cc: Work Package: SasView Bug Fixing


With thanks to Greg Smith for querying this.

The equations and code in mass_surface_fractal are I think correct, however the documentation, default parameters and unit test examples in sasview are wrong!

The “primary particle” radius should be the smaller one (not the larger one) which is envisaged as having a surface fractal (i.e. smooth or rough surfaces) and then sticking together into a larger cluster radius via a mass fractal packing. (The authors also note that polydispersity in the primary particle size might produce similar effects to the surface fractal at high Q, so results might be ambiguous. Though for their fumed silicas microscopy tends to show that the primary particles are pretty monodisperse.)

The problem seems to have arisen in that the Schmidt review paper referenced in the documentation attempted to reformulate the equations and variable names from the original Hurd, Schaefer & Martin paper but in doing so managed to get the description wrong, the equations incomplete, and swap the primary and cluster radii around. The sasview documentation alas followed this erroneous description!

Change History (3)

comment:1 Changed 4 years ago by smk78

Greg Smith adds:
The Hurd paper was the one I was using for inspiration (in general actually, as we are looking at fumed silica scatters), and I have just cross-referenced it against the SasView documentation. The equations seem correct in relation to the Hurd paper, it just seems the default parameters that are transposed. (If indeed in the code, rg_cluster = R_g and rg_primary = r_g, as it seems you suggest as well.) The default parameters appear to be HS-5 Cab-O-Sil from Table 1 in Hurd (Phys Rev A), just with the primary and cluster rg the other way around.

comment:2 Changed 4 years ago by GitHub <noreply@…>

  • Resolution set to fixed
  • Status changed from new to closed

In 02490b34b32c2e133225c89d566c9481f4c53f83/sasmodels:

Merge pull request #68 from SasView?/ticket-1105_mass_surface_fractal

ticket 1105 edit docs and default params for mass_surface_fractal

closes #1105

comment:3 Changed 3 years ago by butler

  • Milestone changed from SasView 4.3.0 to SasView 4.2.0
Note: See TracTickets for help on using tickets.