logo

Matrix Representation of the Harmonic Oscillator Operator 📂Quantum Mechanics

Matrix Representation of the Harmonic Oscillator Operator

Explanation

Harmonic


🚧 이 포스트는 아직 이관 작업이 완료되지 않았습니다 🚧

The Schrödinger equation is a linear equation, so a linear combination of multiple wavefunctions that satisfy the equation also satisfies the equation. If the eigenfunctions of each state of the harmonic oscillator are ψ0\ket{\psi_{0}}, ψ1\ket{\psi_{1}}, \cdots, and ψn>|\psi_{n}>, a linear combination of these eigenfunctions, ψ=c0ψ0+c1ψ1++cnψn>\ket{\psi}=c_{0}\ket{\psi_{0}} + c_{1}\ket{\psi_{1}} + \cdots + c_{n}|\psi_{n}>, is also a solution to the Schrödinger equation. (This does not mean it is an eigenfunction. Generally, a linear combination of eigenfunctions is not an eigenfunction.) Let the basis basis\mathrm{basis} be (ψ0, ψ1,,ψn) (\psi_{0},\ \psi_{1},\cdots , \psi_{n}) and sequence them starting from the ground state. Then the eigenfunctions from the ground state can be represented by the following matrix:

ψ0=(100),ψ1=(010),ψn>=(001) \ket{\psi_{0}} = \begin{pmatrix} 1 \\ 0 \\ \vdots \\ 0 \end{pmatrix}, \ket{\psi_{1}} = \begin{pmatrix} 0 \\ 1 \\ \vdots \\ 0 \end{pmatrix}, |\psi_{n}> = \begin{pmatrix} 0 \\ 0 \\ \vdots \\ 1 \end{pmatrix}

For those unfamiliar with linear algebra, consider how each unit vector is expressed in the Cartesian coordinate system. x^=(1, 0, 0) \hat{\mathbf{x}} = (1,\ 0,\ 0), y^=(0, 1, 0),z^=(0, 0, 1)\hat{\mathbf{y}} = (0,\ 1,\ 0), \hat{\mathbf{z}} = (0,\ 0,\ 1) $

따라서 각 고유함수들의 선형결합인 ψ\ket{\psi}는 아래와 같은 행렬로 표현할 수 있다.

ψ=c0ψ0+c1ψ1++cnψn>=(c0c1cn) \ket{\psi}=c_{0}\ket{\psi_{0}} + c_{1}\ket{\psi_{1}} + \cdots + c_{n}|\psi_{n}> = \begin{pmatrix} c_{0} \\ c_{1} \\ \vdots \\ c_{n} \end{pmatrix}

사다리 연산자의 행렬 표현을 구하기에 앞서 간단한 예시로 그 원리를 설명하겠다.임의의 2×22\times 2행렬인 AA가 있다.이 행렬의 성분을 모르기에 아래와 같이 표현하자.A=(abcd)A= \begin{pmatrix} a & b \\ c & d \end{pmatrix}이 행렬의 1행 1열 성분을 뽑아내는 방법은 앞뒤로 두 행렬 (10)\begin{pmatrix} 1 & 0 \end{pmatrix}, (10)\begin{pmatrix} 1 \\ 0 \end{pmatrix}을 곱하는 것이다.(10)(abcd)(10)=(10)(ac)=a\begin{pmatrix} 1 & 0 \end{pmatrix} \begin{pmatrix} a & b \\ c & d \end{pmatrix} \begin{pmatrix} 1 \\ 0 \end{pmatrix} = \begin{pmatrix} 1 & 0 \end{pmatrix} \begin{pmatrix} a \\ c \end{pmatrix} = a 같은 방법을 써서 각 성분을 구하는 과정은 아래와 같다.

(10)(abcd)(01)=(10)(bd)=b\begin{pmatrix} 1 & 0 \end{pmatrix} \begin{pmatrix} a & b \\ c & d \end{pmatrix} \begin{pmatrix} 0 \\ 1 \end{pmatrix} = \begin{pmatrix} 1 & 0 \end{pmatrix} \begin{pmatrix} b \\ d \end{pmatrix} = b

(01)(abcd)(10)=(01)(ac)=c \begin{pmatrix} 0 & 1 \end{pmatrix} \begin{pmatrix} a & b \\ c & d \end{pmatrix} \begin{pmatrix} 1 \\ 0 \end{pmatrix} = \begin{pmatrix} 0 & 1 \end{pmatrix} \begin{pmatrix} a \\ c \end{pmatrix} = c

