The eigenfunctions of the Laplacian on $SL(2,\mathbb Z)\backslash \mathbb H$ are known to be given by three types: the constant function, the [real analytic Eisenstein series](https://en.wikipedia.org/wiki/Real_analytic_Eisenstein_series) (which come in a continuous family), and the Maass cusp forms (which come in a discrete family -- the first few given e.g. [here](http://www.lmfdb.org/ModularForm/GL2/Q/Maass/1/0/)). Let us denote $E\_s(\tau, \bar\tau)$ as the real analytic Eisenstein series which scales as $y^s + \frac{\Lambda(1-s)}{\Lambda(s)}y^{1-s}$ as $y \rightarrow \infty$, where $y = \text{Im}(\tau)$, $\Lambda(s) = \pi^{-s} \Gamma(s)\zeta(2s)$; and $\nu\_i$ denote the $i^{\text{th}}$ discrete cusp form. 
My question is: Is there a simple expression for $\int\_{SL(2,\mathbb Z)\backslash \mathbb H} \frac{dxdy}{y^2} E\_{s\_1}(\tau, \bar\tau) E\_{s\_2}(\tau, \bar\tau) \overline{\nu\_i(\tau,\bar\tau)}$? (We can take $s\_i \in 1/2 + i \mathbb R$.) For example, if there were only one $E\_s$ in the previous expression, this would vanish due to orthonormality. What about when there is a product of two of them?


The motivation is the following --- it is known that any modular invariant function in $L^2$ can be expanded in Eisensteins on the principal series and Maass cusp forms. If we take the product of two modular invariants in $L^2$ each without support on the discrete cusp forms, does the product necessarily also have no support on the discrete series (let us assume the product remains in $L^2$, or at least has good decomposition..)? 
Similarly are there simple or known expressions for:


$\int\_{SL(2,\mathbb Z)\backslash \mathbb H} \frac{dxdy}{y^2} E\_{s\_1}(\tau, \bar\tau) E\_{s\_2}(\tau, \bar\tau) \overline{E\_{s\_3}(\tau,\bar\tau)}$


$\int\_{SL(2,\mathbb Z)\backslash \mathbb H} \frac{dxdy}{y^2} E\_{s\_1}(\tau, \bar\tau) \nu\_{i\_1}(\tau, \bar\tau) \overline{\nu\_{i\_2}(\tau,\bar\tau)}$


$\int\_{SL(2,\mathbb Z)\backslash \mathbb H} \frac{dxdy}{y^2} \nu\_{i\_1}(\tau, \bar\tau) \nu\_{i\_2}(\tau, \bar\tau) \overline{\nu\_{i\_3}(\tau,\bar\tau)}$
known in the literature?


(The first of these three may not converge but hopefully can be written with delta functions or something.) Thank you! 



Spectral decomposition of product of modular functions
I'm guessing you're a physicist by your unusual notation! The real-analytic Eisenstein series is usually denoted $E(z,s)$ (where $z = \tau$ in your notation), not $E\_s(\tau,\bar{\tau})$ (why does $\bar{\tau}$ appear?). I'll also let $f(z)$ denote a cusp form rather than $\nu\_i(\tau,\bar{\tau})$ (I've never seen this notation used).




---


Anyway, the question asks (among other things) the following. Given two cusp forms $f\_1,f\_2$ and an Eisenstein series $E(z,s)$, determine a closed formula for
$$\int\_{\Gamma \backslash \mathbb{H}} f\_1(z) f\_2(z) E(z,s) \, d\mu(z).$$
Here $\Gamma = \mathrm{SL}\_2(\mathbb{Z})$ and $d\mu(z) = y^{-2} \, dx \, dy$.


