next up previous
Up: Back to the Chemistry 3730 assignment index

Solutions to practice problem set 3

  1. Use Taylor's theorem for dealing with any small quantities:

    displaymath268

    provided x-a is small. In this case, a=0:

    displaymath274

  2. tex2html_wrap_inline276 or tex2html_wrap_inline278
    1. The field splits the states according to the value of tex2html_wrap_inline280 . There are four different values of tex2html_wrap_inline280 ( tex2html_wrap_inline284 ) so there are four states with different energies in a magnetic field.
    2. The energy shift caused by the field is proportional to tex2html_wrap_inline280 . Thus the difference in energy between two states depends only on tex2html_wrap_inline288 (and not on the actual values of the spins) so only one line appears.
    1. Instead of electrons pairing up into orbitals, up to four electrons could share an orbital because particles with a spin of 3/2 can have four different values of tex2html_wrap_inline280 . Thus, there would be 16 main groups.
    2. Particles with integer spin are bosons and are not subject to an exclusion principle. Accordingly, all of the electrons could sit in the same orbital, regardless of how many of them there were. There would be no periodicity as we know it: Every electron would be a valence electron or, to state the matter another way, each element would be in its own ``group''.
  3. If tex2html_wrap_inline294 is a solution of Schrödinger's equation, tex2html_wrap_inline296 . Then we have

    displaymath298

  4. The trial wavefunction is

    > phi := x -> x*(x-L)+c2*x^2*(x-L)^2;

    displaymath252



    For the particle in a box, tex2html_wrap_inline300 . Therefore tex2html_wrap_inline302 is simply


    > pKp:=-hbar^2/(2*m)*int(phi(x)*diff(phi(x),x$2),x=0..L);

    displaymath253



    tex2html_wrap_inline304 is

    > ip := int(phi(x)^2,x=0..L);

    displaymath254

    > Evar := pKp/ip;

    displaymath255



    Once we have the variational energy, we minimize it to find the best value of the variational parameter:

    > s1:=solve(diff(Evar,c2),c2);

    displaymath256



    There are two solutions, only one of which can be right. Recall that we are looking for the value of tex2html_wrap_inline306 which minimizes the variational energy. Substitute the two solutions back into Evar:

    > evalf(subs(c2=s1[1],Evar));

    displaymath257

    > evalf(subs(c2=s1[2],Evar));

    displaymath258

    > c2:=simplify(s1[2]);

    displaymath259

    > hbar:=h/(2*Pi);

    displaymath308

    > simplify(Evar);

    displaymath310

    > evalf(");

    displaymath312



    The exact value of the ground-state energy is tex2html_wrap_inline314 or tex2html_wrap_inline316 . This is very close to the value obtained with our simple variational wavefunction but, in accordance with the variational theorem, slightly lower.

    To normalize the wavefunction, replace tex2html_wrap_inline318 by tex2html_wrap_inline320 , where A is a normalization factor. We want tex2html_wrap_inline324 . This gives tex2html_wrap_inline326 , but we have already computed this inner product:


    > A := simplify(1/sqrt(ip));

    displaymath328




next up previous
Up: Back to the Chemistry 3730 assignment index

Marc Roussel
Wed Nov 5 14:30:22 MST 1997