Self-correlation: a unique measure of SAXS information content?

The evolution of the reduced Chi-squared as a function of the shift along q.
The evolution of the reduced Chi-squared as a function of the shift along q.

A funny point came up last week when we were trying to assess the Q-uncertainty inherent in our measurement. We were trying a method inspired by one from Christian Gollwitzer (part 3.3.2, paper here). What we did is the following:

 

A pattern with sharp features compared to itself shifted by 0.1 1/nm.
Figure 1: A pattern with sharp features compared to itself shifted by 0.1 1/nm.

We take a given measured pattern (absolute scale, complete with uncertainties), and shift it by a tiny bit in q. We then compare the fit of the shifted pattern with its original (including scaling the intensity), and calculate the goodness-of-fit parameter \chi^2_r (Figure 1). In order to calculate this self-similar \chi^2_r, we do need to perform a (linear) interpolation on the intensity and the uncertainties.

A smoothly decaying pattern compared to itself shifted by 0.1 1/nm.
Figure 2: A smoothly decaying pattern compared to itself shifted by 0.1 1/nm.

With this method, we can calculate the dependency of \chi^2_r for a given shift in q. As you may remember, \chi^2_r \leq 1 when the two patterns overlap on average to within the uncertainty of the data. The q shift at which we reach \chi^2_r \approx 1 is our maximum allowed shift q_{\mathrm{maxshift}} (Figure 4).

A mostly flat pattern compared to itself shifted by 0.1 1/nm.
Figure 3: A mostly flat pattern compared to itself shifted by 0.1 1/nm.

This maximum allowed shift will vary depending on the shape in the datasets: the allowed q_{\mathrm{maxshift}} is small (0.005) for patterns with sharp features (Figure 1), broader (0.04) for smooth patterns (Figure 2), and very large (0.1) for flat profiles (Figure 3). These values should be compared with the smallest q value in our data: q_{min} = 0.1. So what we seem to have, then, is some sort of measure for the q-dependent information content in the scattering pattern.

The evolution of the reduced Chi-squared as a function of the shift along q.
Figure 4: The evolution of the reduced Chi-squared as a function of the shift along q.

This could maybe be used to estimate the uncertainty in our q-vector, provided we take a pattern with sharp features, and preferably know the size of the object that it represents. Perhaps something silver behenate would do. In any case, we know that the uncertainty determined in this manner lies well below the width the datapoints themselves describe.

For now, at least, we have an interesting measure and some food for thought.

1 Comment

  1. The uncertainty in the q-vector is an important measure. But a priori it is not clear how it affects the parameter of our, let’s say, fitted model. For large spheres of low dispersity, a shift of the q-vector has a dramatic effect. But if you measure the absolute intensity of water, i.e. its isothermal compressibiliy, a huge shift of the q-vector has no effect. Therefore, I think that it is important to know what you are looking at to estimate the uncertainty of your parameteres.

2 Trackbacks / Pingbacks

  1. Page not found – Looking At Nothing
  2. A complete q-uncertainty story – Looking At Nothing

Leave a Reply

Your email address will not be published.


*


*

This site uses Akismet to reduce spam. Learn how your comment data is processed.