Yes, I got the same result.

When I calculated the derivative of transition dipole moment (\mu_{fi}) respect to the normal coordinate (Q_k), I set the normal displacement 0.1, 0.2, ...... , 1.0. Then x component of the transition dipole moment is like: -4.182, -4.181, 4.180, -4.180, 4.179, ...

Because the sign of transition dipole moment <e| \mu_x |g> is undetermined, I set the result their absolute value: 4.182, 4.181, 4.180, 4.180, 4.179, ... Then I got a smooth curve of \mu_x(Q_k) respect to Q_k.

In such a calculation, I can not find an unique criterion to fix the phase!

Does anyone has some suggestions?