Instantaneous solvent term in a relaxation modulus ($G={\mu }_0{\delta }_{+}+G_0{e}^{-s/\tau }$)

For a linear viscoelastic fluid, write the solvent stress explicitly as $2{\mu }_{0}E(t)$, or use the displayed kernel with a causal endpoint delta satisfying ${\int }_0^{\infty }{\delta }_{+}(s)f(s)ds=f(0)$. If a symmetric delta is integrated with half its weight at the endpoint, its coefficient must instead be $2{\mu }_0$. Mixing those conventions loses a factor of two. The Oldroyd-B model has the displayed exponentially decaying polymer part after linearization.