-- a number of Fermionic modes or spin-orbitals
NF=6
-- a number of Bosonic modes (phonon modes, ...)
NB=0
-- For a p-shell we would like the have 6
-- spinorbitals with the quantum numbers
-- spin down ml=-1,ml=0,ml=1 and
-- spin up with ml=-1, ml=0, ml=1
-- We can group different spin-orbitals into
-- lists and assign meaning to them
IndexDn={0,2,4}
IndexUp={1,3,5}
-- the code knows that a 3 fold degenerate shell
-- has l=1 and ml=-1, 0 and 1 are
-- assigned to them automatically
-- the wave-function with one electron in the
-- px orbital with spin down is created as
psipx = NewWavefunction(NF, NB, {{"100000",math.sqrt(1/2)}, {"000010",math.sqrt(1/2)}})