[ase-users] Can not use NEB with SIESTA?

Jens Jørgen Mortensen jensj at fysik.dtu.dk
Wed Feb 18 09:44:04 CET 2009


On Mon, 2009-02-16 at 04:12 -0800, Nguyen Ngoc Ha wrote:
> Dear Developers and ASE users,
> As I told in my mail sent at our forum
> (https://listserv.fysik.dtu.dk/pipermail/ase-users/2009-January/000177.html), I've never gotten the convergence of NEB with SIESTA. 
> I think that the season is in scf solution of the oder-n in SIESTA :
> To calculate Hamiltonian must use new  charge density from
> previous wavefunctions so calculating Hamiltonian of (i+1)-th image
> will use i-th wavefunctions. To avoid this I think wavefunctions of
> each image should be contained in a separate folder. 
> Sorry if I'm wrong,

It looks like you are using the same Siesta calculator instance for all
images.  Try to use one instance for each image and use different labels
for each of them:

for i in range(3):
    ...
    calc = Siesta('image%d' % i, ...)
    image.set_calculator(calc)

Jens Jørgen





More information about the ase-users mailing list