Hi Zohair, I'm posting this back to the list, because I couldn't solve the problem. I did experiment with your file for a bit; it works fine for the atoms up to 4 angstroms, but ones much past that make feff fail. Maybe someone else can suggest why feff is choking on what appears on inspection to be a valid file. (Note to others on the list: the attached feff.inp file is the one that doesn't work; the pasted one is my original, and does work.) As for atoms being able to work with icosahedral structures, that will probably never be possible. Atoms is designed to work with crystal structures, and icosahedral clusters are not crystals; i.e. they don't have translational symmetry. A separate helper application would have to be written for this kind of cluster. --Scott Calvin Sarah Lawrence College Begin forwarded message:
From:
Date: June 27, 2011 1:12:06 AM PDT To: Subject: Re: [Ifeffit] Icosahedral cluster coordinates Dear Scott,
I came across an e-mail exchange from 2005 between you and Patrick Kluth regarding icosahedral clusters and feff files as below.
I’m trying to fit data collected at the Cu k-edge for Mg65Cu25Gd10 bulk metallic glass(BMG). It has been reported that BMGs have an icosahedral structure. A recent paper (attached) gives a demonstration of that. What I’ve tried to do is to modify the feff file you made for Patrick to my MgCuGd structure. I attach the feff.inp file I made for that.
The problem I’m running into is that when I try and run this modified feff file through artemis I get an error message (file attached). I should point out that when I run your feff file for Fe it runs without any errors.
I’m wondering if you can either tell me what I’m doing wrong or guide me to make a feff file from an atoms input file for the MgCuGd material.
In the long run, it would be good to know how to make Feff files for icosahedral structures starting from an atoms input file (I usually use webatoms for this) as I would be working with various BMG materials (Al, Ti and Mg based) that should potentially have this structure.
I’d be grateful for your help.
Kind Regards,
Zohair Hussain Postdoctoral Fellow
CPSE - CSIRO Private Bag 33 Clayton South MDC Clayton VIC 3169 Australia e-mail: zohair.hussain@csiro.au Phone +61 3 9545 7891 Fax +61 3 9544 1128
[Ifeffit] Icosahedral cluster coordinates
Scott Calvin SCalvin at slc.edu Tue Jul 19 19:06:17 CDT 2005 Previous message: [Ifeffit] Icosahedral cluster coordinates Next message: [Ifeffit] Icosahedral cluster coordinates Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] Hi Patrick,
Hopefully this will do it. A disclaimer: I never actually tested this on an icosahedral material, and I did the coordinates by hand, so I can't guarantee it.
--Scott Calvin Sarah Lawrence College
I am working with metal clusters and would like to do some FEFF calculations on icosahedral clusters. Is there anyone out there who has a FEFF input file (or equivalently a list of coordinates) for an icosahedral cluster (147 atoms is probably sufficient for my purpuses).
Perhaps, some problems with constructing atomic coordinates? It appears that they were obtained by modifying the hcp structure... but I am not sure if that algorithm listed in the input file was correct (from the input file: "*c and a have been chosen to give 5% expansion between layers and 10% within the layer in the first shell as per Farges (1982).This thus roughly simulates a relaxed icosohedron").
In a regular icosahedron all distances in the radial direction are the same, and all distances within the shell are about 5% longer. This program can be used to generate atomic coordinates in an icosahedron, up to the 10th order (for larger icosahedral clusters you may contact me directly).
http://www3.bnl.gov/frenkel/coords.html
FEFF works on these coordinates (it has been tested).
Regards,
Anatoly
________________________________
From: ifeffit-bounces@millenia.cars.aps.anl.gov on behalf of Scott Calvin
Sent: Mon 6/27/2011 1:06 PM
To: XAFS Analysis using Ifeffit
Subject: [Ifeffit] Fwd: Icosahedral cluster coordinates
Hi Zohair,
I'm posting this back to the list, because I couldn't solve the problem. I did experiment with your file for a bit; it works fine for the atoms up to 4 angstroms, but ones much past that make feff fail. Maybe someone else can suggest why feff is choking on what appears on inspection to be a valid file. (Note to others on the list: the attached feff.inp file is the one that doesn't work; the pasted one is my original, and does work.)
As for atoms being able to work with icosahedral structures, that will probably never be possible. Atoms is designed to work with crystal structures, and icosahedral clusters are not crystals; i.e. they don't have translational symmetry. A separate helper application would have to be written for this kind of cluster.
--Scott Calvin
Sarah Lawrence College
Begin forwarded message:
From:
Wouldn't you just put the atoms at the vertices of an icosahedron? That is, how would you know what was beyond the first shell? The feff.inp file below might be a good start. --Matt TITLE Fe at center of an icosahedron HOLE 1 1.0 Fe K edge CONTROL 1 1 1 1 PRINT 1 0 0 3 * set value to change near neighbor distance * from 1.902 (=sqrt(1 + phi^2)) to 2.500 RMULTIPLIER 1.31435 RMAX 3.0 POTENTIALS * ipot z label 0 26 Fe 1 26 Fe ATOMS 0.000000 0.000000 0.000000 0 0.000000 1.000000 1.618034 1 0.000000 -1.000000 1.618034 1 0.000000 1.000000 -1.618034 1 0.000000 -1.000000 -1.618034 1 1.000000 1.618034 0.000000 1 -1.000000 1.618034 0.000000 1 1.000000 -1.618034 0.000000 1 -1.000000 -1.618034 0.000000 1 1.618034 1.000000 0.000000 1 1.618034 -1.000000 0.000000 1 -1.618034 1.000000 0.000000 1 -1.618034 -1.000000 0.000000 1
participants (3)
-
Frenkel, Anatoly
-
Matt Newville
-
Scott Calvin