LEPSILON vs LOPTICS
Posted: Fri Apr 09, 2010 8:07 pm
I'm trying to determine both the static and frequency-dependent dielectrical matrix of a sample of water molecules using, respectively, a calculation with LEPSILON=.TRUE. and another with LOPTICS=.TRUE. In the latter, I set NEDOS=2000 and use a variety of value of NBANDS, ranging between 2 and 3 times the VASP default, as recommended in the VASP manual. All the other parameters are the same for both the LEPSILON and LOPTICS calculations.
The manual also recommends verifying that the zero-frequency value obtained using the LOPTICS calculation should exactly match the static value obtained from the LEPSILON calculation. Here is the problem. They are NOT the same. the LOPTICS calculation gives around 1.5 (with variations only in the second digit by increasing NBANDS), whereas the LEPSILON calculation gives 1.95, which is significantly higher.
Any idea as to the origin of this discrepancy?
The manual also recommends verifying that the zero-frequency value obtained using the LOPTICS calculation should exactly match the static value obtained from the LEPSILON calculation. Here is the problem. They are NOT the same. the LOPTICS calculation gives around 1.5 (with variations only in the second digit by increasing NBANDS), whereas the LEPSILON calculation gives 1.95, which is significantly higher.
Any idea as to the origin of this discrepancy?