Abstract:We study diffusion and butterfly velocity (v B ) in two holographic models, linear axion and axion-dilaton model, with a momentum relaxation parameter (β) at finite density or chemical potential (µ). Axion-dilaton model is particularly interesting since it shows linear-T -resistivity, which may have something to do with the universal bound of diffusion. At finite density, there are two diffusion constants D ± describing the coupled diffusion of charge and energy. By computing D ± exactly, we find that in the incoherent regime (β/T 1, β/µ 1) D + is identified with the charge diffusion constant (D c ) and D − is identified with the energy diffusion constant (D e ). In the coherent regime, at very small density, D ± are 'maximally' mixed in the sense that D + (D − ) is identified with D e (D c ), which is opposite to the case in the incoherent regime. In the incoherent regime D e ∼ C − v 2 B /k B T where C − = 1/2 or 1 so it is universal independently of β and µ. However, D c ∼ C + v 2 B /k B T where C + = 1 or β 2 /16π 2 T 2 so, in general, C + may not saturate to the lower bound in the incoherent regime, which suggests that the characteristic velocity for charge diffusion may not be the butterfly velocity. We find that the finite density does not affect the diffusion property at zero density in the incoherent regime.