A continuación se hallan varios métodos para realizar cálculos cuánticos sencillos. La idea principal de estos es calcular la energía y densidad electrónica de los fenómenos:
- Partícula en una Caja
paricula_en_caja.py - Oscilador Armónico
oscilador_armonico.py - Átomo de Hidrógeno (H)
atomo_h.py - Molécula de Hidrógeno (H2)
hf-scf.py
Cabe resaltar que lo que se halla en cada archivo solo son las funciones para realizar los cálculos; el cálculo como tal y los parámetros deseados deben ser ensamblados por el usuario. A continuación se proponen algunos ejemplos.
# Funcion para el calculo de la onda con nodos n, frecuencia f y la masa del electron
pec = psi(x = 0.5, n = 1, L = 1, m = cnts.m_e)
# Calculo de la energia de la particula en una caja
e_pec = energia(n = 1, L = 1, m = cnts.m_e)
# Grafica de la funcion de onda de la particula en una caja
orbital(n = 2, L = 1, m = cnts.m_e)
# Grafica de las funciones de onda de la particula en una caja en planos Psi[x] * Psi[y]
orbital2D(nx = 1, ny=1, Lx = 1, Ly = 1, mx = cnts.m_e, my = cnts.m_e)# Funcion para el calculo de la onda con nodos n, frecuencia f y la masa del electron
oa = psi(x = 0, n = 1, f = 0.015, m = cnts.m_e)
# Calculo de la energia del oscilador armonico
e_oa = energia(n = 1, f = 0.015)
# Grafica de la funcion de onda del oscilador armonico
orbital(n = 2, f = 0.015, m = cnts.m_e)
# Grafica de las funciones de onda del oscilador armonico en planos Psi[x] * Psi[y]
orbital2D(nx = 1, ny=1, fx = 0.015, fy = 0.015, mx = cnts.m_e, my = cnts.m_e)
# Grafica de las funciones de onda del oscilador armonico para diferentes niveles n
orbitalN(n_max = 6, f = 0.015, m = cnts.m_e):# Funcion para el calculo de la funcion de onda radial (Rho)
rho = Rho(n=1, l=0, r=0.5, Z=1)
# Funcion para el calculo de la funcion de onda azimutal (Theta)
theta = Theta(x = 0.3, y = 0.5, m = 0)
# Funcion para el calculo de la funcion de onda polar (Phi)
phi = Phi(x = 0, y = 0, z = 1, l=1, m=-1)
# Calculo de la energia del atomo de hidrogeno
e_hidrogeno = E_h(n = 1, Z = 1)
# Grafica de la funcion de onda a lo largo del radio atomico
orbital_r(n = 2, l = 1, Z = 1, d = [-0.5, 5, 0.1])
# Grafica de los orbitales del atomo de hidrogeno
orbital2D(n = 2, l = 1, m = -1)# Calculo de las constantes para la molecula de H2
orb = SCF(r = 1.4632, Z=[1,1], b1 = GTO["H"], b2 = GTO["H"], b = 2, vbs=True)
# Calculo de las energias de enlace (negativa) y antienlace (positiva)
E_anti_enlace = ener_orbs(x = orb["X"], f = orb["F"])
# Calculo de la energia electronica y total de la molecula
E_elec = ener_elec(p = orb["P"], h = orb["H"], f = orb["F"])
E_tot = ener_tot(r = 1.4632, Z=[1,1], elec=E_elec)
# Graficas de las funciones de onda a lo largo del eje
orbital(orb["C"], r = 1.4632, b1 = GTO["H"], b2 = GTO["H"], b = 2)
orbital2D(orb["C"], r = 1.4632, b1 = GTO["H"], b2 = GTO["H"], b = 2)