Your answer is a bit confusing, you should try to be more schematic... I mean, I think you got it almost right, but it lacks something, and thinking more schematically can and will help you.
This is what I mean:
1. You are looking at the electric field, so (though useful) you don't actually need the point charge q0.
Now, as you did implicitly, let's make some symmetrical analysis:
2. The situation is axial-symmetrical: so, we can calculate Ey of one charge and then (in this case) double it to consider the effect of the other one.
3. Now your turn: since you need the dependence of the electric field from the y-component, let's calculate Ey(y). Which is the basic equation? How do you get the y-component?
4. How do you calculate Ey depending only from y? (Here is the problem with your answer: your Ey depends on y, but also on teta! This makes taking the derivative a very complex task.)
5. You need to find the maximum... all right then, back to calculus! :D As you said, you need to find the Y so that dEy(Y)=0.
Hope it helps. You were almost there, so if you like you can just use your result (Ey = (2kQ2sin(θ)) / (q0 sqrt(a2+y2))) and consider only points 4 and 5 to go on.
Remember to replace sin(teta) with some y-dependent expression before deriving. This is the main thing to do.