[SIESTA-L] Energy does not converge!! very urgent..please help someone...

Tue, 13 Jul 2010 15:46:08 -0700 

I have made a slab with FCC configuration and I want to calculate the energy
of the slab!! The slab contains 5 layers of Nickel atoms and a carbon atom
on the top of it. I just want to calculate the energy of this system!!!
....I made a supercell having 5X5X5 lattice constants in each direction ( 5a
5a 5a 60 60 60) and then gave a BIG vaccum space so that the nearby
Replicated supercells do not interact.(Is this a correct concept??????) .
The problem is that the energies that i get are very large!!!! about
-50KeV!! Can any body tell what is wrong in my fdf input file....it's very
urgent!!

---------------------------------------------------FDF
FILE----------------------------------------------------------------------------------------------------------------------------------------------

%block ChemicalSpeciesLabel
 1  28  Ni
 2  6   C
%endblock ChemicalSpeciesLabel

SaveRho .true.

LatticeConstant 1.0 Ang

%block LatticeParameters
99.6  99.6  99.6  60.  60.  60.
%endblock LatticeParameters

%block AtomicCoordinatesOrigin
0.000  0.000  0.000
%endblock AtomicCoordinatesOrigin



MD.TypeOfRun CG
MD.NumCGSteps 0

AtomicCoordinatesFormat  Fractional

%block AtomicCoordinatesAndAtomicSpecies
0.000 0.000 0.000 1
0.025 0.000 0.000 1
0.050 0.000 0.000 1
0.075 0.000 0.000 1
0.000 0.025 0.000 1
0.025 0.025 0.000 1
0.050 0.025 0.000 1
0.075 0.025 0.000 1
0.000 0.050 0.000 1
0.025 0.050 0.000 1
0.050 0.050 0.000 1
0.075 0.050 0.000 1
0.000 0.075 0.000 1
0.025 0.075 0.000 1
0.050 0.075 0.000 1
0.075 0.075 0.000 1
0.000 0.000 0.025 1
0.025 0.000 0.025 1
0.050 0.000 0.025 1
0.075 0.000 0.025 1
0.000 0.025 0.025 1
0.025 0.025 0.025 1
0.050 0.025 0.025 1
0.075 0.025 0.025 1
0.000 0.050 0.025 1
0.025 0.050 0.025 1
0.050 0.050 0.025 1
0.075 0.050 0.025 1
0.000 0.075 0.025 1
0.025 0.075 0.025 1
0.050 0.075 0.025 1
0.075 0.075 0.025 1
0.000 0.000 0.050 1
0.025 0.000 0.050 1
0.050 0.000 0.050 1
0.075 0.000 0.050 1
0.000 0.025 0.050 1
0.025 0.025 0.050 1
0.050 0.025 0.050 1
0.075 0.025 0.050 1
0.000 0.050 0.050 1
0.025 0.050 0.050 1
0.050 0.050 0.050 1
0.075 0.050 0.050 1
0.000 0.075 0.050 1
0.025 0.075 0.050 1
0.050 0.075 0.050 1
0.075 0.075 0.050 1
0.000 0.000 0.075 1
0.025 0.000 0.075 1
0.050 0.000 0.075 1
0.075 0.000 0.075 1
0.000 0.025 0.075 1
0.025 0.025 0.075 1
0.050 0.025 0.075 1
0.075 0.025 0.075 1
0.000 0.050 0.075 1
0.025 0.050 0.075 1
0.050 0.050 0.075 1
0.075 0.050 0.075 1
0.000 0.075 0.075 1
0.025 0.075 0.075 1
0.050 0.075 0.075 1
0.075 0.075 0.075 1
0.000 0.000 0.100 1
0.025 0.000 0.100 1
0.050 0.000 0.100 1
0.075 0.000 0.100 1
0.000 0.025 0.100 1
0.025 0.025 0.100 1
0.050 0.025 0.100 1
0.075 0.025 0.100 1
0.000 0.050 0.100 1
0.025 0.050 0.100 1
0.050 0.050 0.100 1
0.075 0.050 0.100 1
0.000 0.075 0.100 1
0.025 0.075 0.100 1
0.050 0.075 0.100 1
0.075 0.075 0.100 1
0.033 0.033 0.125 2
0.100 0.000 0.000 1
0.100 0.025 0.000 1
0.100 0.050 0.000 1
0.100 0.075 0.000 1
0.100 0.100 0.000 1
0.100 0.000 0.025 1
0.100 0.025 0.025 1
0.100 0.050 0.025 1
0.100 0.075 0.025 1
0.100 0.100 0.025 1
0.100 0.000 0.050 1
0.100 0.025 0.050 1
0.100 0.050 0.050 1
0.100 0.075 0.050 1
0.100 0.100 0.050 1
0.100 0.000 0.075 1
0.100 0.025 0.075 1
0.100 0.050 0.075 1
0.100 0.075 0.075 1
0.100 0.100 0.075 1
0.100 0.000 0.100 1
0.100 0.025 0.100 1
0.100 0.050 0.100 1
0.100 0.075 0.100 1
0.100 0.100 0.100 1
0.000 0.100 0.000 1
0.025 0.100 0.000 1
0.050 0.100 0.000 1
0.075 0.100 0.000 1
0.000 0.100 0.025 1
0.025 0.100 0.025 1
0.050 0.100 0.025 1
0.075 0.100 0.025 1
0.000 0.100 0.050 1
0.025 0.100 0.050 1
0.050 0.100 0.050 1
0.075 0.100 0.050 1
0.000 0.100 0.075 1
0.025 0.100 0.075 1
0.050 0.100 0.075 1
0.075 0.100 0.075 1
0.000 0.100 0.100 1
0.025 0.100 0.100 1
0.050 0.100 0.100 1
0.075 0.100 0.100 1
%endblock AtomicCoordinatesAndAtomicSpecies

SolutionMethod diagon





Thnx in advance!!

Responder a