Modele pequenas oscilações em uma molécula de CO
Experimentalmente, uma molécula de CO oscila em torno do seu comprimento de equilíbrio com uma constante efetiva de elasticidade de . As oscilações são governados pela equação de oscilador harmônico quântico. A seguir,
é a massa reduzida de uma molécula,
é a frequência natural,
é o deslocamento da posição de equilíbrio, e
é a constante de Planck reduzida.

qho = -(\[HBar]^2/(2 m)) Laplacian[u[x], {x}] + (m \[Omega]^2)/
2 x^2 u[x];
Calcule os primeiros quatro valores e funções próprios normalizados.

sol = DEigensystem[qho, u[x], {x, -\[Infinity], \[Infinity]}, 4,
Assumptions -> \[HBar] > 0 && m > 0 && \[Omega] > 0,
Method -> "Normalize"]

Supondo que a partícula está em uma superposição igual dos quatro estados, a função de onda terá a forma .

\[Psi][x_, t_] = Total[MapThread[1/2 Exp[I E t #1/\[HBar]] #2 &, sol]]

Calcule os três parâmetros ,
, e
usando unidades de base de unidades de massa atômica, femtossegundos, e picometros, enquanto os valores resultantes estarão perto de unidade de ordem.

m = QuantityMagnitude[(
Entity["Element", "Carbon"][
EntityProperty["Element", "AtomicMass"]] Entity["Element",
"Oxygen"][EntityProperty["Element", "AtomicMass"]])/(
Entity["Element", "Carbon"][
EntityProperty["Element", "AtomicMass"]] +
Entity["Element", "Oxygen"][
EntityProperty["Element", "AtomicMass"]]), "AtomicMassUnits"]


\[Omega] =
Sqrt[QuantityMagnitude[Quantity[1.86, "Kilonewtons"/"Meters"],
"AtomicMassUnit"/"Femtoseconds"^2]/m]


\[HBar] =
QuantityMagnitude[Quantity[1., "ReducedPlanckConstant"],
"AtomicMassUnit"*"Picometers"^2/"Femtoseconds"]

A função densidade de probabilidade do deslocamento é dada por .

\[Rho][x_, t_] =
FullSimplify[ComplexExpand[Conjugate[\[Psi][x, t]] \[Psi][x, t]]]

Como uma distribuição de probabilidade, a integral de sobre os reais é 1 para todo
.

Chop[Integrate[\[Rho][x, t], {x, -\[Infinity], \[Infinity]}]]

Visualize a densidade de probabilidade ao longo do tempo.
