{"id":846,"date":"2013-09-09T12:00:39","date_gmt":"2013-09-09T03:00:39","guid":{"rendered":"http:\/\/www.lookingatnothing.com\/?p=846"},"modified":"2013-12-14T15:37:52","modified_gmt":"2013-12-14T06:37:52","slug":"how-to-fit-a-scattering-pattern-part-1-the-initial-look","status":"publish","type":"post","link":"https:\/\/lookingatnothing.com\/index.php\/archives\/846","title":{"rendered":"How to fit a scattering pattern part 1: the initial look."},"content":{"rendered":"<p><em>By the way, my topical review paper on SAXS data collection and correction has been published and is available open access <a href=\"http:\/\/iopscience.iop.org\/0953-8984\/25\/38\/383201\/\">here<\/a>!<\/em><\/p>\n<p>&nbsp;<\/p>\n<p>Recently, some good colleagues (who are not familiar with scattering) have started asking questions on how to go about fitting a scattering pattern. This was a very good opportunity to think about the process from a layman perspective. How do we get from scattering pattern to morphological information in a straightforward way?<!--more--><\/p>\n<p>For me, this typically takes about three steps: 1) the initial look, 2) some simple fits and 3) the final fit. There is a possible fourth step: advanced fits, but this is not a necessity in some cases. This post will focus on the first step: the initial look.<\/p>\n<p>So imagine you have collected a scattering pattern and subjected it to the applicable corrections to get a background-subtracted, <a href=\"http:\/\/iopscience.iop.org\/0953-8984\/25\/38\/383201\/\">corrected<\/a> scattering pattern. This sets the stage for step 1. In step 1, we check if there is any suitable information in the scattering pattern at all, and we collect information from other techniques.<\/p>\n<figure id=\"attachment_848\" aria-describedby=\"caption-attachment-848\" style=\"width: 300px\" class=\"wp-caption alignright\"><a href=\"http:\/\/www.lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Cyl_intensity.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-medium wp-image-848\" alt=\"Simulated scattering pattern of polydisperse cylinders (with a background), where typical SAXS regions have been identified.\" src=\"http:\/\/www.lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Cyl_intensity-300x188.png\" width=\"300\" height=\"188\" srcset=\"https:\/\/lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Cyl_intensity-300x188.png 300w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Cyl_intensity-1024x644.png 1024w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Cyl_intensity-150x94.png 150w\" sizes=\"auto, (max-width: 300px) 100vw, 300px\" \/><\/a><figcaption id=\"caption-attachment-848\" class=\"wp-caption-text\"><strong>Figure 1:<\/strong> Simulated scattering pattern of polydisperse cylinders (with a background), where typical SAXS regions have been identified. Click to embiggen.<\/figcaption><\/figure>\n<p>What I&#8217;d do is to plot the data on a double-logarithmic plot, in intensity versus scattering vector: <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=%5Clog%28I%29&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"&#92;log(I)\" class=\"latex\" \/> versus <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=%5Clog%28q%29&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"&#92;log(q)\" class=\"latex\" \/>. Bonus points if you can plot error bars on the intensity. You can then try to identify different features in your scattering pattern. At least you should be able to see an asymptotic behaviour towards a constant (flat) background. Furthermore, you can likely identify a straight region of a constant slope which is called the Porod region. A simple example exhibiting many features is shown in Figure 1. We will discuss the features starting from those at highest <em>q<\/em>-values.<\/p>\n<p>The constant background contribution can arise from a plethora of causes. It can come from very small-scale inhomogeneities in one or both of the phases in your sample, or remains from insufficient background subtraction or insufficient correction for detector darkcurrent levels [1,2]. There is a modicum of information in the background level, but usually it is partially or wholly artificial.<\/p>\n<p>The Porod region is another region of limited use. The scattering intensity in this region should be proportional to a <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=q%5E%7B-4%7D&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"q^{-4}\" class=\"latex\" \/> slope, but can also assume values lower or higher than that [1,2,3]. The <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=q%5E%7B-4%7D&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"q^{-4}\" class=\"latex\" \/> scattering arises from large scatterers beyond their characteristic scattering region. In other words, it indicates that there are scattering objects, but it does not tell you what size they are. We need to look at lower <em>q<\/em>-values to get <em>that<\/em> information.<\/p>\n<p>Porod slopes smaller than <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=q%5E%7B-4%7D&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"q^{-4}\" class=\"latex\" \/> could indicate the presence of a fractal-like structure in your sample: a self-similar structure over a certain size range. This could be a rough surface, but can also arise in very polydisperse mixtures [4,5].<\/p>\n<p>At <em>q<\/em>-values lower than the Porod slope, there can be a region which contains the most information in the scattering pattern which I shall call the Debye region (after the early description of this region by Debye, B\u00fcche, Anderson and Brumberger [6,7]). This region typically has a slope between <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=q%5E%7B-1%7D&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"q^{-1}\" class=\"latex\" \/> and <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=q%5E%7B-3%7D&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"q^{-3}\" class=\"latex\" \/>, and &#8212;when analysed&#8212; provides the size distribution of the structure.<\/p>\n<p>Some say that a particular slope identifies the shape of the scattering entities (discs, cylinders or spheres), but this is not completely unambiguous. As will be shown later, the scattering pattern from polydisperse cylinders can be perfectly described using spheres as well (though perhaps the reverse does not hold, I did not test that yet\u2026). All that can be said, then, is that this region is the one carrying the information.<\/p>\n<p>Before that Debye region lies the Guinier region: a flat region independent of q. This contains information on the volume-square weighted size of the scatterers, and also can be used to find the total volume fraction of scatterers in the sample. As it typically lies at very low q values (and so very close to the beamstop), it is also most easily affected by parasitic scattering and other unrelated scattering effects. In short: This region has some information, but is prone to fouling.<\/p>\n<p>Figure 1 shows the best case scenario for data: that in which all regions are clearly visible and distinct from one another. I will provide a few more examples to indicate some alternatives you might find in your scattering patterns.<\/p>\n<figure id=\"attachment_849\" aria-describedby=\"caption-attachment-849\" style=\"width: 300px\" class=\"wp-caption alignright\"><a href=\"http:\/\/www.lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Disc_gauss_s10_m30_bgnd1e1_q4slope_sasfit.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-medium wp-image-849\" alt=\"Simulated scattering of polydisperse discs in the presence of large scatterers and a flat background contribution.\" src=\"http:\/\/www.lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Disc_gauss_s10_m30_bgnd1e1_q4slope_sasfit-300x229.png\" width=\"300\" height=\"229\" srcset=\"https:\/\/lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Disc_gauss_s10_m30_bgnd1e1_q4slope_sasfit-300x229.png 300w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Disc_gauss_s10_m30_bgnd1e1_q4slope_sasfit-1024x784.png 1024w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Disc_gauss_s10_m30_bgnd1e1_q4slope_sasfit-150x114.png 150w\" sizes=\"auto, (max-width: 300px) 100vw, 300px\" \/><\/a><figcaption id=\"caption-attachment-849\" class=\"wp-caption-text\"><strong>Figure 2:<\/strong> Simulated scattering of polydisperse discs in the presence of large scatterers and a flat background contribution.<\/figcaption><\/figure>\n<p>Figure 2 shows another situation where there is a certain amount and distribution of discs (the hump in the middle), and a background contribution. I also simulated the presence of some larger structures (i.e. aggregates or simply large scattering objects), of which only the <img decoding=\"async\" src=\"https:\/\/s0.wp.com\/latex.php?latex=q%5E%7B-4%7D&#038;bg=ffffff&#038;fg=000&#038;s=0&#038;c=20201002\" alt=\"q^{-4}\" class=\"latex\" \/> slope is visible. The presence of this slope makes it more difficult to distinguish the individual regions, but the hump is the area of interest. It is that area which contains the information on the discs.<\/p>\n<figure id=\"attachment_850\" aria-describedby=\"caption-attachment-850\" style=\"width: 300px\" class=\"wp-caption alignright\"><a href=\"http:\/\/www.lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Sphere_gauss_s10_m30_bgnd1e2_q4slope_struct_0p22_sasfit.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-medium wp-image-850\" alt=\"Simulated scattering of densely packed polydisperse spheres in the presence of large scatterers and a flat background contribution\" src=\"http:\/\/www.lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Sphere_gauss_s10_m30_bgnd1e2_q4slope_struct_0p22_sasfit-300x230.png\" width=\"300\" height=\"230\" srcset=\"https:\/\/lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Sphere_gauss_s10_m30_bgnd1e2_q4slope_struct_0p22_sasfit-300x230.png 300w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Sphere_gauss_s10_m30_bgnd1e2_q4slope_struct_0p22_sasfit-1024x786.png 1024w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Sphere_gauss_s10_m30_bgnd1e2_q4slope_struct_0p22_sasfit-150x115.png 150w\" sizes=\"auto, (max-width: 300px) 100vw, 300px\" \/><\/a><figcaption id=\"caption-attachment-850\" class=\"wp-caption-text\"><strong>Figure 3:<\/strong> Simulated scattering of densely packed polydisperse spheres in the presence of large scatterers and a flat background contribution<\/figcaption><\/figure>\n<p>Figure 3 shows an interesting case of scatterers (spheres in this case) which are quite densely packed. This dense packing causes inter particle interference which shows up as a depression in the scattering pattern. It may not necessarily be as clear as this, and its visibility can reduce due to a plethora of effects (polydispersity, background, etc). Structure factors can be a pain in the ass to fit. Diluting the sample to remove this effect is recommended if possible.<\/p>\n<p>A common problem in data correction is that of oversubtraction of the background measurement (because people tend not to correct for transmission). Fortunately, it is easy to recognise this mistake, as shown in Figure 4: there is a marked drop to negative values of the scattering signal at higher q-values.<\/p>\n<figure id=\"attachment_854\" aria-describedby=\"caption-attachment-854\" style=\"width: 300px\" class=\"wp-caption alignright\"><a href=\"http:\/\/www.lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Sphere_gauss_s20_m30_bgndm1e3_q4slope_0p22_sasfit.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-medium wp-image-854\" alt=\"Figure 4: An extreme case of oversubtraction of a background signal. The marked drop at high q is a clear indication.\" src=\"http:\/\/www.lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Sphere_gauss_s20_m30_bgndm1e3_q4slope_0p22_sasfit-300x224.png\" width=\"300\" height=\"224\" srcset=\"https:\/\/lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Sphere_gauss_s20_m30_bgndm1e3_q4slope_0p22_sasfit-300x224.png 300w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Sphere_gauss_s20_m30_bgndm1e3_q4slope_0p22_sasfit-1024x765.png 1024w, https:\/\/lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Sphere_gauss_s20_m30_bgndm1e3_q4slope_0p22_sasfit-150x112.png 150w\" sizes=\"auto, (max-width: 300px) 100vw, 300px\" \/><\/a><figcaption id=\"caption-attachment-854\" class=\"wp-caption-text\"><strong>Figure 4:<\/strong> An extreme case of oversubtraction of a background signal. The marked drop at high q is a clear indication.<\/figcaption><\/figure>\n<p>The simulations shown here are made using <a href=\"http:\/\/sourceforge.net\/projects\/sasfit\/\">SASfit<\/a> by Joachim Kohlbrecher and <a href=\"http:\/\/www.ingobressler.net\">Ingo Bressler<\/a>. I highly recommend downloading that software, to learn how to use it and to simulate a few scattering patterns to get a feel for what small-angle scattering can look like. Another good package is the <a href=\"http:\/\/usaxs.xray.aps.anl.gov\/staff\/ilavsky\/irena.html\">Irena<\/a> package by\u00a0<a href=\"http:\/\/usaxs.xray.aps.anl.gov\/staff\/ilavsky\/\">Jan Ilavsky<\/a>\u00a0for <a href=\"http:\/\/www.wavemetrics.com\/products\/igorpro\/igorpro.htm\">Igor Pro<\/a>.<\/p>\n<p>Next up in this series: part 2: basic fitting.<\/p>\n<p>[1] W. Ruland. Small-angle scattering of two-phase systems: determination and significance of systematic deviations from Porod\u2019s law. J. Appl. Cryst., 4:70\u201373, 1971.<\/p>\n<p>[2] J. T. Koberstein, B. Morra, and R. S. Stein. The determination of diffuse boundary thicknesses of polymers by small-angle x-ray scattering. J. Appl. Cryst., 13:34\u201345, 1980.<\/p>\n<p>[3] R. Sobry, J. Ledent, and F. Fontaine. Application of an extended porod law to the study of the ionic aggregates in telechelic ionomers. J. Appl. Cryst., 24:516\u2013525, 1991.<\/p>\n<p>[4] J. Teixeira. Small-angle scattering by fractal systems. J. Appl. Cryst., 21:781\u2013785, 1988.<\/p>\n<p>[5] G. Beaucage. Small-angle scattering from polymeric mass fractals of arbitrary mass-fractal dimension. J. Appl. Cryst., 29:134\u2013146, Jan 1996.<\/p>\n<p>[6] P. Debye and A. M. Bueche. Scattering by an inhomogeneous solid. J. Appl. Phys., 20:518\u2013525, 1949.<\/p>\n<p>[7] P. Debye, H. R. J. Anderson, and H. Brumberger. Scattering by an inhomogeneous solid. ii. the correlation function and its application. Journal of Applied Physics, 28:679\u2013683, 1957.<\/p>\n","protected":false},"excerpt":{"rendered":"<div class=\"mh-excerpt\"><p>By the way, my topical review paper on SAXS data collection and correction has been published and is available open access here! &nbsp; Recently, some <a class=\"mh-excerpt-more\" href=\"https:\/\/lookingatnothing.com\/index.php\/archives\/846\" title=\"How to fit a scattering pattern part 1: the initial look.\">[&#8230;]<\/a><\/p>\n<\/div>","protected":false},"author":1,"featured_media":848,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"jetpack_post_was_ever_published":false,"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_memberships_contains_paid_content":false,"activitypub_content_warning":"","activitypub_content_visibility":"","activitypub_max_image_attachments":4,"activitypub_interaction_policy_quote":"anyone","activitypub_status":"","footnotes":"","jetpack_publicize_message":"","jetpack_publicize_feature_enabled":true,"jetpack_social_post_already_shared":true,"jetpack_social_options":{"image_generator_settings":{"template":"highway","default_image_id":0,"font":"","enabled":false},"version":2}},"categories":[1,85,55],"tags":[18,25,179,178],"class_list":["post-846","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-uncategorized","category-looking-into-something","category-preteach","tag-analysis","tag-fitting","tag-instructions","tag-scattering-pattern"],"jetpack_publicize_connections":[],"jetpack_featured_media_url":"https:\/\/lookingatnothing.com\/wp-content\/uploads\/2013\/09\/Cyl_intensity.png","jetpack_shortlink":"https:\/\/wp.me\/p1gZ2v-dE","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/posts\/846","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/comments?post=846"}],"version-history":[{"count":8,"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/posts\/846\/revisions"}],"predecessor-version":[{"id":858,"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/posts\/846\/revisions\/858"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/media\/848"}],"wp:attachment":[{"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/media?parent=846"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/categories?post=846"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/lookingatnothing.com\/index.php\/wp-json\/wp\/v2\/tags?post=846"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}