October 2005 - Ifeffit - millenia.cars.aps.anl.gov

horae Debian package
by Carlo Segre 19 Oct '05

19 Oct '05

Hello All: Because of the new dependency on libmath-round-perl. The Debian package of horae 056 will have to be installed manually. It can be found at the deb http://fermi.phys.iit.edu/debian/woody binary/ deb http://fermi.phys.iit.edu/debian/sarge binary/ deb http://fermi.phys.iit.edu/debian/sid binary/ depending on the distribution you are using. Cheers, Carlo -- Carlo U. Segre -- Professor of Physics Associate Dean for Special Projects, Graduate College Illinois Institute of Technology Voice: 312.567.3498 Fax: 312.567.3494 Carlo.Segre(a)iit.edu http://www.iit.edu/~segre

1 0

new windows updater
by Bruce Ravel 14 Oct '05

14 Oct '05

Hi again, The windows updater containing an update to the runner as well as all the new stuff from today's Athena&Artemis source code release is now available. The runner should find it the next time you start one of the programs. Let us know if you have any troubles with the update or with the new runner. Bruce & Matt -- Bruce Ravel ---------------------------------------------- bravel(a)anl.gov Molecular Environmental Science Group, Building 203, Room E-165 MRCAT, Sector 10, Advance Photon Source, Building 433, Room B007 Argonne National Laboratory phone and voice mail: (1) 630 252 5033 Argonne IL 60439, USA fax: (1) 630 252 9793 My homepage: http://cars9.uchicago.edu/~ravel EXAFS software: http://cars9.uchicago.edu/~ravel/software/

1 0

tarball 056 released & new tip of the week
by Bruce Ravel 14 Oct '05

14 Oct '05

Hi folks! The new Tip of the Week is about automated replotting in Athena. http://cars9.uchicago.edu/iffwiki/HoraeTipOfTheWeek I have also released version 056 of the source code tarball, which is available at my web site and at SourceForge. A description of what's new in Athena and Artemis is at http://cars9.uchicago.edu/iffwiki/HoraeVersion056 The major new features are both in Athena: session defaults and improvements to linear combination fitting. Please read the release notes for a complete description of these new features. An updater for windows users should be available within a few days. Unix and Mac OSX users should download and install the new source tarball. Enjoy, B -- Bruce Ravel ---------------------------------------------- bravel(a)anl.gov Molecular Environmental Science Group, Building 203, Room E-165 MRCAT, Sector 10, Advance Photon Source, Building 433, Room B007 Argonne National Laboratory phone and voice mail: (1) 630 252 5033 Argonne IL 60439, USA fax: (1) 630 252 9793 My homepage: http://cars9.uchicago.edu/~ravel EXAFS software: http://cars9.uchicago.edu/~ravel/software/exafs/

1 0

Small clusters in Feff
by Kelly, Shelly D. 07 Oct '05

07 Oct '05

Hi EXAFS friends, I was reading a paper this morning where the authors used a small hydrated cluster in Feff8.2. The cluster has a metal center surrounded by one shell of oxygen atoms. All feff paths in the calculation were used in the model fits. No mention was made in the manuscript of the muffin tin radii used in the Feff8.2 calculation based on this small cluster, hence I assume that these authors were unaware of the potential pit falls in doing so. Using feff7, I have run a comparison of the calculated signal from a bounded and unbounded shell of atoms at the same distance. I found that the unbounded shell of atoms has a muffin tin radii that is really big because it is not bounded by any other atoms giving rise to an error in the calculated signal from these atoms. Has anyone run this experiment with Feff8.2? I have always used the metal-oxide crystal structure in my modeling of hydrated metal ions and have found reasonable agreement with chemical speciation calculations. I believe that the non-expert users of Feff, find it difficult to remember not to use the exact model within Feff because of this "feature". Is there any discussion of adding a hydration sphere option into Feff? I look forward to your comments. Shelly Shelly Kelly Bldg 203 RM E113 Skelly(a)anl.gov Argonne National Laboratory 630-252-7376 9700 S Cass Ave www.mesg.anl.gov Argonne, IL 60440

2 1

