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

How to generate the input information for Atom
by Fenglong Sun 06 Nov '07

06 Nov '07

Hi, Bruce, I am working on an experiment of arsenic adsorption on pyrite particles. I have the data of mu(E)~E. I use Anthena to do the Fourier transform. Now I want to know the cordination information of arsenic atom. How can I get the input file for Atom if I want to do the fitting job by using Iffefit package? Thanks!

2 1

Build lower doping concentration structure in TKAtoms program.
by 沛昌蔡 06 Nov '07

06 Nov '07

Dear Prof: We have some questions about build lower doping concentration structure in TKAtoms program, as following; Q1: The TKATOM program can build up doping structure for lower doping concentration (below 10%) or not? Q2: If it can work, how should we build up the doping site? Q3: For lower doping case, how do these doped atoms choose where to go in TKATOM program, is it have physical description? ps: The attach ppt. file show more detail about our questions. Thanks for helping! Pei-chang Tsai ____________________________________________________________________________________ 怎樣減少垃圾信？只要看到垃圾信，立即按下「這是垃圾信」按鈕，垃圾信剋星幫助你清除你信箱垃圾。

2 1

help: can I use ifeffit to predict chi?
by Tai-Yen Chen 02 Nov '07

02 Nov '07

hi, Can anyone tell me whether I can use ifeffit to do this job? I want to know the % change in chi when the molecule have 3% change in radius? All I have is just the crystal struture of gamma-Fe2O3, which means I can only get the Feff. How can I introduce 3% change in redius so that I can get the difference of chi? Can anyone help me? Tanks a lot. TaiYen Chen 979-739-7772 Department of Chemistry Texas A&M University P.O. Box 30012 College Station, TX 77842-3012 Visit us on the web at http://www.chem.tamu.edu

4 4

Re: [Ifeffit] 1. Re xanes feff8 sulfur pyrite
by s440697＠stud.uni-goettingen.de 26 Oct '07

26 Oct '07

