def makov_payne_correction(defect_cell, geometry, e_r, mk_1_1):
    """ Makov-Payne correction
    :param defect_cell: Cell object of the defect cell calculation
    :param geometry: geometry of the host cell
    :param e_r: relative permittivity
    :param mk_1_1: Value of the first term of the Makov-Payne correction in the case q = 1 & e_r = 1 """

    c_sh_dico = {'sc': -0.369, 'fcc': -0.343, 'bcc': -0.342, 'hcp': -0.478, 'other': -1./3}
    c_sh = c_sh_dico[geometry]

    return (1. + c_sh * (1. - 1./e_r)) * defect_cell.charge ** 2 * mk_1_1 / e_r