Hydrogen ground state integral - Printable Version +- Python Forum (https://python-forum.io) +-- Forum: Python Coding (https://python-forum.io/forum-7.html) +--- Forum: Data Science (https://python-forum.io/forum-44.html) +--- Thread: Hydrogen ground state integral (/thread-13073.html) |
Hydrogen ground state integral - Eduard - Sep-26-2018 The integral is supposed to return one but it returns zero. https://quantummechanics.ucsd.edu/ph130a/130_notes/node233.html The integral I'm doing is shown in the link above in the section that says "The radial wavefunctions should be normalised as below." I'm using R_10 which is the ground state. import scipy from scipy import constants import math import numpy as np from numpy.polynomial import Laguerre from scipy import special import sympy from sympy import * from sympy import var from sympy import integrals from scipy import integrate #constants a_0 = scipy.constants.value('Bohr radius') e = scipy.constants.value('elementary charge') hbar = scipy.constants.value('Planck constant over 2 pi') a = a_0 def int3(x): c1 = (2/a)**(1.5) y1 = c1*np.exp(-x/a) m = np.conjugate(y1)*y1*x**2 return m I3= integrate.quad(int3, 0, np.inf) print(I3)Thank you in advance RE: Hydrogen ground state integral - Larz60+ - Sep-26-2018 is it really zero, or .99999999999999999999999991 ? RE: Hydrogen ground state integral - Eduard - Sep-26-2018 The output should be 1 or close to one like 0.99999999... for all values of 'a' in the code. But as i make a smaller in this equal to the bohr radius a_0 it becomes zero. |