to index              2015_10_15 Previous lecture note             LECTURE NOTE    2015_11_18           2015_11_19 Next lecture note


Light -Atom Interaction - Part 1: Quantum Time Developement


Most focus in Quantum Mechanics courses is on the origin of discrete energy levels - and thus the stationary states.
Schrödinger equation is introduced usually in time dependent form, even with possibility of time-dependent total
energy (hamiltonian operator explicitely time dependent)
but then all this is quickly limited to time-independent hamiltonians which can then lead to stationary states.
Mathematically - this is shown by the possibility of sepeartion of time and space coordinates.

Next step is then often the postulate that "only eigenstates are possible outcomes of measurement", and only in
more advanced or complete courses we see some treatment of
         -   time-dependent perturbation theory
         -   Fermi Golden Rule  -  for decay rate of excited systems
         -   Ramsey or Rabi oscillations

In this part of our course we thus first visit some simple systems, look at the time-dependent Schrödinger equation,
derive "Fermi Golden Rule"  -  for decay rate of excited systems, but concentrate specially on the natural line widths.
(this and next lecture)

In many texts (even in scientific literature) the existence of natural line widths is ascribed to
"time-energy Heisenberg uncertainety relation",
and the Dirac delta function in the "Fermi Golden Rule" as a "proof" of energy conservation.

We shall see how seriously you should take such statements and those textbooks (hint: not much).

We shall start by time developement in a two-state problem - (well known from some courses)
TWO identical potential wells              "competing" to keep bound one particle                (think it to be an electron)
                                   
                                            Instead of the pictures hamiltonian -  for 2 levels   H = H0 + V1  + V2      

Then we continue by a simple treatment of  a general time-dependent total hamiltonian
in a form often used to introduce time-dependent perturbation theory.




WARNING     since the captured text is from different parts   - the notation is not always consistent

          one of the   confused  notations:    sometimes    H = H0 + H'       and sometimes     H = H0 + H1
                                                                                   thus       H'        and    H1    denote both   "the extra"  coupling energy term   
         another problem is the  notation for the "initial state", sometimes index   a,    sometimes   0   or even  

         In all cases the meaning should be clear from the context


Two identical potential wells - for simplicity assume that only the chosen potential can "keep" only one bound state.
When these are brought together, i.e. for example one electron bound in two closely placed attracting potential wells,
this system gets two close lying energy eigenstates.

If we drop electron (prepare the system) in one well, we shall get probability density oscillating with time

      z_0010_two_wells_system_vs_2_single_wells.png

       z_0010_two_wells_system_vs_2_single_wells.png


Well known result illustrated by animation



This animation shows
time developement of
PROBABILITY DENSITY

We know that eigenstates of the whole                       H = H0 + V1  + V2           form a sort of split doublet (if the wells are moderately separated)
More details are below
Physics:   the "electron" is placed in a state at                   t=0                  confined to    ONE WELL ONLY         (1 well eigenstate)

The probability   is seen to be    OSCILLATING    between the two wells    
                                                                               This is seen from the EIGENSTATES of the whole problem
Note that we discuss two "sets" of state, the phi - eigenstates of simplified system, and the psi eigenstates of the whole
system

      z_0020_two_well_eigenfunctions_vs_2_single_wells.png

       z_0020_two_well_eigenfunctions_vs_2_single_wells.png


Here we discuss in several rounds the following:
We can expand any state i the particular domain of wavefunctions according to many different basis sets

For us is this of interest:
                             1.  Eigenstates of a simplified system
                              2. Eigenstates of a "full" system          (see the formal definitions below)   here
                                                                                                                                or with an illustration here

The matrix S we mention in the next plate is actually the matrix of eigenstates - see the matlab demo below

      z_0023_NOTES-on-two_well_eigenfunctions_vs_2_single_wells.png

       z_0023_NOTES-on-two_well_eigenfunctions_vs_2_single_wells.png


The mathematics and an explanation of the animation

      z_0030_two_wells_time_developement_superposition.png

       z_0030_two_wells_time_developement_superposition.png


In next two plates we illustrate the psi and phi, and their mutual transformations;
Then we show how we can obtain the      time developement in two ways   
        1. solve the Schrödinger equation - time dependent using expansion in the phi-set
        2. use the fact that the solutions of the "total problem" have "independent" time-dependent solution (the exponential energy-phase)
                       and expand each of the psi   in terms of phi   -  where one of the phi-set is the initial condition
                                                                                                   (electron in the left well; see next plate; here general formulae)

      z_0050_eigenstates_of_H_vs_H_0__expansions.png

       z_0050_eigenstates_of_H_vs_H_0__expansions.png


The yellow-plate below is a copy of the above one - but with added psi and phi of the TWO WELL PROBLEM (white spots)

      z_0055_eigenstates_of_H_vs_H_0__expansions_TWO-WELL.png

       z_0055_eigenstates_of_H_vs_H_0__expansions_TWO-WELL.png


Now we turn to different two wells:                                     (upper part:  already discussed oscillating probability)

     - in the lower part:    shown are the many eigenstates    -  quasicontinuum
     - and the result - for really many states, quasicontinuum, we will observe a decay of the original density

      z_0100_2_wells_with_quasi-continuum-decay_not_oscillation.png

       z_0100_2_wells_with_quasi-continuum-decay_not_oscillation.png


                     Animations; the two well illustration
                                                              
              two equal wells    oscillations                                                     two wells  -  one well with quasicontinuum - decaying population

              Below:    The same slide as the above  - but this is the original  - with notes   from the lecture        

      z_0105_2_wells_with_quasi-continuum-decay_not_oscillation.png

       z_0105_2_wells_with_quasi-continuum-decay_not_oscillation.png


