[Pw_forum] ion problem
Giovanni Cantele
Giovanni.Cantele at na.infn.it
Wed Nov 14 09:20:17 CET 2007
Cyrille Barreteau wrote:
> Dear all
>
> I have a question related to this problem.
> I am interested in the study of charged clusters.
> If I understand well, due to to the periodicity of the plane wave codes
> one has to substract the correct "background" in order to keep the
> energy from going to infinity..
>
This is automatically done by the code, as you may read in the INPUT_PW
documentation file:
nelec REAL
number of electron in the unit cell
(may be noninteger if you wish)
Default: the same as ionic charge (neutral cell)
A compensating jellium background is inserted
to remove divergences if the cell is not neutral
So, you set nelec in the input file, and that should be enough.
> I would like to compare the energy of clusters with different
> charges -e -2e etc... Do you have some advices to do such calculation?
> Is it necessary to take even larger supercells than for neutral clusters?
>
The problem with charged clusters is that spurious interactions are
introduced by the
periodic boundary conditions, so as convergence would be achieved at
very large supercells.
Neverthelessm the total energy can be corrected, to get rid of such
interactions, as
described by Makov-Payne, PRB (1995). Such corrections are currently
implemented
in Quantum-ESPRESSO, with the only warning that, as far as I know,
corrections are
applied to the total energy but are not to the forces.
If you are using espresso-3.2, you should
- run a self-consistent calculation for you charged cluster
- run post-processing code pp.x, for obtaining the total charge density
(plot_num=0) in
a suitable format
- run dipole.x, you may read the header of PP/dipole.f90 for very short
explanation
You finally end up with a Makov-Payne correction that must be added to
the total energy
WITH its sign.
You may try the same charged cluster with two different supercells, and
compare the
two uncorrected energies, as well as the corrected ones, finding the
latter much more
converged.
Next releases will probably have an easier procedure. In the CVS version
I understand that a new variable
(assume_isolated) has been added to the pw.x input, that, if set to
.TRUE., should allow to get
the MP correction directly within the scf loop (so post-processing is
not mandatory).
Giovanni
--
Dr. Giovanni Cantele
Coherentia CNR-INFM and Dipartimento di Scienze Fisiche
Universita' di Napoli "Federico II"
Complesso Universitario di Monte S. Angelo - Ed. 6
Via Cintia, I-80126, Napoli, Italy
Phone: +39 081 676910
Fax: +39 081 676346
E-mail: Giovanni.Cantele at na.infn.it
Web: http://people.na.infn.it/~cantele
More information about the Pw_forum
mailing list