10

How does one decide whether a wave function is a physically acceptable solution of the Schrödinger equation? For example: $\tan x$ , $\sin x$, $1/x$, and so on.

Qmechanic
  • 201,751
  • 4
    You mean besides sticking it into the Schrodinger equation and seeing if it is a valid solution? – Kyle Kanos Nov 26 '14 at 15:00
  • 4
    http://physicspages.com/2011/01/25/wave-function-borns-conditions/ contains a nice summary of the conditions you're looking for. Most of the conditions are a consequence of interpreting the wavefunction as a probability amplitude. – Endulum Nov 26 '14 at 15:04
  • 4
    Shove it into the Schrodinger equation, and try to solve for the potential $V$; that will give you some idea of what system it is describing. – JamalS Nov 26 '14 at 15:09
  • 1
    Also, note that physically unacceptable solutions of the Schrödinger equation can still be useful. For example, the eigenfunctions of the position operator (delta function) and momentum operator (plane wave) cannot be created experimentally. But it is still very useful to describe the physically realizable states in terms of these eigenfunctions. – jabirali Nov 26 '14 at 16:14

2 Answers2

15

The very minimum that a wavefunction needs to satisfy to be physically acceptable is that it be square-integrable; that is, that its $L_2$ norm, $$ \int |\psi(x)|^2\mathrm d x, $$ be finite. This rules out functions like $\sin(x)$, which have nonzero amplitude all the way into infinity, and functions like $1/x$ and $\tan(x)$, which have non-integrable singularities.

In the most rigorous case, however, one needs to impose additional conditions. The physically preparable states of a particle denote functions which are continuously differentiable to any order, and which have finite expectation value of any power of position and momentum. Thus:

  • $\psi$ must be continuous everywhere.
  • All of $\psi$'s derivatives must exist and they must be continuous everywhere.
  • The expectation value $\int\psi^*(x) \:\hat x^n \hat p^m \psi(x)\:\mathrm dx$ must be finite for all $n$ and $m$.

This rules out discontinuous functions like $\theta(x)$, functions with discontinuous derivatives, and functions like $(1+x^2)^{-1/2}$, which decay too slowly at infinity. States which satisfy these conditions are called 'physical' because they are the states that can be prepared with finite energy in finite time. The way to implement these states rigorously is by using a construction known as a Rigged Hilbert Space (see also Galindo & Pascual's QM textbook).

In everyday practice, most people adopt a bit of a mixed approach. The requirements that a function be continuous is never dropped, and one requires it to be differentiable at least almost everywhere. If the Hamiltonian is not a nice function of position, though, such as with $\delta$-function or square-well potentials, the requirements are sometimes slackened to only those; this is in the understanding that a truly discontinuous potential is not physical, and that any problems brought into the higher derivatives of $\psi$ can be fixed by using a smoother Hamiltonian.

Emilio Pisanty
  • 132,859
  • 33
  • 351
  • 666
  • Realize that for the condition of finite energy, one just needs the restriction that the second momentum moment $\langle p^2\rangle_\psi$ is finite, not all position and momentum moments. – Mateus Sampaio Nov 26 '14 at 15:53
  • @MateusSampaio The statement is about the preparation of the state, using a Hamiltonian of the form $\hat H=H(\hat x,\hat p,t)$, within finite time and starting from a reference physical state. (You would also ask for all moments of $\hat H$ to be finite if the function $H$ was not well-behaved enough.) – Emilio Pisanty Nov 26 '14 at 16:00
  • 1
    I don't see why the preparation of the state would require all moments of $H$ to be finite. – Mateus Sampaio Nov 26 '14 at 16:03
  • Why should a physical state be on the nuclear space of the rigged hilbert space? I can't see any physical reason for that – Hydro Guy Nov 26 '14 at 16:16
  • @MateusSampaio I'm trying to find references and will post them here when I find a proof. I probably first read the statement here and in related papers. – Emilio Pisanty Nov 26 '14 at 22:08
  • Your second rule (continuity of all derivatives) rules out hydrogen S-orbitals... – Ruslan Nov 27 '14 at 04:04
  • @Ruslan You mean at the point where the Hamiltonian is singular? – Emilio Pisanty Nov 27 '14 at 12:05
  • Yes, the point of electron-nucleus collision. – Ruslan Nov 27 '14 at 14:38
  • @Ruslan Real nuclei (even bare protons) have extended charge distributions, so you can repair the problem by using a more physical model of the nuclear charge distribution. Not that is solves the problem for positronium, though. I guess you need to appeal to a more sophisticated framework for that. – dmckee --- ex-moderator kitten Jun 07 '18 at 00:59
  • Isn't $\sin(x)$ a valid solution to the equation, since $\int \sin^2x\ dx$ is finite and bounded? – Aniruddha Deb Dec 02 '20 at 15:06
  • @AniruddhaDeb $\int\sin^2(x) \mathrm dx$ is an unbounded function over $\mathbb R$. If you're confused about the details, ask separately -- the comments thread of a six-year-old answer is not the place for it. – Emilio Pisanty Dec 02 '20 at 15:59
6

If you are talking about the time independent Schrödinger equation, it's not a trivial question as it may seem, as the comments suggest. I will restrict the answer to the one-dimensional case, since multiply connected domains in higher dimensions give some additional problems. Not all functions $\psi$ that are solutions of the equation $$-\frac{\hbar^2}{2m}\psi''+V\psi=E\psi$$ are valid ones. The first condition is that $\psi\in L^2(\Omega)$, where $\Omega\subset \Bbb{R}$ is the domain of the function, since it must be an element of the Hilbert space, otherwise it would not be a quantum state.

More subtle conditions are required when you look at the domain of the Hamiltonian $$H=-\dfrac{\hbar^2}{2m}\dfrac{d^2}{dx^2}.$$ In general, this will depend on the conditions satisfied by the potential $V(x)$. Usually one ends up with subsets of the Sobolev space $\mathcal{H}^2(\Omega)$, which restricts the original space to functions such that their second order (weak-)derivative is in $L^2(\Omega)$. If $\Omega$ is an interval (which is the usual setup) this can also be put in an equivalent way as the functions that are, along with their derivative, absolutely continuous and whose second derivative is also in $L^2(\Omega)$. Also, when the domain $\Omega$ is a proper subset of $\Bbb{R}$, the boundary conditions, which are established via physical arguments, play a decisive role in choosing the right domain $\mathcal{D}(H)$ of self-adjointness, and so they must be considered too. For example, the condition that $\psi(0)=\psi(a)=0$ for the infinite square well, rules out some solutions of the Schrödinger equation that would satisfy the other conditions.

If we look otherwise to the time dependent Schrödinger equation $$H\psi(t)=i\hbar\partial_t\psi(t), \qquad \psi(0)=\psi_0$$ any function $\psi_0\in L^2(\Omega)$ can be an initial condition for the system. But for $\psi_0 \notin\mathcal{D}(H)$ the trajectory given by $\psi(t)=e^{-iHt/\hbar}\psi_0$ is only a weak solution, in the sense that it's not a differentiable path and that the average energy is not defined (can be considered infinite) for all $t$, so these solutions can be considered non-physical.

Ellie
  • 4,089