There is already a huge amount of literature on this topic. As a start I would recommend [1].
As is usually the case with models, the parameters have been tuned for certain locations, so it makes sense to test them at many locations. In general one can say that Erbs usually overestimates the diffuse component under clear sky conditions, and this bias is severe (of the order of 25% using 10 minute data from a location in Western Europe). You won't see that statistic when looking at all the data (all sky), but if you specifically look at clear sky conditions, you will see that basically Erbs overestimates the atmospheric turbidity / aerosol optical depth, leading to a higher diffuse component. Put another way, the sky is clearer than what Erbs models it to be, there is less scattering due to aerosols and more direct light. This error leads to underestimation of the direct component and underestimation of shading losses in PVsyst
It would definitely make sense to include DIRINT [2] in PVsyst, since that model uses more predictors, not just the clearness index. Specifically, DIRINT has a time component that takes into account variability, i.e. the change in clearness index. That helps then to distinguish between different weather conditions, and thus model them differently.
All of that being said, the best is to import diffuse irradiance, which is however not always possible if it is not measured on site.
[1] Gueymard, C. A., & Ruiz-Arias, J. A. (2016). Extensive worldwide validation and climate sensitivity analysis of direct irradiance predictions from 1-min global irradiance. Solar Energy, 128, 1–30. https://doi.org/https://doi.org/10.1016/j.solener.2015.10.010
[2] Perez, R. R., Ineichen, P., Maxwell, E. L., Seals, R. D., & Zelenka, A. (1992). Dynamic global-to-direct irradiance conversion model. ASHRAE Transactions, 354–369.