def phs_correction(z_h, z_e, de_vbm, de_cbm):
    """ Compute the PHS correction
    :param z_h: number of holes in the PHS
    :param z_e: number of electrons in the PHS
    :param de_vbm: correction of the VBM
    :param de_cbm: correction of the CBM"""

    return - z_h * de_vbm, z_e * de_cbm