<HTML>
<HEAD>
<TITLE>Problems with potential overlap for small molecule</TITLE>
</HEAD>
<BODY>
<FONT SIZE="4"><FONT FACE="Calibri, Verdana, Helvetica, Arial"><SPAN STYLE='font-size:11pt'>Dear FEFF users,<BR>
<BR>
I have a FEFF8.4-related question concerning the abortion of the potential module in a XANES calculation for a small molecule in solution. It might not be totally appropriate for the ifeffit mailing list, but the FEFFusers mailing list doesn’t seem to be used very often. I have looked through the archive, but I haven’t found an entry with this specific problem. <BR>
<BR>
I am trying to simulate the L3-edge XANES spectrum for the linear I3- molecule (I-I-I) embedded in a solvent with O-H groups. In order to simulate the scattering from the solvent, I have put two O-H groups at either end of the molecule with the H-atoms pointing towards the I-atoms (expected charge transfer from I3- towards solvent). The I-H and I-O distances are taken from ab initio calculations in the literature, the O-H distance is well-known and amounts to 0.96 A. It is also known from the literature on similar ions in solution that the H-atoms contribute to the XANES and should therefore be included.<BR>
<BR>
I pasted my feff.inp file below.<BR>
<BR>
Running this file results in abortion at the beginning of the potential module:<BR>
</SPAN></FONT><SPAN STYLE='font-size:11pt'><FONT FACE="Consolas, Courier New, Courier"><BR>
...<BR>
Number of processors =            12<BR>
 Calculating potentials ...<BR>
    free atom potential and density for atom type    0<BR>
 MOVRLP-1<BR>
 Abnormal termination on processor             8<BR>
[n7:09539] MPI_ABORT invoked on rank 8 in communicator MPI_COMM_WORLD with errorcode 512<BR>
...<BR>
</FONT><FONT FACE="Calibri, Verdana, Helvetica, Arial"><BR>
It seems to occur when it is trying to overlap the potentials (MOVRLP-1). I tried including the FOLP 2 0.8 card as suggested for H-atoms in the manual, but this does not help.<BR>
<BR>
It seems like the O-H bond <1 A is too small: moving the O atoms further away such that the O-H bond is 1.05 A “solves” the problem, although this is not what I want. But: moving the H atom closer to the I atom while keeping the O positions fixed, which increases the O-H bond as well, does not help! Also, to make things even more complicated, excluding the H-atoms fully, i.e. putting only two O-atoms at 3.5 A from the iodines, results in the WARNING during the SCF loop that the interstitial potential level lies higher than the Fermi energy (moving the O atom closer to ca. 2 A removes this message). In addition, certain combinations of H- and O- positions give the output “negative dens” during the SCF loop, for which I checked that FOLP 2 0.8 did not help...<BR>
<BR>
It is obvious that FEFF has problems with generating the MT potentials for a small and “open” system such as the one here. Did somebody encounter similar problems (abortion, Ei>Ef, negative dens) for systems like this? Is there something wrong with my input cards?<BR>
<BR>
I’d appreciate any help.<BR>
<BR>
Best regards,<BR>
Renske van der Veen<BR>
<BR>
P.s. I am running the parallel version of FEFF8.4 on Linux<BR>
<BR>
***<BR>
<BR>
TITLE  I3- XANES spectrum<BR>
<BR>
*  I L3 edge energy = 4557.0 eV<BR>
EDGE L3<BR>
S02 1.0<BR>
<BR>
CONTROL 1 1 1 1 1 1<BR>
SIG2 0.01<BR>
<BR>
PRINT 1 0 0 0 0 0<BR>
EXCHANGE 0 0.0 0.0 0 <BR>
 <BR>
SCF 10 1 30 0.2 1 <BR>
<BR>
LDOS -30 60 -1 <BR>
<BR>
XANES 6<BR>
<BR>
FMS 10 1     <BR>
UNFR<BR>
<BR>
POTENTIALS<BR>
0 53 I 3 4<BR>
1 53 I 3 4<BR>
2 1 H<BR>
3 8 O 2 3<BR>
<BR>
*FOLP 2 0.8<BR>
<BR>
ATOMS<BR>
0.0000    0.0000    0.0000    0<BR>
-2.91    0.0000    0.0000    1<BR>
2.91    0.0000    0.0000    1<BR>
-5.45    0.0000  0.0000  2<BR>
5.45    0.0000    0.0000    2<BR>
-6.41    0.0000  0.0000  3<BR>
6.41    0.0000    0.0000    3<BR>
END<BR>
</FONT></SPAN></FONT>
</BODY>
</HTML>