SasView: Ticket #217: Verify Correctness of Shulz Zimm Distribution
http://trac.sasview.org/ticket/217
<p>
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.
</p>
en-usTrac 1.0.1butlerTue, 01 Apr 2014 15:01:30 GMTsummary changed
http://trac.sasview.org/ticket/217#comment:1
http://trac.sasview.org/ticket/217#comment:1
<ul>
<li><strong>summary</strong>
changed from <em>Finish code reorg</em> to <em>Verify Correctness of Shulz Zimm Distribution</em>
</li>
</ul>
TicketbutlerFri, 04 Apr 2014 09:59:05 GMTowner, status, milestone changed
http://trac.sasview.org/ticket/217#comment:2
http://trac.sasview.org/ticket/217#comment:2
<ul>
<li><strong>owner</strong>
changed from <em>Gonzalez</em> to <em>butler</em>
</li>
<li><strong>status</strong>
changed from <em>new</em> to <em>assigned</em>
</li>
<li><strong>milestone</strong>
changed from <em>SasView 3.0.0</em> to <em>Admin Tasks</em>
</li>
</ul>
<p>
After lengthy discussions and looking at papers and doing calculations it is clear that:
<a class="missing wiki">SasView?</a> 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.
</p>
<p>
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.
</p>
<p>
Fundamentally this means that there are TWO different distributions which are being called the Schulz distribution. These are unfortunately related through the mass.
</p>
<p>
<a class="missing wiki">SasView?</a>, 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.
</p>
<p>
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 <a class="missing wiki">SasView?</a> version of the distribution, or something else (such as providing different moments of the distribution, or plots of the distribution etc)
</p>
TicketbutlerFri, 04 Apr 2014 10:00:17 GMTowner changed
http://trac.sasview.org/ticket/217#comment:3
http://trac.sasview.org/ticket/217#comment:3
<ul>
<li><strong>owner</strong>
changed from <em>butler</em> to <em>gonzales/butler/ajj/heenan</em>
</li>
</ul>
TicketbutlerFri, 30 Jan 2015 01:47:17 GMTworkpackage set
http://trac.sasview.org/ticket/217#comment:4
http://trac.sasview.org/ticket/217#comment:4
<ul>
<li><strong>workpackage</strong>
set to <em>SasView Admin</em>
</li>
</ul>
TicketbutlerSun, 01 Feb 2015 23:43:42 GMTpriority changed
http://trac.sasview.org/ticket/217#comment:5
http://trac.sasview.org/ticket/217#comment:5
<ul>
<li><strong>priority</strong>
changed from <em>major</em> to <em>trivial</em>
</li>
</ul>
TicketgonzalezmFri, 27 Feb 2015 11:46:56 GMTstatus changed; resolution set
http://trac.sasview.org/ticket/217#comment:6
http://trac.sasview.org/ticket/217#comment:6
<ul>
<li><strong>status</strong>
changed from <em>assigned</em> to <em>closed</em>
</li>
<li><strong>resolution</strong>
set to <em>invalid</em>
</li>
</ul>
<p>
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.
</p>
Ticketsmk78Wed, 04 Mar 2015 13:04:13 GMT
http://trac.sasview.org/ticket/217#comment:7
http://trac.sasview.org/ticket/217#comment:7
<p>
Web FAQ updated accordingly.
</p>
Ticket