Hi Josh, i have made two calculations respectively pyrite. You can see the results in the postscript file. The first (A) plot is from your pyrite.dat, the second plot is based on your feff.inp but recalculated on my Computer. The last (C) plot is based on my own crystallographic data (spacegroup p 2_1 /a -3) and feff.inp (attached as eck-feff.inp) file, maybe the same crystallographic data you used? With my feff.inp i can approve all of your atom-positions with the exception of atom 11 in your list respectively atom Fe3 in my list. But the results in the diagramm are nearly the same. The third file (feff.inp) is again your feff.inp file you sent to me. >Hi Josh, > >thanks for your help and attached files. I will test your feff.inp files >and compare the results. > >cheers, >Eckhard > > Hi Eckhard, > Try the attatched feff.inp files (gypsum.inp and pyrite.inp). > I find that these produce different spectra with edges at > ~2468.5 - gypsum > ~2469.0 - pyrite > > If you use do not get the same results, then there is something wrong with the code and we will have some debugging to do. > > Also, if there are multiple sulfur sites in the material which have different local structure, you will > need to do a configurational average by running with each site as the absorbing atom and averaging the signals. > Cheers, > Josh > > P.S. > I am also sending the result (gypsum.dat and pyrite.dat) so that you can compare. > > > On Tue, 16 Oct 2007, ifeffit-request(a)millenia.cars.aps.anl.gov wrote: > >> Send Ifeffit mailing list submissions to >> ifeffit(a)millenia.cars.aps.anl.gov >> >> To subscribe or unsubscribe via the World Wide Web, visit >> http://millenia.cars.aps.anl.gov/mailman/listinfo/ifeffit >> or, via email, send a message with subject or body 'help' to >> ifeffit-request(a)millenia.cars.aps.anl.gov >> >> You can reach the person managing the list at >> ifeffit-owner(a)millenia.cars.aps.anl.gov >> >> When replying, please edit your Subject line so it is more specific than "Re: Contents of Ifeffit digest..." >> >> >> Today's Topics: >> >> 1. Re xanes feff8 sulfur (Kelly, Shelly D.) >> 2. Re: Re xanes feff8 sulfur (s440697(a)stud.uni-goettingen.de) >> >> >> ---------------------------------------------------------------------- >> >> Message: 1 >> Date: Mon, 15 Oct 2007 15:18:14 -0500 >> From: "Kelly, Shelly D." <SKelly(a)anl.gov> >> Subject: [Ifeffit] Re xanes feff8 sulfur >> To: "XAFS Analysis using Ifeffit" <ifeffit(a)millenia.cars.aps.anl.gov> Cc: "John J. Rehr" <jjr(a)leonardo.phys.washington.edu>, John Katsoudas >> <katsoudas(a)agni.phys.iit.edu> >> Message-ID: >> <0BC5625B1BD9494886EBA2F3CF10E957012C1C90(a)BIOPOST.bio.anl.gov> >> Content-Type: text/plain; charset="us-ascii" >> >> Hi Eckard, >> >> >> >> I realize that you have posted a message every month for about 3 months in a row without any reply: Shame on us... I will try to answer some of your questions and maybe that will help to get some discussion going. The energy scale in feff8 is NOT defined as you would expect it to be for determining the absolute energy of an edge. The energy scale in FEFF8 is related to the difference between the energy level of the intersist (potential level inbetween the muffin tins that represent atoms) and the threshold energy level for the excitation of the electron. >> >> >> >> I realize that this doesn't solve your problem, but hopefully it will spark some discussion. >> >> >> >> Cheers, >> >> Shelly >> >> >> >> >> Mon Aug 13 06:00:19 CDT 2007 >> Does anybody know a person, who has experience respective the topics xanes+feff8 and sulfur K-Edges? >> I would be very thankful for an answer :-) >> >> sincerely yours, >> Eckhard Bosman >> >> >> >> >> Tue Jul 10 10:07:38 CDT 2007 >> >> Hello, >> >> i have problems with feff-xanes calculations related to pyrite. I try to get the k-edge of sulfur (as absorber) in FeS2. I get the edge at ~2472 eV >> and i think this is wrong because i get the same energy with gypsum and i >> think there should be an energyshift related to 2472 eV. My proceeding is >> as follows: >> >> calculate the atomic positions with the asymetric unit from the http://rruff.geo.arizona.edu/AMS/amcsd.php database (Pyrite spacegroup Pa3) respectively >> http://rruff.geo.arizona.edu/AMS/xtal_data/CIFfiles/00596.cif >> with "powdercell 2.4" >> >> feffrun with the feff.inp-file like in the attachment >> (thanks to Josh Kas for his feff-hints) >> >> I started the feff-run on the one hand with 12 atoms (unitcell) and on the >> other hand with 22 atoms (14x Fe and 8x S in the cube) >> >> After this i plotted from xmu.dat the absorbtion mu (column 4) as a function of the energy (column 1). The result is in the attachment too. >> >> I don't know why i get everytime the Sulfur-k-edge at the same position (for pyrite or gypsum) respectively without energyshift related to 2472 eV. >> >> many thanks for your help :-) >> Ciao, >> Eckhard >> >> >> >> Fri Jun 15 08:48:21 CDT 2007 >> >> Hallo @all, >> >> (first at all, sorry for my bad english) >> >> i'm a student working on my diploma thesis (nexafs/xanes) and my job is to >> simulate with FEFF 8.4 spectra, for example absorption as a function of energy (mu_(E)). I'm a newbie in this topic and i tried to get the k-edge >> of Sulfur for the substance gypsum (CaSO4 * 2(H2O)) and pyrite (FeS2). I'm >> interested in the XANES-range of the spectra not the EXAFS. I get everytime my peaks at ~2472 eV for Sulfur as absorber whatever substance i >> use (gypsum, pyrite,...). There are no differences with the peak position >> on the energy-scale or shifts. I'm really confused :-( I think, there should be differences in the spectra... I send you two feff.inp files from >> my feff-runs. My source is the mineralogy database >> (http://rruff.geo.arizona.edu/AMS/amcsd.php) where i get the data for the >> unitcell. Then i put the data in the XtalDraw program to get the positions >> of the Atoms in Angstrom, by the way i get i nice 3D-View of the unitcell >> :-). Then i put the values of the positions into feff... >> I simulated pyrite with 12 atoms and gypsum with 48. >> I'm appreciative for every help :-) >> >> ciao, >> Eckhard >> >> >> >>

