How do you derive the following distribution \[ ln S_u - ln S_t \sim N\big( \mu (u-t), \sigma^2(u-t)\big) \]
It is obtained by solving the SDE for geometric Brownian motion: dSt = St ( (mu + 0.5*sigma^2) dt + sigma dBt) where Bt is standard Brownian motion. Note that this SDE (and hence distribution) is parameterised differently to the SDE for geometric Brownian motion used in conjunction with the Black-Scholes model: dSt = St (mu dt + sigma dBt)
The two SDEs relate to the same basic model (geometric Brownian motion) just with the parameters defined differently. They both give rise to a lognormal distribution, but with different parameters, so I can't quite say they're "equivalent". The SDE: dSt = St ( (mu + 0.5*sigma^2) dt + sigma dBt) appears rarely in CT8 - I quoted it as it answers the question you asked - how that particular form of the lognormal distribution is derived.