I searched quite a lot textbooks to answer this question, but I couldn't find anything satisfactory. This could be more general question, but it will be easier to just restrain to particle (say electron) in the infinite potential well.
When it comes to the general solution of Schroedinger equation, the wavefunction is a superposition of many stationary states:
$$\Psi(x,t) =\sum\limits_{n}{c_n\psi _n(x)e^{\frac{-iEt}{ \hbar } }}$$
And from this point, textbooks tend to continue to method of calculating the values of the $c_n$ coefficients by using Fourier transform and given $\Psi(x,t=0)$.
My question is, is it possible, to calculate these values from first principles, i.e. using some other related theories? I can imagine the task like: Given that electron with energy $E_e$ is in the one-dimensional infinite potential well of length $L$, give the probability that measured energy will correspond to the energy of lowest possible state in this box.
I doubt that when modeling such a electron, it is neccessary to make up all the c-coefficients, and just adjust it to have $\sum\limits_{n}{c_n^2}=1$ and $\sum\limits_{n}{c_n^2E_n}=E_e$. I suppose that there have to be some theory that is capable of providing answers. I was considering Fermi-Dirac distribution, which in principle seems to be created for this, but I couldn't find a straightforward answer (which probably do not exist), and I know too little to make use of it.