By ‘correct’ we mean one that produces autocorrelations that fit empirical autocorrelations, i.e. Mittag-Leffer functions .
Apply to obtain
Integrate the second quantity by parts using the fact that to get
This then still leaves the fractional integration in the first term of (*). So this is partly the problem with using fractional Brownian motions only in the model it would seem.