Opened 3 years ago

Closed 3 years ago

Last modified 3 years ago

#217 closed task (invalid)

Verify Correctness of Shulz Zimm Distribution

Reported by: butler Owned by: gonzales/butler/ajj/heenan
Priority: trivial Milestone: Admin Tasks
Component: SasView Keywords:
Cc: Work Package: SasView Admin

Description

Klaus Huber and Ralph Schweins believe we are effectively using the weight average distribution rather than the number average in the case of the Shulz-Zimm distribution. Need to verify. If wrong should fix. If correct need to document.

Change History (7)

comment:1 Changed 3 years ago by butler

  • Summary changed from Finish code reorg to Verify Correctness of Shulz Zimm Distribution

comment:2 Changed 3 years ago by butler

  • Milestone changed from SasView 3.0.0 to Admin Tasks
  • Owner changed from Gonzalez to butler
  • Status changed from new to assigned

After lengthy discussions and looking at papers and doing calculations it is clear that:
SasView? assumes the underlying distribution is a number distribution and thereby the reported mean from that distribution is a number average parameter. It properly weights that distribution by volume when fitting to the scattering data as is required to fit scattering data.

However, from the perspective of the Schulz distribution, it is clear that the litterature has conflicting views about what the correct form of the number average Schulz should be.

Fundamentally this means that there are TWO different distributions which are being called the Schulz distribution. These are unfortunately related through the mass.

SasView?, like FISH and NIST IGOR macros, take their distribution from the M. Kotlarchyk and S-H. Chen (J. Chem. Phys. 79 (1983) 2461-2469) paper which specifically gives the distribution as being by number. This is the same number Schulz distrbition given by Aragon and Pecora in their 1976 J Chem Phys paperaper as well as in the 1992 J. Chem Phys paper of Bartlett and Ottewill, among others. On the other hand, the early papers on this distribution were specifically applied to degrees of polymerization but were clear that they viewed the equivalent formula to be a weight average degree of polymerization. It seems that at least Welch and Bloomfield in their 1973 J. Pol. Sci. paper and Klaus Huber 2012 now use the same equation as Aragon and Pecora's number average formula but define it as a weight average formula and compute a new number avarage formula therefrom.

At the end these are just numerical distributions and can both be used but the question is whether to add another distribution with the same name (which means making sure we distinguish somehow), leave as is based on the fact the the preponderance of usage in scattering seems to be the current SasView? version of the distribution, or something else (such as providing different moments of the distribution, or plots of the distribution etc)

Last edited 3 years ago by butler (previous) (diff)

comment:3 Changed 3 years ago by butler

  • Owner changed from butler to gonzales/butler/ajj/heenan

comment:4 Changed 3 years ago by butler

  • Work Package set to SasView Admin

comment:5 Changed 3 years ago by butler

  • Priority changed from major to trivial

comment:6 Changed 3 years ago by gonzalezm

  • Resolution set to invalid
  • Status changed from assigned to closed

Finally we have concluded that there is no problem with the Schulz distribution. Using the current definition for the z parameter describing the width of the distribution, it results that the weight distribution using the weight average has exactly the same mathematical form than the number distribution using the number average. Therefore there is no need to add an alternative Schulz distribution. The misunderstanding came from the fact that z is defined in some paper in an alternative way, so that z' = z+1, and in that case the weight and number distributions do not have the same mathematical form. Klaus Huber has checked this and agrees.

comment:7 Changed 3 years ago by smk78

Web FAQ updated accordingly.

Note: See TracTickets for help on using tickets.