Dear Reza,
There is a bug in siesta, it can't calculate a system with only 1 orbital ...
Regards,
Eduardo


On 04/06/2008, at 2:28, reza behnam wrote:

Dear Edan,
Thank you for your advises. According to Eduardo`s and your advises I set the PAO.Basis as DZP and it worked. But why?just because in SZ for H-atom the splitting for 1s orbital has no meaning or other reasons?

Regards,
Reza

Eduardo
----
M. R. Benam
Visiting Professor,
Department of Physics and Astronomy,
University of British Colombia,
Vancouver, Canada.


----- Original Message ----
From: Edan Scriven <[EMAIL PROTECTED]>
To: SIESTA-L@listserv.uam.es
Sent: Tuesday, May 27, 2008 9:21:30 PM
Subject: Re: [SIESTA-L] WARNING: Qtot

Some thoughts, hope they help!

I've crashed SIESTA before while trying to run calculations on one atom.
I don't remember if my error was the same as yours, but my problem was
that I was trying to run it on multiple cpus. The problem went away with
the same input on one cpu.

Take out everything from your .fdf that you don't actually need,
which'll help diagnose the problem. For example, you don't need the band
structure section for an isolated atom. Take out the lattice constant
and let SIESTA calculate its own (it'll pick something that guarantees
your basis functions fall off to zero before the cell boundary). Come to think of it, you can comment out every line from "# MD options" on. Once you've minimised your .fdf, you can mess around with one line at a time
to find the one that's causing trouble.

Try other systems a step up in complexity, e.g. water or methane, and
see if your problem persists. This is also a way to test the sanity of
your pseudopotentials, but won't help you with the one atom, multiple
cpus thing.

Good luck!
Edan.






Reply via email to