source: sasmodels/doc/guide/sesans/sans_to_sesans.rst @ ff31782

core_shell_microgelsmagnetic_modelticket-1257-vesicle-productticket_1156ticket_1265_superballticket_822_more_unit_tests
Last change on this file since ff31782 was f0fc507, checked in by GitHub <noreply@…>, 7 years ago

Clarification of the Hankel transform

I've tried to tidy up the documentation and get it more in line with what's going on with in the actual code.

  • Property mode set to 100644
File size: 2.0 KB

SANS to SESANS conversion

?
.. currentmodule:: sasmodels

The conversion from SANS into SESANS in absolute units is a simple Hankel transformation when all the small-angle scattered neutrons are detected. First we calculate the Hankel transform including the absolute intensities by

G(δ) = 2π0J0(Qδ)(dΣ)/(dΩ)(Q)QdQ, 

in which J0 is the zeroth order Bessel function, δ the spin-echo length, Q the wave vector transfer and (dΣ)/(dΩ)(Q) the scattering cross section in absolute units.

Out of necessity, a 1-dimensional numerical integral approximates the exact Hankel transform. The upper bound of the numerical integral is Qmax , which is calculated from the wavelength and the instrument's maximum acceptance angle, both of which are included in the file. While the true Hankel transform has a lower bound of zero, most scattering models are undefined at :math: Q=0, so the integral requires an effective lower bound. The lower bound of the integral is Qmin = 0.1 × 2π ⁄ Rmax , in which Rmax is the maximum length scale probed by the instrument multiplied by the number of data points. This lower bound is the minimum expected Q value for the given length scale multiplied by a constant. The constant, 0.1, was chosen empirically by integrating multiple curves and finding where the value at which the integral was stable. A constant value of 0.3 gave numerical stability to the integral, so a factor of three safety margin was included to give the final value of 0.1.

From the equation above we can calculate the polarisation that we measure in a SESANS experiment:

P(δ) = et(λ)/(2π)2(G(δ) − G(0)), 

in which t is the thickness of the sample and λ is the wavelength of the neutrons.

Note: See TracBrowser for help on using the repository browser.