Different Energy Pulses - Ideal Qubits

[EdLatham]

Dear Spinnach Experts,

The purpose is to look at ideal qubit systems.
I am considering an experiment demonstrating the loss in the coherence of a population over successive, different energy, EPR pulses.
i.e for a simple Hamiltonian = SDS + g \mu B S, for a S = 3/2 system where energy levels are ordered -3/2, -1/2, +1/2, +3/2, initialised in m_s = -3/2 through a very cold Maxwell-Boltzmann distribution ~ 90% population in the ground state.

Then the experiment consisting of two pulses,
-one of the corresponding energy to move the population to m_s = +3/2, E = E_{ms = +3/2} - E_{ms = -3/2}
-then another to split the population from there to m_s = +1/2, E = E_{ms = +3/2} - E_{ms = +1/2}

I would like the output to measure the intensity of that superposition of states. (ms = +3/2 & ms = +1/2)
I'm not sure if Spinnach has a density matrix plot that would be ideal or possible to implement which may be of use to display their population against the other states.

From here I could see how to change the script to initialise in other states and do similar sequences finishing in different superposed states and introduce more uncertainty and constraints created by real EPR machines. such as noise terms in the Hamiltonian, or a higher temperature initialisation.

A starting script or a point in the direction towards one that looks most similar would be very useful.
Thanks in advance,
Ed.

[kuprov]

I will not assist with a project that deals with "spin quantum computing".