your Ni atoms are too close to get reasonnable results.
The pseudopotential spheres (defined by the maximum core radius) should not overlap.
Your pseudopotential radius being 2.050 a.u. (1.1 Ang.), any calculation with atoms closer than 2.2 Ang gives dubious results.
A distance of 1 Ang is much too small. Siesta probably complained about the overlap in the output file.
The same thing applies for abinit.
Jean-Paul Crocombette
At 20:52 11/06/04 +0200, you wrote:
Dear all,
I have made a simple calculation of the forces between 2 Ni atoms at 1 Ang
distance
with siesta I obtain: -89, +89 eV/Ang
with abinit :-321,+321 eV/Ang
the pseudo potential that I used is obtained from atom with the following
input file:
pg Nickel
tm2 0.00
n=Ni c=car
0.0 0.0 0.0 0.0 0.0 0.0
5 4
4 0 2.000 0.000
4 1 0.000 0.000
3 2 8.000 0.000
4 3 0.000 0.000
2.050 2.050 2.050 2.050
the fdf file is:
SystemName Ni
SystemLabel Ni
NumberOfAtoms 2
NumberOfSpecies 1
%block ChemicalSpeciesLabel
1 28 Ni # Species index, atomic number, species label
%endblock ChemicalSpeciesLabel
AtomicCoordinatesFormat Ang
%block AtomicCoordinatesAndAtomicSpecies
0.000 0.000 0.000 1
1.000 0.000 0.000 1
%endblock AtomicCoordinatesAndAtomicSpecies
anyone knows the reason?
best regards,
Raffaele
_________________________________________________________________
Personalizza MSN Messenger con sfondi e fotografie!
http://www.ilovemessenger.msn.it/
SPAM: ------------------------ Spamnix Spam Report ------------------------- SPAM: Spamnix identified this message as spam. This report shows which SPAM: rules matched the message and how many points each rule contributed. SPAM: SPAM: Content analysis details: (102.9 hits, 5.0 required) SPAM: 0.0 BAYES_50 BODY: Bayesian spam probability is 50 to 56% SPAM: [score: 50.4%; H*c:flowed:0 H*c:format:0] SPAM: [anyone:99 knows:99 MSN:2 abinit:2] SPAM: [HTo:D*LISTSERV.UAM.ES:98 SIESTA-L:98 input:98] SPAM: [distance:98 HTo:D*ES:98 HTo:U*SIESTA-L:98] SPAM: [HTo:D*UAM.ES:98 atomic:3 HPrecedence:list:3] SPAM: [H*F:D*COM:96 H*R:U*SIESTA-L:96 H*R:D*es:96] SPAM: [H*R:D*uam.es:96 calculation:5 atoms:5] SPAM: [N:H*r:NN.NNN.NN:90 obtain:90 comparison:12] SPAM: [car:12 H*M:com:86 reason:85 H*r:pickup:15] SPAM: [H*r:service:15 H*r:hotmail.com:15] SPAM: 100 FROM_REJECTED Sender is on Reject list SPAM: [EMAIL PROTECTED] matches [EMAIL PROTECTED] SPAM: 2.9 FAKED_HOTMAIL_DAV X-Originating-Email header does not match From SPAM: SPAM: Spam level: ************************************************** SPAM: --------------------- End of Spamnix Spam Report ---------------------