Register to reply

Problem with Hermicity of Density matrix

by climbon
Tags: density, hermicity, matrix
Share this thread:
May1-12, 07:07 AM
P: 18
Hi, I am trying to solve a modified Jayne's cummings model using the Von Neumann equation and Wigner function but am having a problem with Density matrix hermicity;

I am trying to solve in Schrodinger picture.

I have my system Hamiltonian as;

H_{0} = \frac{1}{2}\hbar \Omega \sigma_{z} + \hbar \omega aa^{\dagger}

I am having a problem with the second term of this system Hamiltonian...when I follow it through into the Von neumann equation and Wigner function I get a non-Hermitian matrix??

I can't find a problem with the maths...the only thing I can think is to drop the second term and leave it out of the Von Neumann equation...if this is possible, why?

Phys.Org News Partner Physics news on
UCI team is first to capture motion of single molecule in real time
And so they beat on, flagella against the cantilever
Tandem microwave destroys hazmat, disinfects
May1-12, 07:55 AM
Sci Advisor
P: 2,546
In the Schroedinger picture of time evolution the full time dependence is on the states, and thus the (picture-independent) von Neumann equation, translates into

[tex]\frac{\mathrm{d}}{\mathrm{d} t} \hat{R}=\frac{1}{\mathrm{i} \hbar} [\hat{R},\hat{H}]+ \left (\frac{\partial \hat{R}}{\partial t} \right )_{\text{expl}}.[/tex]

Since [itex]\hat{R}[/itex] is Hermitian at the initial time and [itex]\hat{H}[/itex] is Hermitian, [itex]\hat{R}[/itex] is Hermitian at any time.

So there must be some mistake in your calculation, because your Hamiltonian is obviously Hermitian. To find, what's wrong, you should post your concrete problem and your ansatzes for the solution.
May18-12, 11:13 AM
P: 18
Thanks, this is my problem;


H=H_{0} + H_{I}



H_{0}= \frac{1}{2} \hbar \Omega \sigma_{z} + \hbar \omega a^{\dagger}a


H_{I}=g_{1} \sigma_{x} (a+a^{\dagger}) + g_{2} \sigma_{z} (a+a^{\dagger}) +g_{3} \sigma_{z} (a^{2} + aa^{\dagger} + a^{\dagger}a +a^{\dagger 2})

I put this into the von neumann equation;

i\hbar \frac{\partial}{\partial t} \rho = \left[ H_{0} + H_{I} , \rho \right]

Using the density matrix as a 2x2 matrix with elements [itex] \rho_{11} , \rho_{12}, \rho_{21}, \rho_{22} [/itex].

I then put all the equations into the von Neumann equation and multiply out all of the matrices and split the density matrix into its component parts. I then convert this from being a function of a and [itex] a^{\dagger} [/itex] into the Wigner function of x and y (using standard conversion and correspondences). The Wigner function should still be Hermitian. Below are the four matrix elements of the Wigner function;

i \hbar \frac{\partial}{\partial t} W_{11} = \hbar \omega \left(-1 + ix \frac{\partial}{\partial y} - iy \frac{\partial}{\partial x} \right) W_{11}
+g_{1} \left( 2x + \frac{1}{2} i \frac{\partial}{\partial y} \right) W_{21}
-g_{1} \left( 2x - \frac{1}{2} i \frac{\partial}{\partial y} \right) W_{12}
+g_{2} i \frac{\partial}{\partial y} W_{11}
+g_{3} \left( 4ix \frac{\partial}{\partial y} \right) W_{11}

i \hbar \frac{\partial}{\partial t } W_{12} = \hbar \Omega W_{12} + \hbar \omega \left( -1 + ix \frac{\partial}{\partial y} - iy \frac{\partial}{\partial x} \right) W_{12}
+g_{1} \left( 2x + \frac{1}{2} i \frac{\partial}{\partial y} \right) W_{22}
-g_{1} \left( 2x - \frac{1}{2} i \frac{\partial}{\partial y} \right) W_{11}
+g_{3} \left( 8x^{2} - \frac{1}{2} \frac{\partial^{2}}{\partial y^{2}} \right) W_{12}

i \hbar \frac{\partial}{\partial t } W_{21} = - \hbar \Omega W_{21} + \hbar \omega \left( -1 + ix \frac{\partial}{\partial y} - iy \frac{\partial}{\partial x} \right) W_{21}
+g_{1} \left( 2x + \frac{1}{2} i \frac{\partial}{\partial y} \right) W_{11}
-g_{1} \left( 2x - \frac{1}{2} i \frac{\partial}{\partial y} \right) W_{22}
-g_{3} \left( 8x^{2} - \frac{1}{2} \frac{\partial^{2}}{\partial y^{2}} \right) W_{21}

i \hbar \frac{\partial}{\partial t} W_{22} = \hbar \omega \left(-1 + ix \frac{\partial}{\partial y} - iy \frac{\partial}{\partial x} \right) W_{22}
+g_{1} \left( 2x + \frac{1}{2} i \frac{\partial}{\partial y} \right) W_{12}
-g_{1} \left( 2x - \frac{1}{2} i \frac{\partial}{\partial y} \right) W_{21}
-g_{2} i \frac{\partial}{\partial y} W_{22}
-g_{3} \left( 4ix \frac{\partial}{\partial y} \right) W_{22}

Which are non-Hermitian. I have checked the calculation several times and cannot find a mathematical error. I can only think I've omitted something. For example, do I need to apply some unitary operator somewhere, eg;

U(t,t_{0})= e^{i \hbar H_{0}(t_0)}

Am I missing something?

Thanks for any help, is much appreciated.

Register to reply

Related Discussions
Problem with density matrix Quantum Physics 19
Conditions for a density matrix; constructing a density matrix Introductory Physics Homework 5
Charge Density from Density Matrix? Quantum Physics 0
Testing Hermicity Quantum Physics 2
Hermiticity and expectation value Quantum Physics 6