Is it possible to recover Classical Mechanics from Schrödinger's equation? Let me explain in details. Let $\Psi=\Psi(x,t)$ be the wave function of a particle moving in a unidimensional space. Is there a way of writing $\Psi(x,t)$ so that $|\Psi(x,t)|^2$ represents the probability density of finding a particle in classical mechanics (using a Dirac delta function, perhaps)?
 A: @Arnoques
Sorry, but I think there is an error in your answer: 
The spatial extent of the particle wave-function, must be much smaller (and not longer) than the variation length-scale of the potential, to transform $\langle \nabla V(x)\rangle$ turns into $\nabla V\left(\langle x\rangle \right).$
Only in this case, it is possible to make a Taylor series of $V(X)),$ because $V(X)$ is slowly varying in the domain where the wave function is not null, and you can take the mean expectation : 
$$\nabla{\mathrm  i}\, V(X)  = \nabla{\mathrm i}\, V(\langle X\rangle )  + (X_j - \langle X_j\rangle ) \nabla j \,\nabla{\mathrm i} \,V(\langle x\rangle )\, +\,\textrm{negligible higher order terms in}\,\, (X_j - \langle X_j\rangle)$$
So,  $\langle \nabla{\mathrm i}\, V(X)\rangle  = \nabla{\mathrm i} \,V(\langle X\rangle),$ because   $\langle X_j - \langle X_j\rangle \,\rangle = 0\;.$
A: You can recover Schroedingers equation from the path integral formulation of Quantum mechanics by Feynman. In the path integral picture the classical trajectories are the stationary points of the integrand. So in the stationary phase approximation, they are the contribution of $0$-th order in $\hbar$. Of course that is not a direct relation between the Schroedinger equation and classical trajectories. 
A: Sure you can! This is actually a simple but very interesting result, and it is usually shown in quantum mechanics courses. It's called the Ehrenfest theorem, and I won't prove it here but I'll copy the result from Sakurai Modern Quantum Mechanics (1991). You can check the mathematical details there, or in many other books.
If you have a hamiltonian with the form
$$H = \frac{p^2}{2\,m}+V(x)$$
you can prove that, in the Heisenberg picture, 
$$m \frac{\mathrm{d}^2x}{\mathrm{d}t^2} = -\nabla V(x) .$$
If you now take the expectation value of that equation (for certain state kets), you get
$$m \frac{\mathrm{d}^2\langle x \rangle}{\mathrm{d}t^2} = \frac{\mathrm{d}\langle p \rangle}{\mathrm{d}t} = -\langle \nabla V(x) \rangle .$$
This result is valid in both Heisenberg and Schrödinger's picture.
If you want to recover the classical limit, you need to say that the area where the wavefunction is significantly nonzero is much smaller than the scale of variations of the potential. In that case, you can identify the center of the wavefunction with the position of the particle, and $\langle \nabla V(x) \rangle $ turns into $\nabla V(\langle x \rangle) $.
What this means, conceptually, is that the center of the wavefunction will move according to the classical laws if you can't "see" that your object/particle it's not a material point, and if your potential is also classical, in that it doesn't have variations that are comparable to the "size" of the wavefunction.
A: The short answer: No, does not exist any wavefunction in Hilbert space which reproduces classical mechanics.
The classical limit of quantum mechanics is studied with some deep in Ballentine textbook. For instance, section 14.1 is devoted to the Ehrenfest theorem and it is shown that the theorem is neither necessary nor sufficient to define the classical
regime. 
The paper What is the limit $\hbar \rightarrow 0$ of quantum theory? (Accepted for publication in the American Journal of Physics) shows that Schrödinger's equation for a single particle moving in an external potential does not lead to Newton's equation of motion for the particle in the general case. Page 9 of this more recent article precisely deals with the question of why no wavefunction in the Hilbert space can give a classical delta function probability.
A: More intuitive picture is in Arnoques answer, alternative and a bit more formal approach is to note that all QM equations of motion have their classical mechanics equivalent if you formulate them using commutators and then replace commutator with Poisson bracket ($\partial A/\partial t = [H,A]$ $\Rightarrow$ $ \partial a/\partial t = \{ H,a \}_{q,p} $, if you "hide" Planck constant). The commutator itself is of course zero in classical case, when operators reduce to numbers. Accordingly, all general system properties easily map from QM to CM. 
And it may be shown (too much to copy, sorry) that a formal limit $\hbar\to0$ leads to exact equivalence between commutator and Poisson bracket. 
Concerning the wavefunction, classical motion is definite. Instead of probability you have definite correspondence between $t$ and $x$. Indeed, you may formulate it saying that classical $|\Psi(x,t)|^2=\delta(x-x_c(t))$ where $x_c(t)$ is classical tragectory. To write $\Psi$ itself, you have to treat $\hbar\to0$ accurately to avoid divergent integrals. Normally, there is no reason to do this. And technically, there is no guarantee that a limit $\hbar\to0$ of some particular quantum state is a "normal" solution of classical problem.
A: This is exactly what Feynman's path-integral does, many-body quantum effects reveals classical limit at high temperature! Instead of solving the Schrödinger equation, it solves Newton's equation by splitting the atom into beads and perturbing the system effectively to find out the thermal equilibrium state of the system. The method showed the isomorphism of quantum theory and classical statistical mechanics, which leads to an interesting point: the imaginary time is isomorphic to the inverse of temperature (see Wick-rotation).
You may want to read Feynman's thesis.
BTW, I personally started the topic with reading this paper:
Feynman's derivation of the Schrödinger equation 
