Let $f(n)$ be a nonnegative arithmetic function satisfying
- $f(p^l) \leq A_1^l$ for all primes $p$, integers $l\geq 1$, and some constant $A_1$;
- $f(n) \leq A_2 n^\varepsilon$ for all $\varepsilon > 0$ and some constant $A_2 = A_2(\varepsilon)$.
From the work of Shiu, we know that for functions $f$ satisfying 1. and 2., we have the estimate $$ \tag{1} \sum_{x-y<n\leq x} f(n) \ll \frac{y}{\log x} \exp\left(\sum_{p\leq x} \frac{f(p)}{p}\right) $$ uniformly for $x^\alpha < y < x$ and any fixed $\alpha \in (0,1)$.
I'd like to know if such a bound has been proved for functions like $$ f(n) = \sum_{d^2\mid n} d^\theta, $$ where $0 <\theta < 1$. This function fails both conditions 1. and 2. Nevertheless, it is bounded on average. One can use Dirichlet's hyperbola method to establish (for some constant $C = C(\theta)$) $$ \sum_{n\leq x} \sum_{d^2\mid n} d^\theta = \zeta(2-\theta)x + Cx^{\frac{1+\theta}{2}} + O\left(x^{\frac{1+\theta}{3}}\right), $$ which immediately gives a result like (1) in the range $y \geq x^{\frac{1+\theta}{3}}$. Certainly a result like (1) cannot hold in a range like Shiu's work; indeed, for any positive integer $m$, take $x = m^2$ to see that $$ \sum_{x - y < n\leq x} \sum_{d^2\mid n} d^\theta \geq m^\theta = x^{\frac{\theta}{2}}, $$ so (1) fails for $y \ll x^{\frac{\theta}{2}}$. I conjecture that this is the only obstruction, so a result like (1) should hold for $y > x^{\frac{\theta}{2}+\alpha}$ and any $\alpha \in (0,1)$.
This certainly seems like the thing someone might have considered before, but I'm not familiar enough with the literature to know. Any references, ideas, or feedback would be most appreciated.