I came across the notion of complex energy while studying instanton method to study the unstable state. Unstable states are those which have energy with an imaginary part. But as we know Hamiltonian is a Hermitian operator. So how can the energy be imaginary?

  • 1
    $\begingroup$ probably a calculation trick akin to here: physics.stackexchange.com/q/60185 - but I don't know the exact circumstances. $\endgroup$ – Martin May 26 '14 at 16:32
  • $\begingroup$ Are you still looking for further clarification? If so you can always comment on the existing answers. Don't forget to accept an answer once you are satisfied. $\endgroup$ – DanielSank May 29 '14 at 2:36

Since nobody is answering this I'll give it a whack.

A Hermitian Hamiltonian gives unitary evolution. The evolution operator is

$$ U(t) = \exp[-i\Omega t] $$

where $\Omega \equiv H/\hbar$. The thing in the exponent is anti-Hermitian because of the $i$ in front of the Hamiltonian. When you exponentiate an anit-Hermitian operator, you get a unitary one, as I'll explain now.

Let $|\Psi'\rangle= \exp[-ih]|\Psi\rangle$ and Let $|\Phi'\rangle= \exp[-ih]|\Phi\rangle$. We compute\begin{equation} \langle \Phi'|\Psi'\rangle = \langle\Phi|e^{ih}e^{-ih}|\Psi\rangle = \langle\Phi|\Psi\rangle. \end{equation} Since the operator $e^{ih}$ preserves the inner product, it is unitary by definition.

So, as you said Hermitian Hamiltonians give unitary evolution. The thing is, if you consider only sub-systems of a complete physical system, the evolution is not unitary. For example, if we're studying a particle in box, there could be a possibility of the particle leaving the box, in which case our assumption that the particle is in the box is wrong. You can account for this by manually inserting a decreasing function into the evolution operator

$$U(t) = e^{-t\Gamma}e^{-i\Omega t}.$$

Here $1/\Gamma$ is an exponential decay time constant which suppresses the magnitude of the wave function (or state vector, whatever). It's convenient to just write this as a single exponential:

$$U(t) = \exp[-i(\Omega - i\Gamma)t].$$

This can be interpreted as $\Omega\rightarrow\Omega-i\Gamma$ which looks like a non-hermitian Hamiltonian. In fact, if you take the point of view that the Hamiltonian is the thing that generates time evolution then the Hamiltonian is not Hermitian for this subsystem. This does not contradict the idea that the Hamiltonian of a physical system must be Hermitian, because we're doing a special case where we're dealing with a subsystem, and we're encapsulating the effect of the outside environment (into which our particle can escape) by using a non-Hermitian Hamiltonian.

of course, in the end, most systems under study are not completely isolated from other degrees of freedom, so in real life the notion of non-Hermitian Hamiltonian is extremely useful.


Imagine you have two 2-level systems interacting with one another through the following Hamiltonian (at this point I'm carelessly dropping dimensions from things):

$$ H = \left[ \begin{array}{cccc} 0&0&0&0 \\ 0&0&1&0 \\ 0&1&0&0 \\ 0&0&0&0 \end{array} \right]$$ where the order of basis vectors is $\{|00\rangle, |01\rangle, |10\rangle, |11\rangle\}$. Suppose this interaction is in place for a time $t$. Then the evolution operator is

$$U(t) = \left[ \begin{array}{cccc} 0&0&0&0 \\ 0&\cos(t)&-i\sin(t)&0 \\ 0&-i\sin(t)&\cos(t)&0 \\ 0&0&0&0 \end{array} \right].$$

Suppose we now "forget" that the second of the two 2-level systems exists, thus making it an "environment" of the first 2-level system. In this condition our remaining 2-level system cannot be represented by a state vector; we have to use a density matrix. Suppose the initial state is $$|\Psi(0)\rangle = \left[\begin{array}{c} 0\\1\\0\\0 \end{array}\right].$$ Then after time $t$ the state vector is $$|\Psi(t)\rangle = \left[\begin{array}{c} 0\\\cos(t)\\-i\sin(t)\\0 \end{array}\right].$$ The density matrix of the full system is $$\rho = \left[ \begin{array}{cccc} 0&0&0&0 \\ 0&\cos(t)^2&i\sin(t)\cos(t)&0 \\ 0&-i\cos(t)\sin(t)&\sin(t)^2&0 \\ 0&0&0&0 \end{array} \right].$$ To find the density matrix of the subsystem consisting of just the remaining 2-level system, we trace over the "lost" particle's subspace, which leaves $$\rho_{\text{subsystem}} = \left[ \begin{array}{cc} \sin(t)^2&0 \\ 0&\cos(t)^2 \end{array}\right]$$ where here the states are ordered $\{|0\rangle, |1\rangle\}$.

First of all, note that you cannot get this evolution from a normal Hamiltonian. The state is traveling from one pole of the Bloch sphere to the other, going through the center of the sphere. This is impossible for a pure state under unitary evolution. Also note that the density matrix is always diagonal which means that the information we have about our subsystem is always classical. This is typical of cases where we are missing a piece of an entangled system (in fact it is argued that this is the essence of the quantum to classical transition).

This is a demonstration of how you get non-unitary evolution if you look only at a subsystem of the full quantum system.

  • $\begingroup$ the interesting physical question that remains is how Hamiltonians of open, perturbed systems become non-Hermitian $\endgroup$ – diffeomorphism May 27 '14 at 2:17
  • $\begingroup$ @diffeomorphism: So you want an explicit demonstration of how the evolution of a subsystem can be described with a non-Hermitian Hamiltonian? How about the type of situation that leads to Fermi's golden rule. That's a pretty good case where you can compute things explicitly. $\endgroup$ – DanielSank May 27 '14 at 4:50
  • $\begingroup$ I was thinking more at the algebraic Hilbert operator level. To my understanding the Fermi golden rule is kind of heuristic reasoning using the map of experimental eigenvalues to obtain transition rates $\endgroup$ – diffeomorphism May 27 '14 at 4:54
  • 1
    $\begingroup$ @diffeomorphism: Fermi's golden rule can be derived without any hand waving if you assume your environment dissipates information quickly compared to the interaction strength of your system and the environment. Also, this notion you have of representing environments as non-Hermitian components of the Hamiltonian is very well known. Look at the Lindblad master equation, or any of the so-called "quantum noise" literature. $\endgroup$ – DanielSank May 27 '14 at 6:01
  • 1
    $\begingroup$ @diffeomorphism To be explicit, you are asking about the quantum trajectories formalism, also called quantum Monte Carlo or quantum jumps. This allows you to derive, under a specific set of assumptions, a non-Hermitian effective Hamiltonian dynamics, which describes the evolution of the wave function conditioned on whether a transition (e.g. a la Fermi's golden rule) has occurred or not. This is the canonical review article for the topic, but there is a huge literature that you can find easily by Googling. $\endgroup$ – Mark Mitchison May 27 '14 at 9:50

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.