new tip of the week
by Bruce Ravel 07 Oct '05

07 Oct '05

Hi, The Tip Of The Week took last week off to go visit my mom ;-) This week's is about "Selecting paths to plot automatically after a fit" in Artemis. http://cars9.uchicago.edu/iffwiki/HoraeTipOfTheWeek Enjoy! B -- Bruce Ravel ---------------------------------------------- bravel(a)anl.gov Molecular Environmental Science Group, Building 203, Room E-165 MRCAT, Sector 10, Advance Photon Source, Building 433, Room B007 Argonne National Laboratory phone and voice mail: (1) 630 252 5033 Argonne IL 60439, USA fax: (1) 630 252 9793 My homepage: http://cars9.uchicago.edu/~ravel EXAFS software: http://cars9.uchicago.edu/~ravel/software/exafs/

1 0

delR versus thermal expansion A and first cumulant C1*
by Leandro Langie Araujo 06 Oct '05

06 Oct '05

Dear all, I’m writing because I’ve been trying to analyze some thermal dependent EXAFS data for the first shell of a c-Ge sample and ended up in a conceptual confusion. If any of you could give me some help on this I’d be very grateful. Whenever we fit an EXAFS spectrum with Artemis/IFEFFIT, we get a /delR/ value. As far as my (very) limited comprehension goes, /delR/ is the difference between the theoretical (FEFF8) bond length, /R0/, and the result of the best fit to my data, /R/, which gives the experimental bond length. So, this should be the /first cumulant/ of my distance distribution, right? The net linear thermal expansion /A/ (as defined by Anatoly Frenkel and John Rehr in [PRB48, 585, 1993]) from an one-dimensional model seems to have been considered as equal to the temperature variation of the /delR/ values by several people, including Matt Newville (on his thesis). Following this approach, it is possible to relate the EXAFS cumulants (/delR/, /ss/, /C3/ and /C4/) of the distance distribution with the constants of an one-dimensional anharmonic effective potential, following a correlated Einstein model. I was going through this path... … when I came across some papers (P. Fornasini and G. Dalba in [PRL82, 4240, 1999; PRB70, 174301, 2004] and E. A. Stern in [J.Phys.IV 7, C2-137,1997]) where an issue is raised; it is said that the thermal evolution of the first cumulant /C1/* of the real distance distribution in a crystal (which was said to be equivalent to /delR/) is NOT equal to /A/, because in a crystal there are vibrations perpendicular to the bond direction which are not considered in the one-dimensional model. It is argued that the temperature dependence of these vibrations perpendicular to the bond direction give raise to a positive shift of the minimum of the effective pair potential, while the net thermal expansion /A/ accounts only for changes due to the asymmetry of the potential. So, the main questions hammering my head are: - does /delR/ include any contribution from vibrations perpendicular to the bond direction? Can really /delR/ and /C1/* be the same quantity? - is the thermal evolution of /delR/ equal to the net thermal expansion /A/? Or should some correction for perpendicular vibrations be added to relate both quantities? Sorry if I made it too long, but it was something hard to express in a few words… Any help will be greatly appreciated. Cheers, Leandro

3 4

R-factor and reduced chi square
by OWEN N LI 05 Oct '05

05 Oct '05

Hi, I understand R-factor and reduced chi square are statistical ways to see if the model is reasonable (R-factor) and if one model is better than the others (reduced chi square). But how much should I read into it? I have model A and B, I believe A is more theoretical sound than B. but when I fit them using IFEFFIT, the R-factor and reduced chi square of B is slightly better than A (B: 0.03, 28; A: 0.04, 33). What does that tell me? Is my assumption incorrect? Or are they more or less the same? Owen

5 4

RE: [Ifeffit] New Method for Normalization?
by Barton, David (DG) 04 Oct '05

04 Oct '05

Bruce, You mentioned: > Isn't Ifeffit's Cromer-Liberman normalization scheme > substantially similar to the method presented in that paper? > When I use the Cromer-Liberman background subtraction routine in Athena the Normalized data is not flattened. Normalized, but not *flattened* data results in a skewed XANES region which is affected by the amount of curvature in the background. If I choose "Flatten" in Athena, it is my understanding that it applies a linear pre-edge and 2nd order polynomial post-edge regardless of which background was chosen (Autobk or CLNorm). The Clnorm background does appear to "flatten" the chi(k) data which can be useful for EXAFS analysis, but it doesn't flatten the mu(E) data. When one attempts to analyze XANES data it is preferred to have the entire post-edge "flattened" in mu(E) space so that the background curvature is eliminated from the dataset. Ultimately, as Matt stated, it would be nice to be able to combine a Victoreen and CLNorm background to get a reasonably "flattened" mu(E) function. I would also suggest adding the ability to independently set the pre-edge and post-edge functions (n'th order polynomial, Victoreen, CLNorm, or MBACK??) for Normalization and Flattening. These choices would be separate from the choice of background removal method to get chi(k) (Autobk or CLNorm). Dave > -----Original Message----- > From: ifeffit-bounces(a)millenia.cars.aps.anl.gov > [mailto:ifeffit-bounces@millenia.cars.aps.anl.gov] On Behalf > Of Bruce Ravel > Sent: Thursday, September 29, 2005 5:20 PM > To: XAFS Analysis using Ifeffit > Subject: Re: [Ifeffit] New Method for Normalization? > > > On Thursday 29 September 2005 15:48, Barton, David (DG) wrote: > > Has anyone considered using the methods of Penner-Hahn to normalize > > their data to get reliable intensities of near-edge features? The > > reference to their method of normalization is below: > > > > J. Synchrotron Rad. (2005). 12, 506-510 > [doi:10.1107/S0909049504034193] > > "A method for normalization of X-ray absorption spectra" > T.-C. Weng, > > G. S. Waldo and J. E. Penner-Hahn > > > > It is my understanding the Ifeffit uses a linear function to remove > > the pre-edge and a quadratic for the post-edge and that the > Athena's > > flatten is a subtraction of these functions. This is an excellent > > method under most scenarios since the background is usually fairly > > smooth and can be well approximated by a quadratic equation. > > Personally, I have had only a few rare cases where the quadratic > > function on the post-edge was not sufficient to reliably > normalize the > > data and in those cases a third-order polynomial was > sufficient. Does > > anyone else have an opinion on using alternative normalization > > routines? > > Isn't Ifeffit's Cromer-Liberman normalization scheme > substantially similar to the method presented in that paper? > > B > > > -- > Bruce Ravel ---------------------------------------------- > bravel(a)anl.gov > > Molecular Environmental Science Group, Building 203, Room > E-165 MRCAT, Sector 10, Advance Photon Source, Building 433, > Room B007 > > Argonne National Laboratory phone and voice mail: > (1) 630 252 5033 > Argonne IL 60439, USA fax: > (1) 630 252 9793 > > My homepage: http://cars9.uchicago.edu/~ravel > EXAFS software: http://cars9.uchicago.edu/~ravel/software/exafs/ > > > _______________________________________________ > Ifeffit mailing list > Ifeffit(a)millenia.cars.aps.anl.gov > http://millenia.cars.aps.anl.gov/mailman/listinfo/ifeffit >

1 0

RE: [Ifeffit] R-factor and reduced chi square
by Kelly, Shelly D. 04 Oct '05

04 Oct '05

Hi Owen, A test for a significant change in the reduced chi square (RCS). [(RCS1)^2/(RCS2)^2] - 1 greater than 2 * sqrt ( 2 / degrees of freedom) If you can not read that, I have this published in J. Phys Chem B, V105 no27 year (2001). You can down load that paper from my web page www.mesg.anl.gov. I guess it is a reasonable test:). For most model fits a change by a factor of 2 is needed to be statistically significant. Good luck Shelly Shelly Kelly Bldg 203 RM E113 Skelly(a)anl.gov Argonne National Laboratory 630-252-7376 9700 S Cass Ave www.mesg.anl.gov Argonne, IL 60440 > -----Original Message----- > From: ifeffit-bounces(a)millenia.cars.aps.anl.gov [mailto:ifeffit- > bounces(a)millenia.cars.aps.anl.gov] On Behalf Of OWEN N LI > Sent: Tuesday, October 04, 2005 1:08 PM > To: ifeffit(a)millenia.cars.aps.anl.gov > Subject: [Ifeffit] R-factor and reduced chi square > > Hi, > I understand R-factor and reduced chi square are statistical ways to > see if the model is reasonable (R-factor) and if one model is better > than the others (reduced chi square). But how much should I read into it? > > I have model A and B, I believe A is more theoretical sound than B. > but when I fit them using IFEFFIT, the R-factor and reduced chi square > of B is slightly better than A (B: 0.03, 28; A: 0.04, 33). What does > that tell me? Is my assumption incorrect? Or are they more or less the > same? > > Owen > > > _______________________________________________ > Ifeffit mailing list > Ifeffit(a)millenia.cars.aps.anl.gov > http://millenia.cars.aps.anl.gov/mailman/listinfo/ifeffit

1 0

Muffin tin radius adjustment?
by Paul Fons 04 Oct '05

04 Oct '05

I am trying to analyze some XAFS data on a Sb-Te compound and when I run feff(6) I get errors regarding the coordination number being too small implying the muffin tin radius is too big. In this material the first nearest neighbor bond length is about 3 angstroms (slightly more or less depending on the site). I assume this is related to the cause of the problem. Any suggestions as to how to fix/avoid this sort of problem? feff6 output is pasted below. Interestingly enough, feff8 does not output any error messages and the sum of the paths looks identical (at least at first glance). Is this a false error message or (more likely) has feff8 become more intelligent with regards to defining the muffin tin radius? Feff6 output Calculating potentials and phases... free atom potential and density for atom type 0 free atom potential and density for atom type 1 free atom potential and density for atom type 2 overlapped potential and density for unique potential 0 overlapped potential and density for unique potential 1 overlapped potential and density for unique potential 2 muffin tin radii and interstitial parameters phase shifts for unique potential 0 Hard tests failed in fovrg. Muffin-tin radius may be too large; coordination number too small. phase shifts for unique potential 1 Hard tests failed in fovrg. Muffin-tin radius may be too large; coordination number too small. phase shifts for unique potential 2 Hard tests failed in fovrg. Muffin-tin radius may be too large; coordination number too small. Preparing plane wave scattering amplitudes... nncrit in prcrit 9 Searching for paths... WARNING in PATHS Module: rmax > distance to most distant atom. Some paths may be missing. rmax, ratx 15.50000 15.40291 Rmax 15.5000 keep and heap limits 0.0000000 0.0000000 Preparing neighbor table Feff8 output Calculating potentials ... free atom potential and density for atom type 0 free atom potential and density for atom type 1 free atom potential and density for atom type 2 initial state energy overlapped potential and density for unique potential 0 overlapped potential and density for unique potential 1 overlapped potential and density for unique potential 2 muffin tin radii and interstitial parameters iph, rnrm(iph)*bohr, rmt(iph)*bohr, folp(iph) 0 1.79879E+00 1.71274E+00 1.13279E+00 1 1.81342E+00 1.72495E+00 1.13595E+00 2 1.83703E+00 1.74437E+00 1.14149E+00 mu_old= -0.508 Done with module 1: potentials. Calculating LDOS ... LDOS calculation for specified grid Calculating energy and space dependent l-DOS. It takes time ... potential type 0 potential type 1 potential type 2 Calculating chi and rho... 0 Calculating chi and rho... 1 Calculating chi and rho... 2 Done with LDOS. Calculating cross-section and phases... absorption cross section phase shifts for unique potential 0 phase shifts for unique potential 1 phase shifts for unique potential 2 Done with module 2: cross-section and phases... Done with module 3: FMS. Preparing plane wave scattering amplitudes... Searching for paths... WARNING: rmax > distance to most distant atom. Some paths may be missing. rmax, ratx 1.55000E+01 0.00000E+00 Rmax 15.5000 keep and heap limits 0.0000000 0.0000000 Preparing neighbor table 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(a)aist.go.jp The lines below are in a Japanese font 〒３０５−８５６８茨城県つくば市東１−１−1 つくば中央第４近接場光応用工学センターポール・フォンス主任研究官

3 2