  It has been recently conjectured that bridge functions remain nearly
invariant along phase diagram lines of constant excess entropy for the broad
class of R-simple liquids. To test this hypothesis, the bridge functions of
Yukawa systems are computed outside the correlation void with the
Ornstein-Zernike inversion method and structural input from ultra-accurate
molecular dynamics simulations. The effect of statistical, grid, finite-size,
tail and isomorphic errors is quantified. Uncertainty propagation analysis is
complemented with a detailed investigation of the sensitivity of the bridge
function to periodic and aperiodic multiplicative perturbations in the radial
distribution function. In the long and intermediate range, bridge functions are
demonstrated to be approximately isomorph invariant.
