Thank you Bruce and Shelley! I just wanted to thank Matt, Bruce and everyone involved with IFEFFIT/Artemis/Athena... it has been a Godsend.... I couldn't have stumbled on it at a better time (I'm frantically trying to analyze about 4 beamtimes worth of SEXAFS data for papers and my thesis and WINXAS wasn't really working out for me)!!!!! Does anyone have any suggestions for additional literature in regards to better understanding the quality of fits and the parameters involved?.... or perhaps know when we might be able to expect a updated version of the Artemis literature? (What is currently there has helped me sooo much that I've been checking for an update every couple of weeks hoping to get see sections 2.2-3.2..... ) Also.... I guess I've realized that only knowing IFEFFIT with Artemis I kinda take things for granted (like I didn't know that IFEFFIT was having problems, in certain cases since I had only been checking the "Results" and the error messages on the bottom of the window....) which meant I missed errors like this in the IFEFFIT (since I didn't notice anything unusual in the "Results"): ............ # Do the fit! ## fitting difference of Cov_0_06ML_bl1435_111703.chi and one or more paths ... ## === data set #1 of 1 feffit(1-8, group=data0_fit, chi=data0.chi, k=data0.k, kweight=2.00, rmin=1, rmax=4, kmin=2.5, kmax=15, dk=2, kwindow=kaiser-bessel, fit_space=R, do_bkg=no, data_set=1, data_total=1) feffit fitting 1 data sets fitting ... fitting ... fitting ... fitting ... fitting ... fitting ... fitting ... fitting ... no uncertainties estimated! fit gave a warning message: too many fit iterations. try again with better guesses or a simpler problem. ............................................... I guess my question is.... is there anything else I should be looking for in the IFEFFIT printout that I might not otherwise notice by just using Artemis and looking at the error messages and "Results"? Thanks for the input! Sincerely, Kristine Witkowski Northwestern University
On Wednesday 07 January 2004 02:51 pm, k-kupiecki@northwestern.edu wrote:
I just wanted to thank Matt, Bruce and everyone involved with IFEFFIT/Artemis/Athena... it has been a Godsend.... I couldn't have stumbled on it at a better time (I'm frantically trying to analyze about 4 beamtimes worth of SEXAFS data for papers and my thesis and WINXAS wasn't really working out for me)!!!!!
Well, thank you very much, Kristine. Email like this one of the great rewards of working on the codes and one of things that keeps me going. Best of luck with the thesis!
Does anyone have any suggestions for additional literature in regards to better understanding the quality of fits and the parameters involved?....
Have you looked through all the stuff linked to on Matt's XAFS page? http://cars9.uchicago.edu/xafs/
or perhaps know when we might be able to expect a updated version of the Artemis literature? (What is currently there has helped me sooo much that I've been checking for an update every couple of weeks hoping to get see sections 2.2-3.2..... )
Uhhhh.... hmmmm.... okay, everyone, move along.... nothing to see here ...
Also.... I guess I've realized that only knowing IFEFFIT with Artemis I kinda take things for granted (like I didn't know that IFEFFIT was having problems, in certain cases since I had only been checking the "Results" and the error messages on the bottom of the window....) which meant I missed errors like this in the IFEFFIT (since I didn't notice anything unusual in the "Results"):
............. # Do the fit!
## fitting difference of Cov_0_06ML_bl1435_111703.chi and one or more paths ... ## === data set #1 of 1 feffit(1-8, group=data0_fit, chi=data0.chi, k=data0.k, kweight=2.00, rmin=1, rmax=4, kmin=2.5, kmax=15, dk=2, kwindow=kaiser-bessel, fit_space=R, do_bkg=no, data_set=1, data_total=1) feffit fitting 1 data sets fitting ... fitting ... fitting ... fitting ... fitting ... fitting ... fitting ... fitting ... no uncertainties estimated! fit gave a warning message: too many fit iterations. try again with better guesses or a simpler problem.
................................................
I guess my question is.... is there anything else I should be looking for in the IFEFFIT printout that I might not otherwise notice by just using Artemis and looking at the error messages and "Results"?
That's the most important one to keep an eye on. Artemis is pretty good at making sure that all the other ifeffit commands are run in a way that will not make ifeffit unhappy. However, the only way to know if ifeffit is happy fitting is to run the fit. Certainly, if Artemis every starts behaving in some way out of the ordinary, it would be prudent to scroll through the ifeffit buffer messages. Ifeffit may have told you something that will help you solve your own problem. Or it may help you make a bug report that will better help solve the problem. As you no doubt recall, you and I had an exchange on this topic last month here on the mailing list. Here's my response: http://millenia.cars.aps.anl.gov/pipermail/ifeffit/2003-December/000602.html The suggestions I made then still seem appropriate. I started noticing the sort of situation that leads to the error message that you cite before I went on holiday. Artemis really should be looking for that error message and post a warning to the user when it happens. I'll put that on the list of things to do. 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 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/
A simple question. I am have been doing multiple edge fits of different materials for a while in artemis (diana). Now, I would like to analyze a material with inequivalent sites for one of the edges. I understand how to read in the data and have read in the different feff calculations (for each inequivalent site), but what I just want to confirm that artemis/diana/ifeffit will do the "right" thing, namely coherently add the contributions *within* each feff calculation and add the intensities of each feff calculation -- e.g. I don't think there is any interference between the XAFS signal on different crystallographic sites -- or even if there is it is ignorable at least for first neighbor interactions. If there is anyone out there with experience in doing this (or is the right thing done automatically?), drop me (or the list) a line. Thanks a lot! Dr. Paul Fons Senior Researcher National Institute for Advanced Industrial Science & Technology METI Center for Applied Near-Field Optics Research (CANFOR) AIST Central 4, Higashi 1-1-1 Tsukuba, Ibaraki JAPAN 305-8568 tel. +81-298-61-5636 fax. +81-298-61-2939 email: paul-fons@aist.go.jp The lines below are in a Japanese font 〒305−8568 茨城県つくば市東1−1−1 つくば中央第4 近接場光応用工学センター ポール・フォンス主任研究官
Hi Paul, Actually, the inequivalent sites do interfere coherently. EXAFS chi's aren't like quantum mechanical wave functions, for which only relative phase matters. The phase gives information about E0, and that's a physically meaningful quantity. In other words, if one of the sites is producing a peak at 7200 eV and another is producing a trough, they will tend to cancel. So yes, Artemis and Ifeffit do the right thing, but the right thing is to add everything up coherently regardless of which feff calculation it comes from. A few other notes on this kind of fit: you must, of course, weight the amplitudes of each feff calculation by the relative number of sites of that type. Also, it is quite possible that you'll need different e0's for the different sites. --Scott Calvin Sarah Lawrence College
A simple question. I am have been doing multiple edge fits of different materials for a while in artemis (diana). Now, I would like to analyze a material with inequivalent sites for one of the edges. I understand how to read in the data and have read in the different feff calculations (for each inequivalent site), but what I just want to confirm that artemis/diana/ifeffit will do the "right" thing, namely coherently add the contributions *within* each feff calculation and add the intensities of each feff calculation -- e.g. I don't think there is any interference between the XAFS signal on different crystallographic sites -- or even if there is it is ignorable at least for first neighbor interactions. If there is anyone out there with experience in doing this (or is the right thing done automatically?), drop me (or the list) a line. Thanks a lot!
Thanks for the message. I guess what I wanted to say was that from what I understand due to lifetime effects the "effective" radius about the absorbing site for which the outgoing scattered wave "exists" is on the order of 2 nn (I have big lattice constant material). It would seem to me then that the inequivalent sites being physically separated by more than this distance (we are talking about crude approximations here) would not interfere. In reality, with bond lengths on the order of 3 Angstroms and 2nn along the lines of 4.5 Angstroms, my case is somewhere in between the coherent interaction and the dilute dopant -- e.g. every site is by itself -- extreme. Certainly, in the latter case it would be wrong to add the signals together coherently! What is the consensus for this sort of thing? The different E0's is a good point too. Paul On 2005/01/14, at 16:16, Scott Calvin wrote:
Hi Paul,
Actually, the inequivalent sites do interfere coherently. EXAFS chi's aren't like quantum mechanical wave functions, for which only relative phase matters. The phase gives information about E0, and that's a physically meaningful quantity. In other words, if one of the sites is producing a peak at 7200 eV and another is producing a trough, they will tend to cancel.
So yes, Artemis and Ifeffit do the right thing, but the right thing is to add everything up coherently regardless of which feff calculation it comes from.
A few other notes on this kind of fit: you must, of course, weight the amplitudes of each feff calculation by the relative number of sites of that type. Also, it is quite possible that you'll need different e0's for the different sites.
--Scott Calvin Sarah Lawrence College
A simple question. I am have been doing multiple edge fits of different materials for a while in artemis (diana). Now, I would like to analyze a material with inequivalent sites for one of the edges. I understand how to read in the data and have read in the different feff calculations (for each inequivalent site), but what I just want to confirm that artemis/diana/ifeffit will do the "right" thing, namely coherently add the contributions *within* each feff calculation and add the intensities of each feff calculation -- e.g. I don't think there is any interference between the XAFS signal on different crystallographic sites -- or even if there is it is ignorable at least for first neighbor interactions. If there is anyone out there with experience in doing this (or is the right thing done automatically?), drop me (or the list) a line. Thanks a lot!
_______________________________________________ Ifeffit mailing list Ifeffit@millenia.cars.aps.anl.gov http://millenia.cars.aps.anl.gov/mailman/listinfo/ifeffit
Dr. Paul Fons Senior Researcher National Institute for Advanced Industrial Science & Technology METI Center for Applied Near-Field Optics Research (CANFOR) AIST Central 4, Higashi 1-1-1 Tsukuba, Ibaraki JAPAN 305-8568 tel. +81-298-61-5636 fax. +81-298-61-2939 email: paul-fons@aist.go.jp The lines below are in a Japanese font 〒305−8568 茨城県つくば市東1−1−1 つくば中央第4 近接場光応用工学センター ポール・フォンス主任研究官
Hi Paul, I think we're getting confused over what we each mean by "interfere." Whether sites are equivalent or inequivalent, neighboring or far apart, the wave functions do not interfere at all, in part because the odds of two x-ray photons being absorbed by sites close to each other at nearly the same time is very very small--if it wasn't, your material would disintegrate in an instant! The idea behind the path expansion used in EXAFS analysis by ifeffit, however, is that the chi(k) for each possible scattering event by a single electron can be computed separately, and then the results added (I'll defer to John or Matt or whoever for a clear explanation of why this is justified). In other words, the EXAFS spectra can be thought of as if there were a contribution from an electron scattering off a near-neighbor, another scattering off a 2nd-nearest-neighbor, another doing a multiple-scattering thing, etc. (all weighted, of course, by the relative contributions of these scattering events). These are literally just added. Given that system, if there really are multiple sites that are each having their own scattering events, then it is perfectly appropriate to treat them in exactly the same way--just add. So, depending on your terminology, either we are combining all scattering paths coherently or we are adding none of them coherently, regardless of whether they correspond to the same absorbing site or not. It would not be correct to add, say, the magnitudes of the FT's for the separate contributions...the phase of the FT indicates where the peaks and troughs are in k-space (and thus energy-space), and it certainly matters whether a given scattering path has enhanced or suppressed absorption at a given energy! Hope that helps... --Scott Calvin Sarah Lawrence College P.S. I've published several analyses of systems with inequivalent sites. If you'd like a pdf of one of those papers, let me know.
Thanks for the message. I guess what I wanted to say was that from what I understand due to lifetime effects the "effective" radius about the absorbing site for which the outgoing scattered wave "exists" is on the order of 2 nn (I have big lattice constant material). It would seem to me then that the inequivalent sites being physically separated by more than this distance (we are talking about crude approximations here) would not interfere. In reality, with bond lengths on the order of 3 Angstroms and 2nn along the lines of 4.5 Angstroms, my case is somewhere in between the coherent interaction and the dilute dopant -- e.g. every site is by itself -- extreme. Certainly, in the latter case it would be wrong to add the signals together coherently! What is the consensus for this sort of thing? The different E0's is a good point too.
Paul
Re: [Ifeffit] Inequivalent sites and multiple shell fitsPaul, may be it will be more clear what you mean by interference of inequivalent sites if you explain it on the example of equivalent sites. In a bulk fcc metal, contributions to EXAFS from neighboring (equivalent) sites do not interfere - and it is why we can use scattering amplitudes and phases extracted from model compounds (or FEFF theory) with one structure to model unknown compound with a (not too) different structure. Had it not been true, the chemical transferability of scattering phase and amplitude wouldn't be working but it does in most cases. So, if it is working for equivalent sites, why should it not for inequivalent sites? Moreover, this possibility (of interference) is almost always ruled out in FEFF at the state when we assign unique potentials to different atoms, the absorber (0) and its neighbors (1). Sometimes we assign a unique potential 0 to the central atom, and, another unique potential 1 to its first nearest neighbors, to correct for the muffin tin approximation, while the rest of the atoms have the same potential 2 (if there is just one element in the sample). For fcc copper, both models give very similar results for FEFF paths and for the fits. The reason the central atom has different potential than the 1NN is because it is ionized. Had its 1NN been ionized also, its potential would be different too, but we know that we can get away just fine by assigning a neutral atom potential to all the neighbors. Therefore, if FEFF does not worry about interference, how can it be taken into account in principle, and why? Regards, Anatoly Anatoly Frenkel Yeshiva University [Anatoly Frenkel] -----Original Message----- From: ifeffit-bounces@millenia.cars.aps.anl.gov [mailto:ifeffit-bounces@millenia.cars.aps.anl.gov]On Behalf Of Scott Calvin Sent: Friday, January 14, 2005 12:49 PM To: XAFS Analysis using Ifeffit Subject: Re: [Ifeffit] Inequivalent sites and multiple shell fits Hi Paul, I think we're getting confused over what we each mean by "interfere." Whether sites are equivalent or inequivalent, neighboring or far apart, the wave functions do not interfere at all, in part because the odds of two x-ray photons being absorbed by sites close to each other at nearly the same time is very very small--if it wasn't, your material would disintegrate in an instant! The idea behind the path expansion used in EXAFS analysis by ifeffit, however, is that the chi(k) for each possible scattering event by a single electron can be computed separately, and then the results added (I'll defer to John or Matt or whoever for a clear explanation of why this is justified). In other words, the EXAFS spectra can be thought of as if there were a contribution from an electron scattering off a near-neighbor, another scattering off a 2nd-nearest-neighbor, another doing a multiple-scattering thing, etc. (all weighted, of course, by the relative contributions of these scattering events). These are literally just added. Given that system, if there really are multiple sites that are each having their own scattering events, then it is perfectly appropriate to treat them in exactly the same way--just add. So, depending on your terminology, either we are combining all scattering paths coherently or we are adding none of them coherently, regardless of whether they correspond to the same absorbing site or not. It would not be correct to add, say, the magnitudes of the FT's for the separate contributions...the phase of the FT indicates where the peaks and troughs are in k-space (and thus energy-space), and it certainly matters whether a given scattering path has enhanced or suppressed absorption at a given energy! Hope that helps... --Scott Calvin Sarah Lawrence College P.S. I've published several analyses of systems with inequivalent sites. If you'd like a pdf of one of those papers, let me know. Thanks for the message. I guess what I wanted to say was that from what I understand due to lifetime effects the "effective" radius about the absorbing site for which the outgoing scattered wave "exists" is on the order of 2 nn (I have big lattice constant material). It would seem to me then that the inequivalent sites being physically separated by more than this distance (we are talking about crude approximations here) would not interfere. In reality, with bond lengths on the order of 3 Angstroms and 2nn along the lines of 4.5 Angstroms, my case is somewhere in between the coherent interaction and the dilute dopant -- e.g. every site is by itself -- extreme. Certainly, in the latter case it would be wrong to add the signals together coherently! What is the consensus for this sort of thing? The different E0's is a good point too. Paul
I think this echos Scott's and Anatoly's answer, and some of this may be obvious, but I'll take a stab at this. Each absorption of an x-ray is a completely independent event from other absorption events. The excited state lasts for something like 1 femtosecond. We never get more than 10^15 photons/sec, so if the x-ray beam was continuous, we'd never have more than one excited state in the sample at any one time. Even with a very bright and pulsed source (and synchrotrons are pulsed, of course), we can safely assume that having two excited atoms of interest (that is, other atoms can absorb too) at any one time is very rare. When this does happen, the excited atoms will be randomly distributed in the sample, and so very far from each other. (There is very interesting work on doubly-excited atoms, but as far as I know they're still using x-ray + laser, not two x-rays, and I'm getting off topic). The point is that we can very safely assume that EXAFS is a sum over independent absorption events, each happening without any other excited atoms nearby. Since all atoms with atomic number Z have an equal absorption probability and since the full spectra is made up from billions of absorption events (yeah, it's sometimes less), all atoms of number Z are well sampled. The total EXAFS really is then just a sum over the absorbing atoms: chi = Sum_absorbing_atoms[Sum_paths(chi_absorber_path)] If all the absorbing atoms look the same (eg, fcc metal, or the A or B site in a perovskite structure), the first sum is trivial. If there are multiple sites that the central atom can occupy, you just do the weighted sum of sites, where the weights are the occupation probabilities for that site. In a glass with a range of absorbing atom "sites", you could turn the sum into an integral or use a histogram, but it's usually a pretty simple sum. Of course, each site would have its own set of scattering paths (ie, it's own Feff calculation), and each path with that absorbing atom needs the weighting appropriate for that absorber. Now, Feff/Feffit/Ifeffit/Artemis/Diana/SixPack (my goodness, that's getting long!) does only the sum over paths: chi = Sum_paths(chi_absorber_path) To handle a system with multiple absorbing sites, you then need to multiply each path amplitude by the fractional probability of the corresponding central atom: chi = Sum_paths ( weight_absorber * chi_absorber_path) where the weights are fractional, so that Sum (weight_absorber) = 1 For example, if you have 2 different sites, one that's twice as common as the other (say, Fe3O4), you'd have a bunch of path statements like # set of paths for Absorber A: probability = 2x Path 101 Asite/feff0001.dat Amp 101 (2/3) * S02 Path 102 Asite/feff0002.dat Amp 102 (2/3) * S02 Path 103 Asite/feff0003.dat Amp 103 (2/3) * S02 # set of paths for Absorber B: probability = 1x Path 201 Bsite/feff0001.dat Amp 201 (1/3) * S02 Path 202 Bsite/feff0002.dat Amp 202 (1/3) * S02 I think that's exactly what Scott and Anatoly are getting at. Hope that helps, --Matt
Wow. An interesting discussion has started. I am afraid my original question was not very clear (my fault). I just wanted to be sure that athena (diana) was doing the "right thing", e.g adding the contributions from separate feff calculations, e.g. treating them as incoherent sums where the phase information is not preserved. In retrospect, it seems obvious that it must be doing so, but in a moment of doubt I dashed off a quick email. Thanks for the help.
participants (6)
-
Anatoly Frenkel
-
Bruce Ravel
-
k-kupiecki@northwestern.edu
-
Matt Newville
-
Paul Fons
-
Scott Calvin