 This thread is locked.Only browsing is available.
 Top Page > Browsing uniform charged background Date: 2005/11/30 23:13 Name: peter Dear Prof. Ozaki,I am trying to utilize openmx for small organic molecules calculations. I have a few questions regarding molecular ions.Since recently you have scf.system.charge variable in your .dat file. The comment in the manual suggests that there is a substraction of the uniform compensating background in your scf calculation. I could not find a line in the source code where such substraction is done. Do you really compensate for the extra charge (doping)? The only reference to sfc.system.charge I could find is in the SetPara_DFT.c. It's used there to defind the Fermi level (to the best of my understanding, of course).The thing is that i would like to compare the energy of two different charged states of the same molecule and therefore want to know exactly what is already included in the Uele.Thank you for the great product and your attentionPeter Page: Re: uniform charged background ( No.1 ) Date: 2005/12/01 13:36 Name: T.Ozaki Hi,The charge compensation is implicitly made. The trick relies on the fact that the FFT is used for solving Poisson's equation in OpenMX. In OpenMX the charge density is decomponsed into twocontributions (see PRB 72, 045121): superposition of atomic charge and the difference chage, and the sum of two is exactly the same as that calculated from KS wave functions. Then, the difference charge is Fourier-transformed to find the Hartree potential arising from the difference charge. If the system is neutral, the Fourier component with q=0 is zero, since the Fourier component with q=0 corresponds tothe integration of the difference charge over the unit cell.On the other hand, if the system is charged up negativelyor positively, the Fourier component with q=0 is not zero.This brings a problem in finding the Hartree potential becauseof the division by zero (1/|q|^2, where q=0). Thus, the Fourier component with q=0 is set in zero, which meansthat a uniform charge density with the opposite charge is assumed.This is the meaning for "The charge compensation is implicitly made"(see also Rev.Mod.Phys. 64, 1045).However, it should be noted that in this 'implicit compensation' an artificial Coulomb interaction between charged atoms in the periodically copied cells could be a source of errors as discussed in PRB 51, 4014., while the remedy is not implemented in OpenMX.In a practical sense, for molecular systems you will need to check how the energy diffience between different charged systems behaves as a function of the cell size.Best regards,T.Ozaki Page: