PHYS 860:  Atomic and Molecular Dynamics
Spring 2007
TU 11:30, CW 143


    Discuss them with your classmates,
    Email me at, or
    Stop by my office, CW 329.

Guidelines for homework:

  I encourage you to discuss the homework assignments with your classmates, but write your homework solutions on your own!

  As a scientist in training, you need to learn to communicate scientific information in an effective, efficient manner.  You should consider homework assignments as practice in this art.  It is your responsibility to present your homework solutions in a readable, logical manner not mine to decipher and interpret them.

  In an incorrect solution, I will reward statements in homework and on exams that show you know it is incorrect and why.  I will also reward any effort above and beyond what is explicitly asked for in a problem.

bulletK-matrix for HW #2 as a function of total collision energy for coupling strength 0.005

bulletK-matrix for HW #2 as a function of total collision energy for coupling strength 0.001

Computer codes:
bulletSimple fortran code for solving single-channel scattering problems with a partial wave decomposition;
bulletCalculates phase shifts and cross sections;
bulletUses the method discussed in class;
bulletIincludes Bessel function subroutines that might be useful in your own codes.
bulletSimple fortran code for finding bound states of a one-dimensional potential;
bulletUses finite differencing on a uniform grid for the kinetic energy;
bulletAssumes wave function at each boundary is zero.
bulletInput file for FiniteDiffs.f that is hopefully self-explanatory.
bulletUse it under Unix/linux as:
        > FiniteDiffs.x < FiniteDiffs.inp
assuming FiniteDiffs.x is the executable.
bulletFortran code for solving multichannel Schrodinger equation;
bulletCalculates NumOpenChannels (specified in the input file) linearly independent solutions;
bulletTo turn this into a scattering code, you would have to add the matching step to extract the S- or K-matrix;
bulletTo change the coupling strength, edit V(1,2) and V(2,1) in the subroutine CalcHamiltonian
bulletUse it under Unix/linux as:
        > MultichannelDiffs.x < MultichannelDiffs.inp
assuming MultichannelDiffs.x is the executable.
bulletInput file for MultichannelDiffs.f;
bulletExplanation included in the file.

   Due April 5: Mid-term
   Problem 1: Cross section, K-matrix, S-matrix

Homework assignments:
    Due Jan. 25:  HW 1
    Due Feb. 15:  HW 2