This question has a nice answer, as does the same question where one allows any triple product involving various cusp forms and Eisenstein series. When at least one form is an Eisenstein series, this is the Rankin-Selberg method; when all three are cusp forms, this is known as the Watson-Ichino triple product formula. In Section 2.2 of [this paper of mine](https://doi.org/10.1007/s00208-018-1677-9), I write out the exact formulae for these in several cases, and provide some references.




---


Let me give a more detailed explanation of this in the case written above. It is known, going back to work of Hecke, that any Maass cusp form $f$ can be written as a finite linear combination of Hecke-Maass cusp forms; these are eigenfunctions of the Hecke operators and have Fourier expansions of the form
$$f(z) = \rho\_f \sum\_{n = 1}^{\infty} \frac{\lambda\_f(n)}{\sqrt{n}} W\_{0,it\_f}(4\pi ny) (e^{2\pi i nx} + \epsilon\_f e^{-2\pi i nx}).$$
Here $\rho\_f$ is a normalising constant, $\lambda\_f(n)$ is the $n$-th Hecke eigenvalue (which is multiplicative as an arithmetic function and conjecturally satisfies $\lambda\_f(p) \in [-2,2]$ for any prime $p$), $t\_f$ is the spectral parameter of $f$ (so that $f$ has Laplacian eigenvalue $\lambda\_f = 1/4 + t\_f^2$), $W\_{\alpha,\beta}$ denotes the Whittaker function (in this case, it is basically just a $K$-Bessel function), and $\epsilon\_f \in \{1,-1\}$ is the parity of $f$ (so that $f(-\overline{z}) = \epsilon\_f f(z)$).


So we may assume that $f\_1$ and $f\_2$ are of this form. One then *unfolds* the integral by inserting the identity
$$E(z,s) = \sum\_{\gamma \in \Gamma\_{\infty} \backslash \Gamma} \Im(\gamma z)^s,$$
using the fact that $f\_1(\gamma z) f\_2(\gamma z) = f\_1(z) f\_2(z)$, and that $\Gamma\_{\infty} \backslash \mathbb{H} \cong \{z : x \in [0,1], \ y > 0\}$. One then inserts the Fourier expansions of $f\_1$ and $f\_2$ and evaluates the $x$ integral, then the $y$ integral; this latter integral is essentially the Mellin transform of the product of two $K$-Bessel functions, which has an explicit expression in terms of products of gamma functions.


All in all, one finds that the triple product is equal to
$$\frac{\rho\_{f\_1} \rho\_{f\_2} \pi^{-s}}{\Gamma(s)} \prod\_{\pm\_1, \pm\_2} \Gamma\left(\frac{s \pm\_1 it\_{f\_1} \pm\_2 it\_{f\_2}}{2}\right) \sum\_{n = 1}^{\infty} \frac{\lambda\_{f\_1}(n) \lambda\_{f\_2}(n)}{n^s}.$$
(At least if $\epsilon\_{f\_1} = \epsilon\_{f\_2}$; otherwise, the gamma functions are changed slightly.) Here we are initially assuming that $\Re(s) > 1$; however, by analytic continuation, this extends meromorphically to all of $\mathbb{C}$ (with a simple pole at $s = 1$ if $f\_1 = f\_2$ and no poles otherwise). Let $\zeta(s) = \sum\_{n = 1}^{\infty} n^{-s}$ denote the Riemann zeta function. Then we *define* the Rankin-Selberg $L$-function $L(s,f\_1 \times f\_2)$ to be
$$L(s,f\_1 \times f\_2) = \zeta(2s) \sum\_{n = 1}^{\infty} \frac{\lambda\_{f\_1}(n) \lambda\_{f\_2}(n)}{n^s}.$$
Again, this is initially defined for $\Re(s) > 1$ in this way, but extends meromorphically to $\mathbb{C}$. Moreover, this has a nice functional equation, which can be proven using the functional equation for $E(z,s)$.




---


Next, let me discuss the same problem for
$$\int\_{\Gamma \backslash \mathbb{H}} f(z) E\left(z,\frac{1}{2} + it\right) E(z,s) \, d\mu(z).$$
The key point is that $E(z,1/2 + it)$ has a nice Fourier expansion; $t$ replaces $t\_f$, while $\lambda\_f(n)$ is replaced by $\lambda(n,t) = \sum\_{ab = n} a^{it} b^{-it}$. The exact same method works, except that the ensuing Rankin-Selberg $L$-function factorises as the product of two $L$-functions:
$$\zeta(2s) \sum\_{n = 1}^{\infty} \frac{\lambda\_f(n) \lambda(n,t)}{n^s} = \sum\_{m = 1}^{\infty} \frac{\lambda\_f(m)}{m^{s + it}} \sum\_{n = 1}^{\infty} \frac{\lambda\_f(n)}{n^{s - it}} = L(s + it,f) L(s - it,f).$$




---


Now we consider
$$\int\_{\Gamma \backslash \mathbb{H}}^{\mathrm{reg}} E\left(z,\frac{1}{2} + it\_1\right) E\left(z,\frac{1}{2} + it\_2\right) E(z,s) \, d\mu(z).$$
I have written $\int^{\mathrm{reg}}$ to denote the fact that we must *regularise* (or *renormalise*) this integral in order for it to converge. This goes back to [work of Zagier](https://doi.org/10.15083/00039589). Here not only does the Rankin-Selberg $L$-function factorise into the product of two $L$-functions, but each of these $L$-functions factorises further into the product of two copies of the Riemann zeta function.




---


Finally, let me consider
$$\int\_{\Gamma \backslash \mathbb{H}} f\_1(z) f\_2(z) f\_3(z) \, d\mu(z).$$
Here the Rankin-Selberg method does not apply, since we cannot unfold as there is no Eisenstein series present. Nonetheless, this triple product can be shown (via methods coming from representation theory) to be associated to $L$-functions. This was done in [Watson's thesis](https://arxiv.org/abs/0810.0425v3) and was further expanded upon by [Ichino](https://doi.org/10.1215/00127094-2008-052). What they show is an identity relating the *square of the absolute value* of this integral to a certain $L$-function, namely $L(1/2, f\_1 \times f\_2 \times f\_3)$.


