Spherical Harmonics expansion In the contex of $L^2$ space, it is usually stated that any square-integrable function can be expanded as a linear combination of Spherical Harmonics:
$$
f(\theta,\varphi)=\sum_{\ell=0}^\infty \sum_{m=-\ell}^\ell f_\ell^m \, Y_\ell^m(\theta,\varphi)\tag 2
$$
where $Y_\ell^m( \theta , \varphi )$  are the Laplace spherical harmonics.
The context here is important because this equality holds only in the sense of the $L^2$-norm.
This expansion holds in the sense of mean-square convergence — convergence in [[Lp space|L2]] of the sphere — which is to say that
$$\lim_{N\to\infty} \int_0^{2\pi}\int_0^\pi \left|f(\theta,\varphi)-\sum_{\ell=0}^N \sum_{m=- \ell}^\ell f_\ell^m Y_\ell^m(\theta,\varphi)\right|^2\sin\theta\, d\theta \,d\varphi = 0.$$
So in general this limit is NOT pointwise? So I can't say that the value at a point of the function equals the value of the expansion at the same point?
If so, why it's usually stated out of the context of the structure of Hilbert space, that a bounded function or a square integrable function on the unit sphere can be expanded with Spherical Harmonics if it's not pointwise? I mean in some context outside the Hilbert space, where I am not interested in their square integral.
Furthermore, if it's not pointwise, but only in the norm, am I allowed to sum term by term two different functions, with two different expansions, like in the quantum scattering problem?
 A: As for the reason this expansion is usually not stated in the Hilbert space context, I suspect it is largely cultural. Spherical harmonics are most commonly encountered by physicists and engineers in order to solve specific problems in three-dimensional space, and introductory physics and engineering textbooks typically do not devote a whole lot of time to develop Hilbert space theory. Mathematics textbooks, on the other hand, usually develop Hilbert space theory in full generality, not just for $L^2(\mathbb{R}^3)$  or $L^2(S^2)$. After all, once you have the full machinery of Hilbert spaces, there isn't much to say about spherical harmonics: you just say that they form an orthonormal basis of $L^2(\mathbb{R}^3)$ or $L^2(S^2)$ and you know everything about them from the Hilbert space point of view.
Yes, in general this $L^2$ limit is not pointwise, just like how Fourier series expansions on $L^2([0,1])$ are in general not pointwise limits. However, for the vast majority of functions that you run across in practice, you should be able to show pointwise or even uniform convergence, because most of the time you are going to take spherical harmonic expansions of fairly nice functions. If the coefficients $f_\ell^m$ decay sufficiently rapidly, then you will have uniform convergence, and often you will have ways of knowing the decay of the coefficients (for one, it's tied to the regularity of the functions you're expanding: smoother functions have more rapidly decay coefficients). That's not to say you'll never run into this issue in a physics or engineering context, but often when you do there are still ways around it (whether the physics or engineering literature handles those ways rigorously is a different story).
Lastly, one does not need to be so careful with term by term summation, so long as you're ok with interpreting the results in the Hilbert space context. This is just a general fact about Hilbert spaces: if $f = \sum_n a_n\phi_n$ and $g = \sum_n b_n\phi_n$ are two expansions in a Hilbert space $H$ in terms of an orthonormal basis $(\phi_n)$, then $\sum_n (a_n+b_n)\phi_n$ is a well-defined object in $H$ and can be shown to agree with $f+g$. However, the moment you try to think about pointwise convergence, you return to the previous issue.
