def over_the_barrier_intensity(i_p, z):
    """ i_p is the ionization potental and z in the charge of the ion (after ionization).
    i_p in ev

    return unit (W/cm^2)
    """
    return 4*10**9 * i_p**4 / z**2