Two wells:                       one with many eigenstates    -  quasicontinuum

          and the result - for really many states, quasicontinuum, we will observe a decay of the initial probability density
          also
         another illustration of the important aspect of "two sets of eigenstates" - those of  H0 and those of H = H0 + H'  (or H = H0 + H1  ) 

      z_0110_2_wells_with_quasi-continuum-qutoioniz_example.png

       z_0110_2_wells_with_quasi-continuum-qutoioniz_example.png


Reminder - how do find the eigenstates; expansion in eigenstates of "simpler" system

      z_0120_finding_eigenstates_expansion.png

       z_0120_finding_eigenstates_expansion.png


               Two wells:                       one with many eigenstates    -  quasicontinuum
               and
               Energies in a matrix version of the model

      z_0150_finding_eigenstates_expansion.png

       z_0150_finding_eigenstates_expansion.png

 
  Snapshots from folk.uib.no/nfylk/PHYSTOYS/golden/
  This webpage is a javascript program solving the TDSE (each new parameter set takes a few secs)
                                              (   try the N=2, N=4, .... N=20 clicks )

      z_0200_solving_time-dep-Schrodinger.png

       z_0200_solving_time-dep-Schrodinger.png


The Fermi golden rule simulator illustrated

      z_0220_Golden_Rule_Simulator_solving_time-dep-Schrodinger.png

       z_0220_Golden_Rule_Simulator_solving_time-dep-Schrodinger.png


    Finding eigenvalues in matlab   - a very small matrix of the above H  ( N=8) type

      z_0230_Matlab_H_Eigenvalues.png

       z_0230_Matlab_H_Eigenvalues.png


The above examples should have shown the situation when the system is "forced" to start from a state which is an eigenstate
of a "simplified system (e.g. only one well - or better two completely isolated levels without penetration possibility)

            BUT NOT an EIGENSTATE of the total system - Then we can have a decay

In the next part we recorded the solution of time dependent Schrödinger equation using an APPROXIMATION

                                                                                                Time dependent perturbation theory

      z_0500_Time_dep_Schroedinger.png

       z_0500_Time_dep_Schroedinger.png


Time - dependent Schrödinger equation    - expansion in phi - before applying    Time dependent perturbation theory

      z_0510_Time_dep_Schroedinger-MATRIX_FORM.png

       z_0510_Time_dep_Schroedinger-MATRIX_FORM.png

 
   Time dependent perturbation theory

      z_0520_Time_dep_Perturbation_Theory.png

       z_0520_Time_dep_Perturbation_Theory.png


    Time dependent perturbation theory    
                         assumptions:   The change is small, the original expansion coefficient remains close to one ( = 1 ) 

      z_0530_Perturbation_Theory_assumptions.png

       z_0530_Perturbation_Theory_assumptions.png

 
   Time dependent perturbation theory           Eliminating the time exponential phase factor
                                                                           only "the coupling" hamiltonian remains

      z_0540_Matrix_to_Column_decouple_equations.png

       z_0540_Matrix_to_Column_decouple_equations.png

 
   Time dependent perturbation theory            Decoupling of the coupled equations - see above

      z_0550_Perturbation_decouple_equations.png

       z_0550_Perturbation_decouple_equations.png

 
   Time dependent perturbation theory           the decoupled equations are easily solved
   
      z_0560_Evaluate_time_integrals_perturbation.png

       z_0560_Evaluate_time_integrals_perturbation.png

 
   Time dependent perturbation theory  - rearranging the "probability expression" - as function of omega times time

      z_0580_Evaluate_time_integrals.png

       z_0580_Evaluate_time_integrals.png

 
   Time dependent perturbation theory              the rearranged "probability"  will look like a delta function of omega for large t
                                                                               see below

      z_0600_Dirac_delta_origin.png

       z_0600_Dirac_delta_origin.png


   Time dependent perturbation theory              the rearranged "probability"  as a  function of  x = omega t

      z_0610_Dirac_delta_Sine_step.png

       z_0610_Dirac_delta_Sine_step.png


   Time dependent perturbation theory              the rearranged "probability"  looks like a delta function of omega for large t

      z_0620_Dirac_delta_Approaching.png

       z_0620_Dirac_delta_Approaching.png

 
   Time dependent perturbation theory           summing the probabilities over "all the relevant final states - integral in continuum limit
   
      z_0650_Probability_of_change.png

       z_0650_Probability_of_change.png

 
   Time dependent perturbation theory           summing the probabilities over "all the relevant final states - integral in continuum limit
   
      z_0680_Summing_Probability_of_change--Fermi_Golden.png

       z_0680_Summing_Probability_of_change--Fermi_Golden.png

 
   Time dependent perturbation theory           summing the probabilities over "all the relevant final states - integral in continuum limit
   
    Fermi Golden rule

      z_0700_Fermi_Golden_Rule.png

       z_0700_Fermi_Golden_Rule.png



     Fermi Golden rule:   constant rate of probability change  from Perturbation theory   (note that the delta function disappears)

                                            The density of states appeared when we replaced the summation by an integral
                                                                                       This is always so:     summation ---> integral     allways needs "density of states"
                                                                                                                           ( see also density of states in the NEXT LECTURE )
  NEXT LECTURE:     Exponential decay from constant rate. Natural line width   and further details of the electromagnetic field


to index              2015_10_15 Previous lecture note             LECTURE NOTE    2015_11_18           2015_11_19 Next lecture note