We eventually settled for the model discussed here because it produces better translations during decoding.
since formula (4) overestimates the numerator and denominator equally, the approximation works well in practice.