1 0

scattering from H atoms
by Yael Politi 24 Oct '07

24 Oct '07

Hallo everyone, I have data of an hydrated mineral (a standard material, with known crystallographic structure), where the nearest neighbors are hydrogen atoms at about 0.8 angs from the absorber. my question, if anybody has had similar problem, is how do I differentiate between noise and contribution from H scattering in my exafs data. If I want to, I can get to a stage where I have almost no intensify in the FT data at r<1. but I don't know if this is a better interpretation. the feff input file suggest strong scattering from H atoms and I wouldn't want to just ignore these paths in the fit. Any suggestions would be most appreciated. Yael

1 0

1. Re xanes feff8 sulfur
by joshua jason kas 21 Oct '07

21 Oct '07

Hi Eckhard, Try the attatched feff.inp files (gypsum.inp and pyrite.inp). I find that these produce different spectra with edges at ~2468.5 - gypsum ~2469.0 - pyrite If you use do not get the same results, then there is something wrong with the code and we will have some debugging to do. Also, if there are multiple sulfur sites in the material which have different local structure, you will need to do a configurational average by running with each site as the absorbing atom and averaging the signals. Cheers, Josh P.S. I am also sending the result (gypsum.dat and pyrite.dat) so that you can compare. On Tue, 16 Oct 2007, ifeffit-request(a)millenia.cars.aps.anl.gov wrote: > Send Ifeffit mailing list submissions to > ifeffit(a)millenia.cars.aps.anl.gov > > To subscribe or unsubscribe via the World Wide Web, visit > http://millenia.cars.aps.anl.gov/mailman/listinfo/ifeffit > or, via email, send a message with subject or body 'help' to > ifeffit-request(a)millenia.cars.aps.anl.gov > > You can reach the person managing the list at > ifeffit-owner(a)millenia.cars.aps.anl.gov > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Ifeffit digest..." > > > Today's Topics: > > 1. Re xanes feff8 sulfur (Kelly, Shelly D.) > 2. Re: Re xanes feff8 sulfur (s440697(a)stud.uni-goettingen.de) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 15 Oct 2007 15:18:14 -0500 > From: "Kelly, Shelly D." <SKelly(a)anl.gov> > Subject: [Ifeffit] Re xanes feff8 sulfur > To: "XAFS Analysis using Ifeffit" <ifeffit(a)millenia.cars.aps.anl.gov> > Cc: "John J. Rehr" <jjr(a)leonardo.phys.washington.edu>, John Katsoudas > <katsoudas(a)agni.phys.iit.edu> > Message-ID: > <0BC5625B1BD9494886EBA2F3CF10E957012C1C90(a)BIOPOST.bio.anl.gov> > Content-Type: text/plain; charset="us-ascii" > > Hi Eckard, > > > > I realize that you have posted a message every month for about 3 months > in a row without any reply: Shame on us... I will try to answer some of > your questions and maybe that will help to get some discussion going. > The energy scale in feff8 is NOT defined as you would expect it to be > for determining the absolute energy of an edge. The energy scale in > FEFF8 is related to the difference between the energy level of the > intersist (potential level inbetween the muffin tins that represent > atoms) and the threshold energy level for the excitation of the > electron. > > > > I realize that this doesn't solve your problem, but hopefully it will > spark some discussion. > > > > Cheers, > > Shelly > > > > > Mon Aug 13 06:00:19 CDT 2007 > Does anybody know a person, who has experience respective the topics > xanes+feff8 and sulfur K-Edges? > I would be very thankful for an answer :-) > > sincerely yours, > Eckhard Bosman > > > > > Tue Jul 10 10:07:38 CDT 2007 > > Hello, > > i have problems with feff-xanes calculations related to pyrite. I try to > get the k-edge of sulfur (as absorber) in FeS2. I get the edge at ~2472 > eV > and i think this is wrong because i get the same energy with gypsum and > i > think there should be an energyshift related to 2472 eV. My proceeding > is > as follows: > > calculate the atomic positions with the asymetric unit from the > http://rruff.geo.arizona.edu/AMS/amcsd.php database (Pyrite spacegroup > Pa3) respectively > http://rruff.geo.arizona.edu/AMS/xtal_data/CIFfiles/00596.cif > with "powdercell 2.4" > > feffrun with the feff.inp-file like in the attachment > (thanks to Josh Kas for his feff-hints) > > I started the feff-run on the one hand with 12 atoms (unitcell) and on > the > other hand with 22 atoms (14x Fe and 8x S in the cube) > > After this i plotted from xmu.dat the absorbtion mu (column 4) as a > function of the energy (column 1). The result is in the attachment too. > > I don't know why i get everytime the Sulfur-k-edge at the same position > (for pyrite or gypsum) respectively without energyshift related to 2472 > eV. > > many thanks for your help :-) > Ciao, > Eckhard > > > > Fri Jun 15 08:48:21 CDT 2007 > > Hallo @all, > > (first at all, sorry for my bad english) > > i'm a student working on my diploma thesis (nexafs/xanes) and my job is > to > simulate with FEFF 8.4 spectra, for example absorption as a function of > energy (mu_(E)). I'm a newbie in this topic and i tried to get the > k-edge > of Sulfur for the substance gypsum (CaSO4 * 2(H2O)) and pyrite (FeS2). > I'm > interested in the XANES-range of the spectra not the EXAFS. I get > everytime my peaks at ~2472 eV for Sulfur as absorber whatever substance > i > use (gypsum, pyrite,...). There are no differences with the peak > position > on the energy-scale or shifts. I'm really confused :-( I think, there > should be differences in the spectra... I send you two feff.inp files > from > my feff-runs. My source is the mineralogy database > (http://rruff.geo.arizona.edu/AMS/amcsd.php) where i get the data for > the > unitcell. Then i put the data in the XtalDraw program to get the > positions > of the Atoms in Angstrom, by the way i get i nice 3D-View of the > unitcell > :-). Then i put the values of the positions into feff... > I simulated pyrite with 12 atoms and gypsum with 48. > I'm appreciative for every help :-) > > ciao, > Eckhard > > > >

