Bond-orbital model for second-order susceptibilities