Dear Ian, On Thu, 10 Apr 2003, Ian Drake wrote:
... I would like to calculate and plot the fourier transform of the phase corrected chi function. Can I do this using the UWXAFS software suit? If not, can you recommend other software. I have Mathematica....do you recommend its use in the analysis. My curiosity on this issue stems from the work of Lee and Beni in PHYS. REV. B, 15, 1977, pg.2862.
Thanks for the question. Yes, phase-corrected Fourier transforms can be done. I hope it's OK if I CC this reply to the Ifeffit mailing list, as this is a fairly common question about XAFS analysis... Also, since you asked about UWXAFS software, I'll point out that Ifeffit is the successor to the UWXAFS programs, and includes those programs and quite a bit more capability, including a few GUIs. See: http://cars.uchicago.edu/ifeffit/ for more details. The longer answer about Phase-corrected FTs: Phase-corrected FTs are not always as good as you might hope. They don't give real radial distribution functions. They do give somewhat cleaner XAFS FTs and put peaks in |chi(R)| closer to the corresponding distance. Of course, to do these, you need to know the phase. Also, since you can do only one phase-shift correction, you can't correct all shells or multiple-scattering contributions. Anyway, lots of people like phase-corrected FTs, especially for making plots to show non-XAFS people, and it usually work fine for the first-shell. With ifeffit, there are a couple ways to do phase-corrections. The simplest way is to name the central atom/edge (eg, 'Cu K') to the fftf() command. In ifeffit commands, this looks like: read_data(file= 'my_chi.dat', type='chi') fftf(my_chi.chi, kmin=2, kmax=15, dk=2, kweight=2, pc_caps = 'Cu K') ^^^^^^^^^^^^^^^^ (pc_caps is for Phase-Correction, Central Atom Phase Shift), This will do only the central-atom portion of the phase-shift, so will be incomplete, but is often most of the correction, will work for the whole spectra, and doesn't require knowing much about the system. Athena (the excellent GUI to ifeffit from Bruce Ravel) can do these central-atom phase-corrections very easily (ie, no typing commands). I think Sam Webb's SixPack doesn't do phase-corrected FTs, but I could be wrong. If you know have a Feff calculation for a scattering path, in the form of a feffNNNN.dat file, you can use this for a more-complete phase-corrected FT, like this: path(1, feffNNNN.dat) fftf(my_chi.chi, kmin=2, kmax=15, dk=2, kweight=2, pc_feff_path = 1) ^^^^^^^^^^^^^^^^ This version includes both central-atom and scattering-atom phase-shifts, which will give a much better phase-correction, but assumes you know quite a bit about the system. I think athena doesn't do this form of phase-correction, and I don't remember if artemis does this (yet??). Feffit will aslo write out phase-corrected FTS (the full kind, using the first feffNNNN.dat file listed in feffit.inp) if the keyword 'pcout=1' in in feffit.inp. Finally, some people use phase-corrected FTs for analysis, especially to pick R_NN at a place where Im[chi(R)] = 0. I'm not sure what the accuracies of this are, but pelase note that Ifeffit (and feffit) swaps real and imaginary parts for its FT, so that it is actually Re[chi(R)] that is ~=0 at R=R_NN for these FTs. Hope that helps, but don't hesitate to ask again if it's not clear enough. --Matt Newville
On Sunday 13 April 2003 03:06 pm, Matt Newville wrote:
This will do only the central-atom portion of the phase-shift, so will be incomplete, but is often most of the correction, will work for the whole spectra, and doesn't require knowing much about the system. Athena (the excellent GUI to ifeffit from Bruce Ravel) can do these central-atom phase-corrections very easily (ie, no typing commands).
Hi, I just wanted to quickly chime in on this topic. For what it's worth -- I am not a huge fan of phase corrections. While they do, in many cases, make the chi(R) look more like a radial distribution function, they do require some more explanation when showing your data to a non-xafs person. And, while a pc-chi(R) looks more like the corresponding RDF, it most certainly is not an RDF. Athena does only a central atom correction for the reason that Matt gave. Athena assumes only that your data is XAFS data and does not assume any knowledge of the system. Thus, the only thing that Athena is sure to know is what the absorbing atom is. That is determined from the edge energy. Artemis, on the other hand, has more information to work with. When you choose to do phase corrections in Artemis, you choose one of the feffNNNN.dat files included in the fit as the phase correction standard. Thus, in Artemis, phase correction includes the scattering portion as well as the central atom portion. In this screenshot: http://leonardo.phys.washington.edu/~ravel/software/exafs/images/artemis/art... the "Path to use for phase corrections" menu consists of "None" (i.e. no pc), and each of feff0001.dat through feff0005.dat. You might then choose to use feff0001.dat (i.e. the first shell in the example shown) for phase corrections. Note that the corrections are only for plotting. Fitting in artemis is done with the non-phase-corrected data and theory. Upcoming news flash: I have spent the last week getting multiple data set fits working in Artemis. Ifeffit has supported them for a long time, but now Artemis will provide an interface to that extremely useful tool. I hope to have the new release out soon. B -- Bruce Ravel ----------------------------------- ravel@phys.washington.edu Code 6134, Building 3, Room 222 Naval Research Laboratory phone: (1) 202 767 5947 Washington DC 20375, USA fax: (1) 202 767 1697 NRL Synchrotron Radiation Consortium (NRL-SRC) Beamlines X11a, X11b, X23b, X24c, U4b National Synchrotron Light Source Brookhaven National Laboratory, Upton, NY 11973 My homepage: http://feff.phys.washington.edu/~ravel EXAFS software: http://feff.phys.washington.edu/~ravel/software/exafs/
participants (2)
-
Bruce Ravel
-
Matt Newville