(01)(abcd)(01)=(01)(bd)=d \begin{pmatrix} 0 & 1 \end{pmatrix} \begin{pmatrix} a & b \\ c & d \end{pmatrix} \begin{pmatrix} 0 \\ 1 \end{pmatrix} = \begin{pmatrix} 0 & 1 \end{pmatrix} \begin{pmatrix} b \\ d \end{pmatrix} = d

이 원리를 이해했다면 이제 본격적으로 조화진동자의 사다리 연산자의 행렬표현을 구해보자.이제부터 조화 진동자의 고유 함수를 간단하게 ψn>n>\psi_{n}> \equiv |n>으로 표시하겠다.그러면 고유함수와 사다리 연산자의 관계식을 다음과 같이 쓸 수 있다.a+n>=n+1n>a_{+}|n>= \sqrt{n+1}|n>이를 이용하여 a+a_{+}의 각 성분을 구할 수 있다.조화 진동자의 상태는 n=0n=0부터 시작함에 주의하자.또한 mn=δmn\braket{m | n}=\delta_{mn}를 이용하면 n+1a+n\braket{n+1 | a_{+}|n}일 때만 행렬의 성분이 존재함을 알 수 있다.다른 성분은 전부 00이다.예를 들어 1a+1=122=212=0\braket{1 | a_{+}|1}=\braket{1 | \sqrt{2}|2}=\sqrt{2}\braket{1 | 2}=0 2a+1=222=222=2 \braket{2 | a_{+}|1}=\braket{2 | \sqrt{2}|2}=\sqrt{2}\braket{2 | 2}=\sqrt{2} 00If we only find the non-zero components1a+0=1\braket{1 | a_{+}|0}=1 2a+1=2 \braket{2 | a_{+}|1}=\sqrt{2} 3a+2=3 \braket{3 | a_{+}|2}=\sqrt{3} \vdots n+1a+n=n+1 \braket{n+1 | a_{+}|n}=\sqrt{n+1} Thus, the matrix of a+a_{+} has non-zero value n+1\sqrt{n+1} at row (n+2)(n+2) and column (n+1)(n+1) (Remember again that the state of the harmonic oscillator starts from n=0n=0).a+=(0000010000020000030000040)a_{+}=\begin{pmatrix} 0 & 0 & 0 & 0 & 0& \cdots \\ \sqrt{1} & 0 & 0 & 0 &0 & \cdots \\ 0 & \sqrt{2} &0 & 0 & 0 & \cdots \\ 0 & 0 & \sqrt{3} &0 & 0 & \cdots \\ 0 & 0& 0& \sqrt{4} &0 & \cdots \\ \vdots & \vdots & \vdots & \vdots & \vdots \end{pmatrix} In the same way, the matrix representation of aa_{-} can also be obtained. The process is omitted; try it yourself. It can be seen that the matrix representationa=(0100000200000300000400000)a_{-}=\begin{pmatrix} 0 & \sqrt{1} & 0 & 0 & 0& \cdots \\ 0 & 0 & \sqrt{2} & 0 &0 & \cdots \\ 0 & 0 & 0 & \sqrt{3} & 0 & \cdots \\ 0 & 0 & 0 & 0 & \sqrt{4} & \cdots \\ 0 & 0& 0& 0& 0 & \cdots \\ \vdots & \vdots & \vdots & \vdots & \vdots \end{pmatrix} also satisfies (a+)=a(a_{+})^{\ast}=a_{-}. The Hamiltonian of the harmonic oscillator HH can also be represented as a matrix. Knowing Hn>=(n+12)wn>H|n>=(n+\frac{1}{2})\hbar w|n>, if we use the same method as aboveH=w(120000032000005200000720000092)H=\hbar w\begin{pmatrix} \frac{1}{2} & 0 & 0 & 0 & 0& \cdots \\ 0 & \frac{3}{2} & 0 & 0 &0 & \cdots \\ 0 & 0 & \frac{5}{2} & 0 & 0 & \cdots \\ 0 & 0 & 0 & \frac{7}{2} & 0 & \cdots \\ 0 & 0& 0& 0 & \frac{9}{2} & \cdots \\ \vdots & \vdots & \vdots & \vdots & \vdots \end{pmatrix}, again, the process is omitted, so try it yourself.