Code covered by the BSD License  

Highlights from
Hydrogenic Wavefunction Visulization Tool

5.0

5.0 | 1 rating Rate this file 15 Downloads (last 30 days) File Size: 56.59 KB File ID: #26218
image thumbnail

Hydrogenic Wavefunction Visulization Tool

by Jos Martin

 

29 Dec 2009 (Updated 04 Jan 2010)

Visualize hydrogenic wave-functions in 3D and watch animations of atomic state changes

| Watch this File

File Information
Description

The code provides base functionality to calculate the wave-functions for hydrogen and visualise those wave-functions using iso-surfaces in 3D. This can be used to explore the spatial nature of the wave-functions.

In addition it is easy to calculate something approximating the wave-function of the atom whist undergoing a state transition (to make it easier to understand we assume the Rabi frequency is very slow). Then the wave-function of the atom can be visualised as a function of time whilst a transition is simulated. The transition GUI allows the user to program in any desired transition (even non-electric dipole transitions) and watch the evolution of the wave-function.

The animations clearly show why pi-transition emit linearly polarized light and sigma-transitions circular polarized light.

Developed in conjunction with Oxford Physics. Many thanks to Dr. David Lucas and Mr. Chris Ballance. The code in this submission was written by Chris Ballance whilst undertaking an internship at The MathWorks.

MATLAB release MATLAB 7.8 (R2009a)
Tags for This File  
Everyone's Tags
Tags I've Applied
Add New Tags Please login to tag files.
Comments and Ratings (3)
14 Aug 2011 Oren

Hi, thanks for this contribution.
I am also using Hydrogen eigenfunction \psi=Ylm(\theta,\phi)*Rnl(r) and encountered a problem which I detect in this tool as well, regarding the file: "calculateHydrogenicWavefunction.m". It appears that your wavefunctions are not orthonormal... for example the inner-product of the wavefunctions for (N,L,M)=(2,1,0) and (3,1,0) does not vanish.
If psi1 is the matrix for (2,1,0) and psi2 is for (3,1,0) then sum(sum(sum(psi1*conj(psi2)))) is -2.7754. I suspect it has something to do with the implementation of the radial functions.
Is this a numerical by product? How can this be avoided?

Thanks a lot!
Oren

14 Aug 2011 Oren

Hi, looking at your code again I noticed that in calculating psiNorm in line 168 in: "calculateHydrogenicWavefunction.m" you are missing a power of 3 over factorial(N+L), that is according to :introductory quantum mechanics" by Liboff, table 10.3

08 Oct 2011 Mr Smart  
Please login to add a comment or rating.
Updates
04 Jan 2010

Copyright notices added.

04 Jan 2010

Copyright updated.

Tag Activity for this File
Tag Applied By Date/Time
physics simulation Jos Martin 04 Jan 2010 10:17:01
physics Jos Martin 04 Jan 2010 10:17:01
hydrogen Jos Martin 04 Jan 2010 10:17:01
animation Jos Martin 04 Jan 2010 10:17:01
wavefunction Jos Martin 04 Jan 2010 10:17:01
hydrogen anas daniel 22 May 2011 04:08:31

Contact us at files@mathworks.com