I understand that the Mellin transform of a modular form is expected to satisfy RH when it is an eigenform of all Hecke operators, in which case it has an Euler product. Now about when the form is not an eigenform: Is it known a case where the zeros are all in the critical strip?



On the location of zeros of L functions from modular forms
When it is not a Hecke-Eigen form, the Hecke L-series connected with the modular form does not have an Euler product. However it can still be written as a linear combination of Hecke L-series that have Euler-products. Thus the situation will resemble the case of linear combinations of Dirichlet L-series. In particular we can use joint Voronin universality to obtain $\gg T$ zeroes in any strip $1/2<\sigma\_1<\Re(s)<\sigma\_2 < 1$, $-T < \Im(s)< T.$ These results follows from the analytic properties of the Rankin-Selberg zeta-function (which gives "independence results" for Hecke eigenvalues of primes attached to different cusp forms). You should be able to find results of this kind in Steuding's SLN "Value Distribution of L-functions".


Also the results of Davenport-Heilbronn for the Hurwitz zeta-function can be proved, i.e. For any $\epsilon>0$ there exists $\gg T$ zeroes with $1<\Re(s)<1+\epsilon$ and $-T < \Im(s) < T $.


Stronger results corresponding to results of Karatsuba, Bombieri, Hejhal and Selberg for Dirichlet L-function that holds close to the critical line likely also holds. I think the russian school (Irina Rezvyakova) has proved results in this direction.


