[CCP14 Home: (Frames | No Frames)] CCP14 Mirrors: [UK] | [CA] | [AU]

## Lp correction in GSAS

The CCP14 Homepage is at http://www.ccp14.ac.uk

[Back to Problems and Solutions] [Back to GSAS Hints/Resources]

[The reference to use for GSAS in any resulting publications is: A.C. Larson and R.B. Von Dreele, "General Structure Analysis System (GSAS)", Los Alamos National Laboratory Report LAUR 86-748 (1994).]

 ```Date: Sun, 9 Jul 2000 14:51:58 +1000 (EST) From: Marek Schmidt To: Rietveld_l@ill.fr Dear All, I have two questions concerning the polarization factor in GSAS. I am analyzing synchrotron x-ray diffraction patterns and to work out the diffractometer constants I refined a silicon pattern. My questions are: 1.What is the formula for the RBVD Lp correction? Is it the third formula on page 125 of the manual? 2. The pattern was refined using the RBVD correction. The refined polarization fraction POLA reaches a value of 1.14 (the initial value was set to 0.95) is it possible? The experiment was carried out in Debye-Scherrer geometry with an incident beam monochromator and the patterns were collected using imaging plates. The camera radius is 573mm, the sample was placed in 0.3mm glass capillary. The refinement includes absorption correction and corrections to the silicon scattering factor (f' and f"). Wavelength= 1.796A, the angular range is 5-150deg. Best regards M. Schmidt ====================================================== Marek W. Schmidt Department of Applied Mathematics Research School of Physical Sciences and Engineering Australian National University Canberra ACT 0200, AUSTRALIA ====================================================== ```

 ```Date: Mon, 10 Jul 2000 16:28:35 -0400 From: "Brian H. Toby" [Brian.Toby@nist.gov] Organization: NIST Center for Neutron Research X-Accept-Language: en To: Marek Schmidt , Rietveld_l@ill.fr Subject: Re: Lp correction in GSAS Marek, There are three polarization functions listed in the manual, but only 2 when I try to select the function in EXPEDT. I seem to remember that the 1st two functions are equivalent, other than scaling, but I am not going to work this mystery today. As I recall, with IPOLA=0 (RBVD Lp correction) POLA sets the relative amount of in-plane vs out-of-plane polarization coming off the incident monochromator. Thus for you, POLA should be slightly less than 1.0 (synchrotrons tend to be highly polarized, but the actual amount depends on many factors). I have used 0.95 to 0.99 for my work at X7A. The bad (or good) issue here is that the P part of the LP correction is virtually indistinguishable from the Debye-Waller 2theta curve, which means if you refine it, the displacement parameters (aka "thermal" factors) couple to POLA and you can refine to an unrealistic value, like >1. If you fix POLA to a wrong value, your displacement parameters shift to compensate. For that matter, absorption corrections have ~ the same form, too, particularly since you have data over such a small Q range, so that will also interact with POLA and displacement parameters. Pick a value of POLA and fix it. Someone at the beamline you used must have a rough idea of the right value. If you are really adventurous, pick a few values and compare the results of the different refinements. That will give you a feeling for the impact of POLA on your results. Brian ******************************************************************** Brian H. Toby, Ph.D. Leader, Crystallography Team Brian.Toby@NIST.gov NIST Center for Neutron Research, Stop 8562 voice: 301-975-4297 National Institute of Standards & Technology FAX: 301-921-9847 Gaithersburg, MD 20899-8562 http://www.ncnr.nist.gov/xtal ********************************************************************```

 ```Date: Tue, 11 Jul 2000 11:44:20 +0100 (BST) From: Jon Wright [jpw22@cus.cam.ac.uk] To: Rietveld_l@ill.fr Subject: Re: Lp correction in GSAS X-Loop: rietveld_l@ill.fr X-Sequence: 12 On Mon, 10 Jul 2000, Brian H. Toby wrote: > There are three polarization functions listed in the manual, but only 2 > when I try to select the function in EXPEDT. I seem to remember that the > 1st two functions are equivalent, other than scaling, but I am not going > to work this mystery today. My apologies for such a public display of ignorance, I have a related mystery to ask about. Having had a quick look at the manual, all three corrections are for parallel geometry, I'm guessing "no monochromator", then "incident beam mono." and then "diffracted beam mono.". OK. The question is what about flat plate transmission (STOE) type data versus Bragg-Brentano versus Debye-Sherrer? Doesn't a factor of cos(theta) at least turns up in the transmission data. This is where I start to feel rather stupid, but how can you specify these cases? It doesn't appear to come into the instrument parameter file or the lorentz corrections dicussed here? Any hints or advice would be appreciated. Probably Bob has the best idea of what GSAS is doing, and how it's doing it. I've heard rumours of manual "corrections" being applied to data before presenting it for refinement, which sounds ugly. Fullprof has a flag for transmission geometry, hence my question. Thanks in advance, Jon Wright PS : If someone could refer me to simple derivations of the LP corrections for all the various geometries, I'd be very grateful. ============================================================================ Dept. of Chemistry, Lensfield Road, Cambridge, CB2 1EW, UK```

 ```Date: Tue, 11 Jul 2000 08:50:25 -0600 To: Jon Wright [jpw22@cus.cam.ac.uk] From: vondreele@lanl.gov (Bob Von Dreele) Subject: Re: Lp correction in GSAS Cc: Rietveld_l@ill.fr X-Loop: rietveld_l@ill.fr X-Sequence: 13 Dear Jon (and others), The "standard" reference for X-ray Lp corrections is the paper by L. Azaroff, Acta Cryst., 8,701-704 (1955). It applies equally well for both incident beam and diffracted beam monochromator cases and is independent of choice of Debye-Scherrer, Bragg-Brentano or for that matter single crystal diffraction (hence also applies to transmission geometry - but there may be a separate "foot print" effect that depends on the optics. See below.). The L (Lorentz) part is applied to CW neutron data (usually Debye-Scherrer) since there is no neutron polarization effect. There are three "options" in GSAS for x-ray Lp corrections. The first two are equivalent except for scaling and the third is no correction at all; some folks wanted this as they apply the correction beforehand. However, they must do proper error propagation for this correction and use the "ESD" or "ALT" form for the raw data input to GSAS. The third equation in the GSAS Manual was incorrect, it was deleted sometime ago. The crucial thing to establish for transmission geometry is whether the illuminated (and seen!) sample volume is independent of scattering angle. If it is independent then either GSAS Lp correction will work just fine. If not, then an additional correction must be applied for the change in sampling volume - this correction is not in GSAS as I don't know what it is. It should also be noted that the sampling volume error can occur in Bragg-Brentano geometry when the sample isn't "infinitely thick" or doesn't cover completely the beam foot print at all measured scattering angles. This is a particular problem for diffraction from organic solids with low absorption or for thin "smears" of sample on a glass slide or samples that are too small for complete beam coverage at low angles. GSAS doesn't have anything to deal with these elementary experimental errors. Bob Von Dreele```

[Back to Problems and Solutions] [Back to GSAS Hints/Resources]

 [CCP14 Home: (Frames | No Frames)] CCP14 Mirrors: [UK] | [CA] | [AU]

If you have any queries or comments, please feel free to contact the CCP14