diff --git a/sumpy/expansion/multipole.py b/sumpy/expansion/multipole.py
index 92ae80f3308b4bf8820fbc00da9229725058e3e1..47e33a8c76ae8ad5f85e72efcab7415fbd648253 100644
--- a/sumpy/expansion/multipole.py
+++ b/sumpy/expansion/multipole.py
@@ -88,7 +88,7 @@ class VolumeTaylorMultipoleExpansionBase(MultipoleExpansionBase):
             for i, mi in enumerate(coeff_identifiers):
                 result[i] /= (mi_factorial(mi) * rscale ** sum(mi))
         else:
-            avec = avec * rscale**-1
+            avec = [sym.UnevaluatedExpr(a * rscale**-1) for a in avec]
 
             result = [
                     mi_power(avec, mi) / mi_factorial(mi)