2 1

Fe-Co help
by Dominik Samuelis 21 Oct '07

21 Oct '07

Dear Chachi, first of all, you need to work on your crystallographic description of the system. I'm sure you are aware of the ordering transistion in Fe-Co at around 730 °C. Your description currently uses parts of the ordered phase, and other parts of the unordered. Above 730 °C, your description as a bcc structure is correct, however your atomic parameters and the neighbor shells are not. The system then has a bcc structure (Im-3m), with mixed Fe/Co occupancy on the single crystallographic position 2a. It should be modelled as such, with mixed occupancy on the neighboring sites. There's excellent tutorials how to deal with mixed sites. However, be aware that the scattering properties of iron and cobalt is very similar, you you won't really be able to fit your site occupancy. Below 730 °C, the system starts to order. The ordered phase is NOT bcc anymore, but consists of two intersecting primitive cubic lattices, shifted by (1/2 1/2 1/2). The structure prototype is often called B2, after its Strukturberichte symbol (http://cst- www.nrl.navy.mil/lattice/struk/index.html) Now, the rest of your crystallographic description matches. The neighbor shells will be N1=8Co and N2=6Fe. However, the symmetry is not bcc (Im-3m), but primitive (Pm-3m). Looking at your fit and feff calculation, it seems, that you have used the fully ordered phase for modelling your alloy. As you have pointed out yourself, it's not clear that the alloy really is fully ordered. Your alloy might be partially ordered, as also indicated by the Mossbauer data. In this case, you need to mix several feff calculations, with each cobalt and iron on the first and second neighbor position. You could start by assuming an ordering parameter 'op'. The first shell would then be N1=8(op Co+(1-op)Fe), the second is N2=6(op Fe+(1-op)Co). Cheers, dominik > Dear all > I?m trying to deal with exafs data of a Fe-Co alloy sample, using > Artemis. > The bcc structure is simple to modeling, Fe at 0,0,0 and Co at > 1/2,1/2,1/2 > (at Fe-edge, N1= 8Co, N2= 6Fe, theorically), > But in my fits N1=5.4 and N2=4.9, that is correct? > amp=0.811293 value I used from Iron calibration data, that is > correct too? > A possibility (I think), is that the sample is inhomogeneous, > because, by Mossbauer > spectroscopy, show two Fe magnetic sites, one attributed to ordered > phase, and the > other to disordered phase. -- Dipl.-Ing. Dominik Samuelis samuelis(a)pc.rwth-aachen.de Physikalische Chemie I RWTH Aachen University Phone: +49-241-80-94840 Fax: +49-241-80-92128 Landoltweg 2 52056 Aachen Germany

