About the rotating wave Approximation

To run, you need to save to your directory:
(Naturally, you should know a bit about Matlab!)
  1. rabi.m
  2. rabieq.m
  3. rotwav.m
  4. consta.m

Step by step:

  1. Make a directory where to run
  2. click on each of the above sources (rabi.m,rabieq.m,rotwav.m , consta.m) and save them with these names as 'source' in your directory
  3. Start matlab from the directory
  4. type help rabi, help rotwav ....etc
  5. type
    rabi
    to run the programme.
  6. there are dialogues, but you can answer by simple return
  7. Read the notes about the Rotating Wave Approximation
  8. Change the parameters to obtain various versions.
  9. In particular, observe what happens when the detuning gets larger!


Exact solution, using rabieq.m, answer 'e'


Rotating Wave solution, using rotwav.m, answer 'r'



Help from rabi.m : %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% rabi.m solves the diff. equations for the two state problem exact Rabi rotating wave approximation in normal mode, uses a trick to plot step by step for export of pictures, one must replot (because of the tricks) there is a series of dialogs in this version. needs rabieq.m rotwav.m and consta.m Atomic Physics, Bergen. L. Kocbach Nov. 21 1995 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Models: -------- Help from: rabieq.m % RABI EQ Equation for the Rabi Frequency demo % % A = [ -EPSI V0*cos(omega*t) % V0*cos(omega*t) EPSI ]; % % ydot = -i*A*y; % consta.m % A = [ 0.5*omega- EPSI 0.5*V0 % 0.5*V0 -0.5*omega+ EPSI ]; % % ydot = -i*A*y; % rotwav.m % ROTating WAVe for the Rabi Frequency demo % % A = [ 0 0.5*V0*exp(i*(omega-2*EPSI)*t) % 0.5*V0*exp(-i*(omega-2*EPSI)*t) 0 ]; % % ydot = -i*A*y; %