4.1. Crystal data

The crystal cell data – including lattice constants and angles, the space group symbol, and the elements of the shift vector – are placed in text boxes for easy editing. The coordinates of the unique sites are listed in the grid at the bottom of the window. The absorber is chosen by clicking one of the boxes on the Core column.

Remember that FEFF considers numbers with 5 digits of precision after the decimal point. 0.333 is not the same thing as 0.33333. You may, however, enter things like 1/3 and avoid the precision issue entirely.

As a new feature compared to earlier versions of ATOMS, there are two radial distances. The cluster size determines the extent of the cluster expanded into the feff.inp. This should usually be set to something rather large, 9 Å is often a good default. This probably (but not always!) assures that the cluster in the feff.inp file is adequately large to include all unique potentials and has all atom types sufficiently well bounded that the muffin tin potentials are likely to be be computed reasonably well.

The second distance will set the value of RMAX in the feff.inp file. In general, you do not want this to be much larger than the extent of the data you intend to analyze. 5 Å or 6 Å is usually the largest sensible value for longest path. The reason for this is that the pathfinder part of FEFF has been rewritten for this version of ARTEMIS. While the new pathfinder implementation offers a number of useful new features, it is substantially slower than FEFF's native pathfinder. In any case, there is no benefit to computing paths that you will never use in your fit.

The absorption edge for the calculation is chosen from the menu to the left of the lattice constant area. This is usually determined from the input data, but may need to be explicitly selected. If not specified in the atoms.inp file, the edge will be set to K for element lighter than Ce (Z=58), and to LIII for heavier elements.

The style menu is another new feature in this version of ATOMS. It is used to set how the list of unique potentials is determined from the elements in the atoms list. The choices are

elements
Each unique element species is assigned a potential number.
tags
Each unique tag is assigned a potential number.
sites
Each crystallographic site is assigned a potential number.

Remember that FEFF only allows for 7 unique potentials other than the absorber. The tags and sites options can often result in more than 7 potnatials, which will result in an unrunnable feff.inp file. Specifying unique potentials by tags is a good way of differentiating between dissimilar atoms of the same species. For example, in an oxygenyl species, it is often useful to give the axial oxygen atoms a different potential from the remaining oxygens by using the tags option.

Here is an example of an atoms.inp file for sodium uranyl acetate, which contains two very short axial oxygen atoms double bonded to the uranium atoms at about 1.8 Å and a number of equatorial oxygen atoms at a much longer distance. The axial and equatorial oxygen positions are distinguished by their tags and will given separate unique potentials when using the tags style.

title        Templeton et al.
title        Redetermination and Absolute configuration of Sodium Uranyl(VI) triacetate.
title        Acta Cryst 1985 C41 1439-1441
space = P 21 3
a =   10.6890        b =      10.6890        c =      10.6890
alpha =       90.0   beta =   90.0   gamma =  90.0
core =       U       edge =  L3
atoms
! elem   x          y          z       tag        occ
  U     0.42940    0.42940    0.42940  U             1.00000
  Na    0.82860    0.82860    0.82860  Na            1.00000
  O     0.33430    0.33430    0.33430  Oax           1.00000
  O     0.52420    0.52420    0.52420  Oax           1.00000
  O     0.38340    0.29450    0.61100  Oeq           1.00000
  O     0.54640    0.24430    0.50070  Oeq           1.00000
  C     0.47860    0.22600    0.59500  C             1.00000
  C     0.50880    0.12400    0.68620  C             1.00000

The assignment of potential indeces is explained in detail and with examples in a later chapter.

4.1.1. Polarization

A FEFF calculation considering linear polarization can be triggered by setting one or more non-zero values for the polarization vector.

This vector sets the value of the POLARIZATION keyword in the resulting feff.inp. The value written in the feff.inp file is the value that will be used in the pathfinder and when computing the path contributions. That is, if you edit the POLARIZATION in the feff.inp file, the edited value will take precedence over the value specified here.

Caution

FEFF's ELLIPTICITY keyword is not supported at this time. That means the trick of modeling “polarization in the plane” is not yet supported by ARTEMIS.




DEMETER is copyright © 2009-2016 Bruce Ravel – This document is copyright © 2016 Bruce Ravel

This document is licensed under The Creative Commons Attribution-ShareAlike License.

If DEMETER and this document are useful to you, please consider supporting The Creative Commons.