I implemented the approximation and the calibration seems to work fine now. However, if I use constant values for sigma_0 and beta, I get a non-stationary forward skew. In other words, the n months forward skew with, say, maturity 4 months seems to be a (linearly) decreasing function of n (see chart...