The variance inflation factor (VIF) in an ordinary least squares linear regression coefficient is calculated using the $R^2$ of a linear model that uses the other features to predict the feature to which the coefficient corresponds.
If the feature for which we want to calculate the variance inflation factor is binary, then such a model is a linear probability model, which has several known problems, among them being that such a model can make predictions that are literally impossible (e.g., predicted probability above $1$).
I suppose the math of the working with the feature covariances leads to such a model, but it seems bizarre that such a characteristic (VIF) would relate to a model that can make impossible predictions. Is there any resolution beyond a somewhat unsatisfying, "Yep, that's just what happens when you work with the feature covariance matrix"?