A simple bridge function for strongly coupled Yukawa systems is developed b
ased upon a form previously extracted for a one-component plasma [H. Iyetom
i, S. Ogata, and S. Ichimaru. Phys. Rev. A 46, 1051 (1992)]. Using the prop
osed bridge function in the modified hypernetted chain theory, excellent ag
reement is obtained with molecular dynamics simulations and the compressibi
lity sum rule is satisfied to within a few percent. This result offers a si
mple and very accurate method to quickly compute correlation functions for
Yukawa systems.