def birch_murnaghan(volumes, volume0, bulk_modulus0, bulk_deriv0):
    """
    This evaluates the Birch Murnaghan equation of states
    """
    PV = []
    EV = []
    v0 = volume0
    bm = bulk_modulus0
    dbm = bulk_deriv0
    for vol in volumes:
        pv_val = 3 * bm / 2.0 * ((v0 / vol) ** (7 / 3.0) - (v0 / vol) ** (5 /
            3.0)) * (1 + 3 / 4.0 * (dbm - 4) * ((v0 / vol) ** (2 / 3.0) - 1))
        PV.append(pv_val)
        ev_val = 9 * bm * v0 / 16.0 * ((dbm * (v0 / vol) ** (2 / 3.0) - 1) **
            3 * ((v0 / vol) ** (2 / 3.0) - 1) ** 2 * (6 - 4 * (v0 / vol) **
            (2 / 3.0)))
        EV.append(ev_val)
    return EV, PV
