(Kinetic Energy) Operator as a Matrix

ylem
Messages
30
Reaction score
1
Hey! I am completely confused...

In quantum mechanics at University we are learning about operators. I don't understand what you actually have to do to write an operator as a matrix...

In my notes, I'm looking at the kinetic energy operator in the Hamiltonian...

I can see for something basic - like an operator that makes things the number 3 when you have an orthonormal basis set that the matrix is simply

3 0 0
0 3 0
0 0 3

But I don't know how to adapt this to something more complex...

Thanks for your help! I hope I posted this in the right place :-s

Sam x
 
Physics news on Phys.org
Sometimes, you can write an operator as a matrix. Note that, just like in ordinary linear algebra, the entries of the matrix depend on the basis you are in. You shouldn't forget, a state is a vector in a (Hilbert) space, but the vector looks different in different bases you choose (in analogy with |R^3: the vector (1, 1, 0) in Cartesian coordinates and (1, pi/2, 0) (I think) in spherical coordinates are the same vector, though the components look very differently).

So for example, if I write the state of a spin system in the basis {spin up, spin down} along the z-axis, I can write the matrix of the spin z-component as
\hat S_z = \begin{pmatrix}\hbar/2 & 0 \\ 0 & -\hbar/2 \end{pmatrix}.
But when I take another basis (for example, the spin along the x axis) the matrix of the same operator will look different.

Also, some operators are infinite. For example, the position operator in one dimension has infinitely many eigenvalues (namely, one for each position x). So if we apply the \hat x operator to a state \Psi, it will give us the position. Now we can again choose a basis for \Psi, consisting of all positions x on the real line. The "components" of \Psi are usually written \langle x | \Psi \rangle or \Psi(x) (it is a vector, but it's components are not labelled by some i = 1, 2, ..., n; but by a continuous index x. Though we're looking at it as a vector, such a thing is more commonly called a function). We can write the operator as a(n infinite) "matrix", namely the diagonal matrix with all the (uncountably many!) eigenvalues x on the diagonal. But if we choose another basis, this "matrix" will look quite different. For example, if we choose as basis the momentum p, the x-operator will look like i\hbar \frac{\partial}{\partial p} in that basis and the matrix will no longer be diagonal (in fact, it is not possible to choose any basis such that the x and p operator matrices are diagonal at the same time, which have severe consequences with physical relevance).

So I'm really tired, but I hope I didn't write too much nonsense here :smile:
 
Last edited:
Matrix Operators

OK, so let me first just state that the energy eigenvalues for a particle in an infinite potential well with width 2a is as follows:

E_n = \frac{\hbar^2 \pi^2 n^2}{8ma^2}

And we know that the associated eigenwavefunctions for the system are:

\psi_n(x) = \frac{1}{\sqrt{a}} \sin (\frac{n \pi x}{2a})

Now if we have a state that contains different amounts of these wavefunctions we can write the state as a coloumn vector of eigenwavefunctions:

\Psi_A(x) = \left(\begin{array}{cc}1\\0\\-1\end{array}\right)
\Psi_A(x) = \frac{1}{\sqrt{2}}\left(\frac{1}{\sqrt{a}} \sin (\frac{\pi x}{2a}) -\frac{1}{\sqrt{a}} \sin (\frac{3 \pi x}{2a})\right)

So now recall the definition of an operator:
\hat{O}\psi = \lambda \psi

So if I use the Hamiltonian operator, I'd expect it to pull out the energy eigenvalues for the states present. So I can write this as an arbitrary dimension matrix:

\frac{\hbar^2 \pi^2}{8ma^2} \left(\begin{array}{cccc}1&0&0&0\\0&4&0&0\\0&0&9&0\\0&0&0&n^2\end{array}\right)

So that we can now write something that is easy to work with:

\hat{H}\left(\begin{array}{c}\psi_1\\\psi_2\\\psi_3\\\psi_n\end{array}\right) = \frac{\hbar^2 \pi^2}{8ma^2} \left(\begin{array}{cccc}1&0&0&0\\0&4&0&0\\0&0&9&0\\0&0&0&n^2\end{array}\right)\left(\begin{array}{c}\psi_1\\\psi_2\\\psi_3\\\psi_n\end{array}\right)

Does that help?? I'm not sure which part of the writing the operator as a matrix you didn't understand so I was very vague on details (just gave an example as an addition to CompuChips post, but worked on the other side of the equation). I think that writing matrices in calculations just simplifies equations that are in many orthogonal states. Phew, LATEX is fun :S
 
BerryBoy said:
Phew, LATEX is fun :S

It is. You can even just do something like
\begin{pmatrix}<br /> a &amp; b &amp; c \\ <br /> d &amp; e &amp; f \\<br /> \vdots &amp; \ddots &amp; \vdots \\<br /> x &amp; \cdots &amp; z<br /> \end{pmatrix}
(note that there is no array, no column specification, and no manual bracing)
Oh and \vdots and \ddots will make your matrices even more beautiful :smile:

On topic, I think you are right the TS should try to ask a more specific question now for more details.
 
Last edited:
Thanks so much for your help! I think I understand it a bit more now! :-) x
 
Hi, I had an exam and I completely messed up a problem. Especially one part which was necessary for the rest of the problem. Basically, I have a wormhole metric: $$(ds)^2 = -(dt)^2 + (dr)^2 + (r^2 + b^2)( (d\theta)^2 + sin^2 \theta (d\phi)^2 )$$ Where ##b=1## with an orbit only in the equatorial plane. We also know from the question that the orbit must satisfy this relationship: $$\varepsilon = \frac{1}{2} (\frac{dr}{d\tau})^2 + V_{eff}(r)$$ Ultimately, I was tasked to find the initial...
The value of H equals ## 10^{3}## in natural units, According to : https://en.wikipedia.org/wiki/Natural_units, ## t \sim 10^{-21} sec = 10^{21} Hz ##, and since ## \text{GeV} \sim 10^{24} \text{Hz } ##, ## GeV \sim 10^{24} \times 10^{-21} = 10^3 ## in natural units. So is this conversion correct? Also in the above formula, can I convert H to that natural units , since it’s a constant, while keeping k in Hz ?
Back
Top