Hi Shelly, I'll look up the reference you provide when I get back to a computer that has access to Physica, but for now a question that might be of interest to the list: The MS scheme you outline below seems to treat the E0 due to the scattering species (I'll call that the "bonus" E0 for lack of a better term) as being the average of the bonus E0's for all the scattering atoms in the path. Thus you give the Cu-O1-Cu1-Cu path as using the average of the O1 and the Cu1 bonus E0's, while the Cu-O1-Cu-O1 path uses just the O1 bonus E0. I would naively expect we should be using the sum, rather than the average. After all, if this bonus E0 is a correction to the shift caused by the muffin-tin potential FEFF has calculated, shouldn't a MS path that bounces off the same atom twice have to use twice the correction? I feel like there's something basic I'm missing here... Also, what is your experience as to the actual differences in bonus E0's that you get in ionic compounds? Right now one of my students is working with Sr(NO3)2... should we really expect the difference in bonus E0's between the Sr2+ and the oxygen atoms to be more than 12 eV, or does some kind of charge overlap from the nitrate ions actually make the strontium ions behave as if they are less than +2? Thanks in advance for any comments you have! --Scott Calvin Sarah Lawrence College <fontfamily><param>Arial</param><color><param>0000,0000,ffff</param><smaller> </smaller></color></fontfamily><<<<<<<< Shelly wrote:
<fontfamily><param>Arial</param><color><param>0000,0000,ffff</param><smaller>I wanted to add another comment about using multiple E0's. One of the best formulas that I have found involves assigning different atom types different E0's. This method seems to work well for M.S. paths. Here is an example of how it goes. I think that Daniel Haskel was the first to publish this type of E0 stuff. </smaller></color></fontfamily><color><param>0000,0000,ffff</param><fontfamily><param>Times New Roman</param>"Single and multiple scattering XAFS in BaZrO3: a comparison between theory and experiment" <italic>D. Haskel, B. Ravel, M. Newville and E. A. Stern </italic>Physica B 208&209 (1995) 151 </fontfamily></color> <fontfamily><param>Arial</param><color><param>0000,0000,ffff</param><smaller>Lets pretend that we have a copper oxide that we want to model. First shell Oxygen called O1, Second shell Cu called Cu1, Third shell Oxygen called O2. This method would use three E0's. One for the first shell oxygen atoms, Eo1. One for all the Cu atoms Ecu. and one for all other oxygen atoms Eo2. The Cu-O1 path has E0=Eo1. The Cu-Cu1 path has E0=Ecu. The Cu-O2 path has E0=Eo2. The Cu-O1-Cu-O1 path has E0=Eo1, (only scattering from O1 atom). The Cu-O1-Cu1-Cu path has E0=0.5*Eo1+0.5*Ecu. (scattering once from O1 and once from Cu1) ... And so forth and so on. </smaller></color></fontfamily> <fontfamily><param>Arial</param><color><param>0000,0000,ffff</param><smaller>Now if you set something like this up. You need to look at the results carefully and see if it works. Generally, scattering from an atom that has lost an electron due to an ionic bond will have a positive deltaE, and scattering from an atom that has gained an electron will have a negative deltaE. Of course this negative/positive deltaE is added to the alignment deltaE. If the values for Eo2 and Ecu are the same within the uncertainties of the measurement then you need to make them the same variable, so then you will only have two energy shifts, but you can still use it. This is most often the case as Matt mentioned but I always like to test it out. In addition, one electron transfer is about a shift by 6 eV, so if these values differ by more than say 10 eV something is not right. Since the shift is generally smallish, say 2 eV, your fit has to be really accurate to "see" these shifts. In particular they seem to be most important when you have all the deltaR values constrained in some way so that the error in deltaR is really really small. For example all deltaR values are related to unit cell parameters and atomic positions or and overall expansion term. </smaller></color></fontfamily>