how to add charge in the system during the calculation

[xingbin]

I am study the acidity for the zeolite, and some literatures define the deprotonation energy (DPE) as DPE = EZ− + EH+ − EZH (1) where EZ−, EH+, and EZH are the respective energies of the zeolite anion after deprotonation, of the gaseous proton, and of the neutral H-form of the zeolite. But I don't know how to calculate the energy of bare Proton. And is there some Parameters in INCAR could perform calculation of zeolite with Charge which is may at different sites based on the different T site and not only for the whole System?
Thanks for your help.

[admin]

Note that works on deprotonation energies, which require the calculation of charged particles,
use molecular codes. For a zeolite such works use a cluster representing the zeolite structure.
In periodical codes, like vasp, the calculation of charged particles is possible, but corrections
for charged systems must be considered, cf. e.g.
(1) Phys. Rev. B 89,045116 (2014)
(2) Phys. Rev. B 91,024107 (2015)

[xingbin]

Thanks for your Response.
Is there a Parameter in INCAR could add the correction during calculation?

[admin]

In distributed version of vasp these corrections are not implemented.
Contact the authors of those papers if their code can be used by
other research groups.