#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 11 years ago by butler
- Summary changed from Finish code reorg to Verify Correctness of Shulz Zimm Distribution
comment:2 Changed 11 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
comment:3 Changed 11 years ago by butler
- Owner changed from butler to gonzales/butler/ajj/heenan
comment:4 Changed 10 years ago by butler
- Work Package set to SasView Admin
comment:5 Changed 10 years ago by butler
- Priority changed from major to trivial
comment:6 Changed 10 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 10 years ago by smk78
Web FAQ updated accordingly.
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 Shulz 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 Shulz 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 Aragon and Pecora's number average formula but call it 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)