Re: [SIESTA-L] Consulta sobre NetCharge

2009-06-11 Thread David Strubbe
Mariana,

Cada átomo de sodio tiene un electrón, entonces NetCharge 0.8 remueve 0.8
electrones, y quedan 0.2 electrones.  SIESTA añade una densidad uniforme a
la supercelda para compensar la carga neta del sistema, para hacerlo
neutral.  Por eso, hay que converger el tamaño de la supercelda en la
dirección perpendicular a la lámina, para asegurar que esta densidad
uniforme sea bastante pequeña y no tenga una influencia espuria en los
resultados.

(Aunque SIESTA es un proyecto español, el idioma estándar de esta lista de
correo es inglés.)

David

2009/6/11 Mariana Isabel Rojas 

> Hola a todos:
>
> Yo estoy trabajando con una lámina de Na que la describo con un único átomo
> y condiciones periódicas, con lo cual es como si tuviera un plano infinito.
> A este sistema le asigno una carga NetCharge 0.8 para emular iones de Na+.
> Parece que siesta realiza los cálculos con 0.2 electrones, porque el
> momento
> total de espín es 0.2 y cuando integro la densidad electrónica me da
> también
> 0.2. En el perfil de densidad electrónica en la dirección perpendicular a
> la
> lámina veo un máximo de densidad electrónica en la posición del átomo.
>
> Me gustaría saber:
> ¿Cómo SIESTA compensa la carga que uno le asigna al sistema con el comando
> NetCharge?
>
> Dra. Mariana Isabel Rojas
> Profesora Adjunta
> Investigadora Adjunta de CONICET
> Unidad de Matemática y Física
> Facultad de Ciencias Químicas
> Universidad Nacional de Córdoba
> Ciudad Universitaria, 5000 Córdoba
> Tel. 54-351-4344972
>


[SIESTA-L] Consulta sobre NetCharge

2009-06-11 Thread Mariana Isabel Rojas
Hola a todos:

Yo estoy trabajando con una lámina de Na que la describo con un único átomo
y condiciones periódicas, con lo cual es como si tuviera un plano infinito.
A este sistema le asigno una carga NetCharge 0.8 para emular iones de Na+.
Parece que siesta realiza los cálculos con 0.2 electrones, porque el momento
total de espín es 0.2 y cuando integro la densidad electrónica me da también
0.2. En el perfil de densidad electrónica en la dirección perpendicular a la
lámina veo un máximo de densidad electrónica en la posición del átomo.

Me gustaría saber:
¿Cómo SIESTA compensa la carga que uno le asigna al sistema con el comando
NetCharge?

Dra. Mariana Isabel Rojas
Profesora Adjunta
Investigadora Adjunta de CONICET
Unidad de Matemática y Física
Facultad de Ciencias Químicas
Universidad Nacional de Córdoba
Ciudad Universitaria, 5000 Córdoba
Tel. 54-351-4344972