How to plot the different vibrational wave functions for the ground state of 1D harmonic oscillator?
14 views (last 30 days)
%Here in this code I want to plot the different vibrational wave functions for the ground electronic state of Iodine molecule.
%I took the potential to be harmonic. I put all the parameters for the iodine molecule from literature. but, unable to plot the different wave functions.
%for every 'n' in this code, the plot is showing same as the ground vibrational state which is a gaussian. kindly help and suggest!
%The form of the wave function is attached. check the attached SS.
%%%%%%All the calculations are done in AU (atomic unit)
% define various conversion constants
dr = r(2) - r(1); % Coordinate step.
% mass of I2
% X (ground) state of I2---not used here for dynamics but useful for
V1= 1/2*mass*(we_X*(r-re_X).^2); %ground state potential
% V1= De_X*(1-exp(-be_X*(r-re_X))).^2;
% % initial ground state wavefunction (ground state of harmonic oscillator---an
% % approximation here, but a good one.)
n=2; %quantum no
norm= 1/sqrt(dot(psi_init, psi_init));