Dear Prof. Ravel, Thanks so much for your answer, which remove the cloud partially above me. I actually tried many models, from crystal to the opposite 6O only (different distance, second shell C,P ...). The example shown in the email are two extreme cases, which single P has contribution like 6 P. I also tried 6O+6(PO3), which is close to a crystal structure(6 O are not same distance from center, same for 6(PO3)) ---- when I decease # of PO3 from 6 to 0, FEFF calculations of O change gradually --- not like ' single P has contribution like 6 P'. It make me wonder, 1, FEFF takes Z atoms at equal distance as same. It does not care some have neighbor A, some have neighbor B. Should I worry about it? and when? 2, Does FEFF consider about valence state (like Ur 3 or Ur5 in center)? Model closer to sample, better result. Of course crystal structure is closer to real sample than models like 6O do. In practice, If I only consider single scattering and fit to n shell, can I say that it is good enough to have a model of n+1 shell? For instance I want to fit sample with O and P shells, it is good enough to have a model with Cd -O-P-O3, since O3 shell will influence P(like P atoms influence 1st shell O). Best, Jack Song On 2011-04-11 15:05, Bruce Ravel wrote:
Jack,
You have rediscovered the conventional wisdom about properly bounding a cluster for calculation with Feff. When there are not enough atoms on the periphery of the cluster to constrain the size of the muffin tins, the muffin tins tend to be much too large. As you can see this has consequence for the scattering function calculated by Feff.
In your case, the muffin tins of the oxygen atoms are about 0.73 Angstroms when including even a single P atom on the periphery of the cluster. For the first calculation Rmt is 1.06 -- substantially larger.
What's interesting is that a single P atom is sufficient to make the calculation of the O contribution work out reasonably well.
Have you read Shelly's book chapter review of EXAFS? In her example on determining a second shell atom type, she runs Feff calculations on known crystal structures with various second shell atoms. She then takes just the kinds of paths she needs to model different contributions to the EXAFS. The reason she adopted this strategy rather than trying a small, notional cluster of the sort the you have made in your example is because constraining the muffin tins is quite important. Indeed, the potentially counter-intuitive strategy of snarfing paths from some related crystal structure works well precisely because it is a way of guaranteeing sensible muffin tin radii.
HTH, B
P.S. You can see the muffin tin radii by opening the .apj file as a zip file (which is what it is!) and looking at the misc.dat file in the dataX.feffY folders.
On Monday, April 11, 2011 02:33:54 pm goodhei8@gmail.com wrote:
Dear IFEFFIT Members,
I have a strange problem with FEFF calculation. I am not sure if I made some mistakes in FEFF of input file, or FEFF did not handle it well.
For Oxygen paths calculated with FEFF in Artemis by 3 input files, the results for O are supposed to be close, but I got a big difference, when I transform them into r space(k 2.5~10,kweight=2 in attachment).
There are three input files for FEFF: #1, 6O only on xyz axis #2, + 1P besides 6O #3, + 6P besides 6O
For O, #2 and #3 are close(almost same), but #1 is much bigger. (It would make sense to me if #1 #2 are close.)
Thank a lot.
Best,
Jack Song