c2 = sp.exp(sp.I * sp.log(5))
2_standard = sp.cos(sp.log(5)) + sp.I * sp.sin(sp.log(5))

You have to tell sympy explicitly what to evaluate. expand(complex=true) doesn't automatically recognize non exponential forms as complex exponentials