As the term "summatory function" is somewhat ambiguous, the function I mean is $\Phi_2(n):=\sum _{d|n}\phi(d)^2$, where $\phi$ is Euler's totient. It is well-known that $\Phi_1(n):=\sum _{d|n}\phi(d)=n$. All powers of $\phi$ are multiplicative, hence their summatory functions (over divisors) $\Phi_k(n)$ are also multiplicative, and using that one can get a product formula for them after computing their value on prime powers. I am curious if anything is known beyond that, like Ramanujan's averaged large $n$ asymptotics for the squared divisor sum $\sigma_2(n)$. For a different kind of summatory function, $\sum_{k=1}^n \phi(k)^2$, the asymptotics are derived in Average order of Eulers totient function squared.
$\Phi_2(n)$ comes up as the dimension of a certain algebra, but explaining what it is would take too much space. Was $\Phi_2(n)$ studied, does it 'oscillate' around $cn^2$ for some constant $c$? References to literature are welcome.