1 0

negative e0 value
by mauricio.maps 20 Oct '07

20 Oct '07

Dear all, I´m trying to deal with exafs data of a Ni-phosphate compound, using Athena/Artemis. This is a simple structure (I think), with the hexacoordinated ion surrounded by water molecules on the first shell. The question I have concerns the negative value for e0 I´m having. Is this a bad thing? Maybe I´m doing something wrong? I´ve tried some issues, as those proposed on the Shelly Kelly´s tutorial, but without success. Can you help me, please? Attached are an example of the Athena and Artemis projects i´m using. Thanks a lot. Mauricio

3 3

Debye model and Debye Waller factor
by Cammelli Sebastiano 19 Oct '07

19 Oct '07

Dear Bruce In your "Introduction to FEFFIT" (http://feff.phys.washington.edu/~ravel/) you suggest in the beginning to use the DW factor according to the Debye model [debye(temp, thetad)] for the simple Cu foil. Then you ask to try with the Einstein model. I gave a look at the theory but I did not understand how to apply it to the fit. Anyway performing some FEFF calculation of this Cu foil, it looks that the Debye model works quite well but I am a very beginner feffusers and I am not an expert about IFFEFIT too. Anyway the samples I have to analyze are metallic alloy mainly based on Fe with some diluted element such as (Cu, Ni, Mn each less than 1at%). Is the Debye model good approximation for transition metal alloy? Is it better just to leave all Debye Waller factors free? In FEFF calculation, does it have any meaning to use the NOHOLE card and the S02=0.9? Thank you for your help Best regards, Sebastiano Cammelli

2 1

Re: [Ifeffit] [FEFF8.4] Re(x) is .le. zero in besjn
by joshua jason kas 18 Oct '07

18 Oct '07

Hi Mauro, This is probably a bug in feff. There is probably a workaround. Can you send your feff.inp file and I will try to help find the right cards to make your file run. Also I will add this to a list of bugs to look at. Cheers, Josh > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 01 Oct 2007 18:17:50 +0200 > From: Mauro Rovezzi <rovezzi(a)esrf.fr> > Subject: [Ifeffit] [FEFF8.4] Re(x) is .le. zero in besjn ? > To: Ifeffit-ML <ifeffit(a)millenia.cars.aps.anl.gov> > Message-ID: <47011DAE.2010508(a)esrf.fr> > Content-Type: text/plain; charset=ISO-8859-1 > > Hi all, > > A little question about a FEFF8 error message I have and I cannot > understand what it means : > > Re(x) is .le. zero in besjn > > This error is generated when I try to add the SCF to an EXAFS > calculation on MnO (manganosite) and the SCF iteration is at step 9 of 30. > > Seen that the cards in the feff.inp file are properly used (I normally > run SCF in EXAFS calculations) do you have any idea of what is going on? > > Thanks for your help, as usual! > > Mauro > > -- > Mauro Rovezzi <rovezzi(a)esrf.fr> > PhD student - GILDA BM08 c/o ESRF (Grenoble, FR) > > > ------------------------------ > > _______________________________________________ > Ifeffit mailing list > Ifeffit(a)millenia.cars.aps.anl.gov > http://millenia.cars.aps.anl.gov/mailman/listinfo/ifeffit > > > End of Ifeffit Digest, Vol 56, Issue 1 > ************************************** >

2 1