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
Vibrational motion of a single molecule measured in real time
Researchers demonstrate ultra low-field nuclear magnetic resonance using Earth's magnetic field
Bubbling down: Discovery suggests surprising uses for common bubbles
May1-12, 07:55 AM
Sci Advisor
P: 2,438
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