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

Trace in statistical mechanics

  1. Aug 20, 2005 #1
    Hi, im just starting a 3rd year course in Statistical Mechanics, and am a bit confused about the operator trace, Tr. I know there is a trace for quantum operators, as well as one in classical physics, but i am not sure how to calculate either, or their physical meaning. Any help would be great, thanks.

  2. jcsd
  3. Aug 20, 2005 #2
    Hello rayveldkamp,

    The trace of an operator [tex]\hat{B}[/tex] is defined as:

    [tex]Tr(\hat{B}) = \sum_{i=1}^{N} <\Psi_{i}| \hat{B} |\Psi_{i}>[/tex]

    where the [tex] |\Psi_{i} >[/tex] are a complete basis.
    You so to say take the sum of the diagonal matrix elements.

    Traces are used for example in context with density operators.
    You have for example many particles that are in different states.
    Take for example an ensemble of electrons (ensemble = many particles)
    or atoms in a hot oven.

    This ensemble can be described by the density operator:

    [tex] \hat{\rho} = \sum_{i=1}^{N} p_{i} |\Psi_{i}> <\Psi_{i}| [/tex]

    where [tex] p_{i} [/tex] is the probability (classical probability) of having a particle in the state [tex] |\Psi_{i}> [/tex].

    Ensemble of electrons, [tex]p_{1} =3/4[/tex] with spin up and [tex]p_{2}=1/4[/tex] spin down:

    [tex] \hat{\rho}} = \frac{3}{4} |\uparrow\rangle \langle\uparrow| + \frac{1}{4} |\downarrow \rangle \langle\downarrow| [/tex]

    Suppose you want to know what's the mean value of the ensemble for a certain operator [tex] \hat{A}[/tex] (for example spin, momentum etc).

    This can be done by calculating the trace:

    [tex] <\hat{A}> = Tr(\rho \hat{A}) = \sum_{i=1}^{N} < \Psi_{i}| \hat{\rho} \hat{A} |\Psi_{i}>[/tex].
  4. Aug 20, 2005 #3
    Thankyou for that, that clears up a lot of my misunderstanding. In lectures it's been mentioned that trace is the same as summing the diagonal elements of a matrix, however what matrix are we referring to here? So the density operater describes the probabilities of certain states, how is this related to the density matrix, and what is it exactly?

  5. Aug 20, 2005 #4
    Let's take for example a matrix [tex] M [/tex]

    [tex] M = \left( \begin{array}{cc}
    a & b \\
    c & d \\
    \end{array} \right) [/tex]

    Often the matrix is just denoted as [tex]M_{lm}[/tex],
    with [tex] M_{11} = a[/tex], [tex] M_{12} = b[/tex], [tex] M_{21} = c[/tex] and [tex] M_{22} = d[/tex].
    (l=1 or 2, m=1 or 2)

    Now what's a matrix of an operator in QM?
    The matrix is just defined as

    [tex]A_{lm} = \langle \Phi_{l}|\hat{A}|\Phi_{m} \rangle [/tex]

    So if someone asks you to write down the matrix representation
    of the operator [tex] \hat{A} [/tex], you just have to calculate the matrix elements [tex]A_{lm}[/tex] from above and then write them in a matrix form,
    for example:

    [tex] A_{lm} = \left( \begin{array}{cc}
    \langle \Phi_{1}|\hat{A}|\Phi_{1} \rangle & \langle \Phi_{1}|\hat{A}|\Phi_{2} \rangle \\
    \langle \Phi_{2}|\hat{A}|\Phi_{1} \rangle & \langle \Phi_{2}|\hat{A}|\Phi_{2} \rangle \\
    \end{array} \right) [/tex]

    An example is the matrix form of the [tex] \hat{S}_{z}[/tex] operator, which describes a spin easurement along the z-axis. You surely have heard of the Pauli spin matrices. (Try to calculate the matrix representation for[tex] S_{z}[/tex] or look it up in the Quantum Mechanics book by Cohen Tannoudji, where I found it well described).

    The density matrix is the matrix form of the density operator, namely

    [tex] \rho_{lm} = \langle \Phi_{l}|\hat{\rho}|\Phi_{m} \rangle [/tex]

    We know that

    [tex] B_{lm} = \langle \Psi_{l}| \hat{B} |\Psi_{m} \rangle [/tex]

    What's [tex] B_{ii} [/tex]?
    Answer: [tex] B_{ii} = \langle \Psi_{i}| \hat{B} |\Psi_{i} \rangle [/tex]

    Next step:
    The definition of the trace of an operator [tex]\hat{B}[/tex]

    [tex]Tr(\hat{B}) = \sum_{i=1}^{N} \langle \Psi_{i}| \hat{B} |\Psi_{i} \rangle[/tex]

    From this it follows:

    [tex] Tr(\hat{B}) = \sum_{i=1}^{N} B_{ii} [/tex]

    What's [tex]\sum_{i=1}^{N} B_{ii} [/tex] ?

    Well, it's the the sum of the diagonal elements, namely [tex] B_{11} + B_{22} + B_{33} ...[/tex]

    Therefore taking the trace means summing up the diagonal elements.
  6. Aug 21, 2005 #5
    Thanks, i fully understand now, using the derivation of Pauli spin matrices from last semester i can see how trace works, its just that our lecturer this semester has not really written anything out in matrix form, he's just mentioned it in passing.

  7. Jun 27, 2008 #6
    this explanation helped me too. but understand that the rule: tr(AB)=tr(BA) works for operators as well and I can't see why


  8. Jun 27, 2008 #7


    User Avatar
    Staff Emeritus
    Science Advisor
    Gold Member

    [tex]Tr(XY)=\sum_{a'}\langle a'|XY|a'\rangle=\sum_{a'}\sum_{a''}\langle a'|X|a''\rangle\langle a''|Y|a'\rangle =\sum_{a''}\sum_{a'}\langle a''|Y|a'\rangle\langle a'|X|a''\rangle[/tex]

    [tex]=\sum_{a''}\langle a''|YX|a''\rangle=Tr(YX)[/tex]
  9. Jun 27, 2008 #8
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook