Physics Stack Exchange is a question and answer site for active researchers, academics and students of physics. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

Suppose I have a wave function $\psi(x)$ in position basis. I can make a density function by simply multiplying $\psi(x)$ and its conjugate $\psi^*(x)$. If I operate the density matrix $\rho(x,y)=\psi(x)\psi^*(y)$ with Hamiltonian $\hat{H}$ I will get, \begin{equation} \hat{H}\rho(x,y)=\hat{H}\psi(x)\psi^*(y)=\psi^*(y)\hat{H}\psi(x)+\psi(x)\hat{H}\psi^*(y) \end{equation} since the hamiltonian is a diferential operator. At the same time, If I take discrete basis representation, I will get $\langle n|\psi\rangle=V_{1}$ and $\langle\psi| n\rangle=V_{1}^{\dagger}$ which are the column vector and row vector respectively. Now if I operate the Hamitonian matrix $\hat{H}$ on the density matrix I will get, \begin{equation} \hat{H}\rho=H (V_{1} V_{1}^{\dagger})=(H V_{1}) V_{1}^{\dagger} \end{equation} That means, I need to apply the hamiltonian operator only once on a vector. Why it is not distributed like the differential operator.

\begin{equation} \hat{H}\rho= V_{1} H V_{1}^{\dagger}+(H V_{1}) V_{1}^{\dagger} \end{equation} Why this equation is not valid in the case of discrete basis ?

share|cite|improve this question
I am not sure if your $\hat{H} \rho(x,y)$ and acting $\hat{H}$ on both $\psi(x)$ and $\psi^*(y)$ is a proper manipulation. Take a pure state density operator for instance, $\rho = | \psi \rangle \langle \psi |$. $\hat{H} \rho = \hat{H} | \psi \rangle \langle \psi |$. The coordinate representation is $ \langle x |\hat{H} | \psi \rangle \langle \psi | x' \rangle $, which means the Hamiltonian only acts on $\psi(x)$, not the conjugated one. – user26143 Nov 22 '13 at 16:30
user26143 is directing you in good direction. In order to apply $\hat{H}$ to functions $\psi$, we need to have explicit definition of $\hat H$. For harmonic oscillator, $\hat H = -\frac{\hbar^2}{2m}\frac{ \partial^2}{\partial x^2} + \frac{1}{2}m\omega^2 x^2$. So when we have function of two variables like $\psi(x)\psi(y)$, the Hamiltonian operator will act only on the part that depends on $x$, treating the rest as a constant: $\hat{H} \left[\psi(x)\psi(y)\right] = \left[\hat{H}\psi(x)\right]\psi(y)$. – Ján Lalinský Nov 22 '13 at 18:05
Ok I got it ...Thank you user26143 and Jan Lalinsky – Sijo Joseph Nov 24 '13 at 18:06

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.