Atomic Density Matrix Tutorial | Functions »|Tutorials »|More About » |

In quantum-mechanical matrix notation, the expansion coefficients of a general ket in a particular basis are notated as a column vector, expansion coefficients of a bra are notated as a row vector, and operators are notated as a square matrix. This notation is adapted for use in the AtomicDensityMatrix package.

There is one ambiguous case: A 1×1 matrix satisfies the form of a covariant or contravariant tensor, or an operator. However, since in any case the object transforms as a scalar, this does not cause much difficulty.

Checking for various structures in matrix notation.

In[1]:= |

In[2]:= |

Out[2]= |

In[3]:= |

Out[3]= |

Here is a contravariant state vector (a 31 array). In MatrixForm it is displayed as a column vector.

In[4]:= |

Out[4]//MatrixForm= | |

In[5]:= |

Out[5]= |

In[6]:= |

Out[6]//MatrixForm= | |

In[7]:= |

Out[7]= |

In[8]:= |

Out[8]//MatrixForm= | |

In[9]:= |

Out[9]= |

Here are the spherical components of the J=1/2 representation of the angular momentum operator, which make up a covariant rank-1 tensor. MatrixForm displays it as a row vector of matrix operators.

In[10]:= |

Out[10]//MatrixForm= | |

In[11]:= |

Out[11]= |

The contravariant components are obtained by taking the Hermitian conjugate. MatrixForm displays the resulting expression as a column vector of matrix operators.

In[12]:= |

Out[12]//MatrixForm= | |

In[13]:= |

Out[13]= |

The list of Cartesian components of a vector operator is not displayed in the proper format by MatrixForm; what we want is each component displayed as a matrix. To do this we can use the TensorForm function, which is equivalent in this case to applying MatrixForm to each component. TensorForm could also have been used in the above cases, where it would have produced the same result as MatrixForm.

Here are the Cartesian components of the operator (proportional to the Pauli matrices), displayed using TensorForm.

In[14]:= |

Out[14]= |

In[15]:= |

Out[15]= |

Any operator representation on a space with angular momentum *J* can be decomposed into a sum of irreducible tensor operators. The resulting decomposition is represented as a list of spherical tensors with ranks zero through 2J.

Here is the tensor decomposition of an operator on a J=1 space (e.g., the density matrix for a J=1 state). It consists of a rank-zero, rank-one and rank-two operator.

In[16]:= |

Out[16]= |

In[17]:= |

Out[17]= |

An expression is considered to be a scalar if it does not match any of the tensor forms listed above.

In[18]:= |

Out[18]= |