The average $\langle E\rangle$ you have on the right is not to be understood in the sense of statistical mechanics (and thus as a temperature).
The wavefunction $\psi$ can be a linear combination of states $\psi_n$ of definite energy of your molecule. The $\psi_n$’s are solutions of the Schrodinger equation. Thus,
\begin{align}
\psi(x,t)&=\sum_i a_i \psi_i(x,t)\, ,\qquad \sum_i \vert a_i\vert^2=1\, , \tag{1} \\
\hat H\psi_i(x,t)&=E_i\psi_i(x,t)\, .
\end{align}
With this wavefunction the average energy
$$
\langle E\rangle = \sum_i \vert a_i\vert^2 E_i \tag{2}
$$
is a weighted average of the possible energies of your molecule, with the modulus square $\vert a_i\vert^2$ of the complex amplitude $a_i$ functioning as statistical weight. $\psi(x,t)$ given in (1) is an example of a pure state.
For instance, if you have a hydrogen wavefunction of the form
$$
\psi(r)=a_1\psi_{100}(r)+a_2\psi_{200}(r)
$$
then the average energy of this system is
$$
\vert a_1\vert^2\times (-13.6) + \vert a_2\vert^2 \times (-13.6/4). \tag{3}
$$
Measuring the energy, you will sometimes get the value $-13.6$ (this outcome will occur $\vert a_1\vert^2$ of the time) and sometimes get the value $-13.6/4$ (this outcome will occur $\vert a_2\vert^2$ of the time). The average energy is exactly given by (3).
In quantum statistical mechanics one introduces the concept of a mixted state. For mixed states one cannot define a wavefunction for the system as in (1). The system is described using a matrix that represents a statistical mixture of wavefunction that is not of the form given in (1); the average energy for this statistical mixture is also not of the form (2) since the latter comes from (1), which does not exist for mixted states.