Sungyun Kim
Sungyun Kim
In `m_fit_polynomial_coeff.F90`, the parameter `HaBohr_meVAng` is the conversion factor from `Ha/Bohr` to `eV/Ang` not `meV/Ang`.
I got unphysical values (negative values from the optimized structure) of elecstic tensor when using DFPT + PAW + PBESol. When using lda or nc psp, there are no problem....
Is it possible to include the dipole-dipole interaction term described in Zhong, Vanderbilt, Rabe's [1994 PRL](https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.73.1861) and [1995 PRB](https://link.aps.org/doi/10.1103/PhysRevB.52.6301) papers in CASMcode?