Hamiltoniana cuántica 

En mecánica cuántica estándar, los sistemas evolucionan de acuerdo con la ecuación de Schrödinger , donde es una matriz Hermitaniana llamada Hamiltoniana. Los siguientes son posibles Hamiltonianas.

In[1]:=

X \[ScriptCapitalH]1 = {{0, -I}, {I, 0}}; \[ScriptCapitalH]2 = {{1, 0, 1 + I}, {0, 1, 0}, {1 - I, 0, 2}};

In[2]:=

X {HermitianMatrixQ[\[ScriptCapitalH]1], HermitianMatrixQ[\[ScriptCapitalH]2]}

Out[2]=

Las siguientes matrices no pueden ser Hamiltonianas cuánticas porque no son Hermitianas.

In[3]:=

X bad1 = {{0, 1}, {-1, 0}}; bad2 = {{1, 0, I}, {0, 1 - I, 2}, {-I, 1 + I, 5}};

In[4]:=

X {HermitianMatrixQ[bad1], HermitianMatrixQ[bad2]}

Out[4]=

La matriz exponencial es llamada el operador de evolución de tiempo y es siempre una matriz unitaria (asumiendo que el tiempo y que la constante de Planck ℏ son reales).

In[5]:=

X u1 = MatrixExp[-((I \[ScriptCapitalH]1)/\[HBar]) t]

Out[5]=

In[6]:=

X u2 = MatrixExp[-((I \[ScriptCapitalH]2)/\[HBar]) t]

Out[6]=

In[7]:=

X Function[m, UnitaryMatrixQ[m, SameTest -> (Simplify[#1 == #2, {t, \[HBar]} \[Element] Reals] &)]] /@ {u1, u2}

Out[7]=