Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Probability Density and Current of Dirac Equation

  1. Feb 4, 2013 #1
    Hey,

    I'm trying to determine the probability density and current of the Dirac equation by comparison to the general continuity equation. The form of the Dirac equation I have is

    [tex]i\frac{\partial \psi}{\partial t}=(-i\underline{\alpha}\cdot\underline{\nabla}+\beta m)\psi[/tex]

    According to my notes I am supposed to determine the following sum to make the relevant comparisons to the continuity equation and therefore determine the probability density/current

    [tex]\psi(Dirac)^{\dagger}+\psi^{\dagger}(Dirac)[/tex]

    Where 'Dirac' refers to the above equation. However I have tried this and I can only get it to work if I multiply one term by 'i' and the other by '-i' in the above.

    [tex]\psi(i\frac{\partial \psi}{\partial t})^{\dagger}+\psi^{\dagger}(i\frac{\partial \psi}{\partial t})=-i\psi\frac{\partial \psi^{*}}{\partial t}+i\psi^{*}\frac{\partial \psi}{\partial t}\neq i\frac{\partial (\psi^{*}\psi)}{\partial t}[/tex]

    Any help is appreciated!

    Thanks,
    SK
     
  2. jcsd
  3. Feb 4, 2013 #2

    dextercioby

    User Avatar
    Science Advisor
    Homework Helper

    Your notes should use the gamma matrices. It's the modern treatment on the Dirac equation/field and make the special relativity invariance easier to see.
     
  4. Feb 4, 2013 #3
    Perhaps they should but I'm reckoning these are introduced later, so considering we don't 'know' these yet this appears to be the simplest way of demonstrating the probability density/current of the Dirac equation. I'm confused though, I perhaps maybe taking the adjoint of the Dirac equation incorrectly.
     
  5. Feb 4, 2013 #4
    Have you considered that taking the hermitian conjugate is not only taking the complex conjugate but also the transposition?
     
  6. Feb 4, 2013 #5
    I think so, taking the adjoint of ψ doesn't bring out a minus sign does it? With regards to the RHS of the Dirac equation I think β is diagonal and so the transposition doesn't affect it, though I'm a bit confused as how I'd go about doing the hermitian conjugate on the dot product of the alpha matrix with the ∇...

    I wouldn't be surprised though if this transposition is the issue, I'll keep looking at it!
     
  7. Feb 4, 2013 #6
    Also, [itex]\psi^{\dagger}\psi[/itex] is a number, while [itex]\psi\psi^{\dagger}[/itex] is a matrix, so I don't really quite get the whole thing.
     
  8. Feb 4, 2013 #7

    stevendaryl

    User Avatar
    Staff Emeritus
    Science Advisor

    The complicated thing is correctly taking the Hermitian conjugate of the Dirac equation. What I think is true is this:

    1. [itex] i \dfrac{d\Psi}{dt} = -i \alpha \cdot (\nabla \Psi) + \beta m \Psi[/itex]
    2. [itex] -i \dfrac{d\Psi^\dagger}{dt} = +i (\nabla \Psi^\dagger \cdot \alpha) + \Psi^\dagger \beta m[/itex]
    Taking the conjugate reverses the order of matrices. So if you multiply the top equation on the left by [itex]-i \Psi^\dagger[/itex] and multiply the bottom equation on the right by [itex]+i \Psi[/itex] and add them, you get:
    [itex]\Psi^\dagger \dfrac{d\Psi}{dt} + \dfrac{d\Psi^\dagger}{dt}\Psi= - \Psi^\dagger \alpha \cdot (\nabla \Psi) - (\nabla \Psi^\dagger) \cdot \alpha \Psi[/itex]

    (the terms involving [itex]\beta[/itex] cancel). You can rewrite this as (I think):

    [itex]\dfrac{d}{dt} (\Psi^\dagger \Psi) = - \nabla \cdot (\Psi^\dagger \alpha \Psi)[/itex]

    This can be rearranged as a continuity equation for probability.

    There's a different continuity equation for electric charge, but I've forgotten what that is.
     
  9. Feb 5, 2013 #8
    Thanks Stevendaryl!

    This is exactly what I wanted, I can see what was incorrect now - brilliant!

    Thanks again,
    SK
     
  10. Feb 5, 2013 #9

    stevendaryl

    User Avatar
    Staff Emeritus
    Science Advisor

    My last comment was stupid. If you multiply probability current by the electron charge, you get the charge current. It's that simple. For some reason, though, that's not the case with the solutions of the Klein Gordon equation.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: Probability Density and Current of Dirac Equation
Loading...