Semi-classical systems

Semi-classical systems

The framework offers an implementation of semi-classical systems and their dynamics. By semi-classical systems we mean systems that are described by a quantum state (a ket or density matrix) and additionally a set of classical variables. Such a state can be implemented using semiclassical.State.

For example, under certain conditions the motion of a two-level atom can be described classically, while its internal degrees of freedom are still of quantum nature. Say the atom can move in one dimension, e.g. $x$. The atom's state would then be fully described by a ket $|\psi\rangle$ that gives the internal state of the atom (meaning whether it is excited or not) and the classical variables for the motion $x$ and $p$, where $x$ is the position in space and $p$ is the momentum.

Using the framework's implementation, we can define this state by:

b_atom = SpinBasis(1//2)

ψ = (spinup(b_atom) + spindown(b_atom))/sqrt(2)
ρ = dm(ψ)

x = 0.2
p = 7.5

u = ComplexF64[x, p] # Complex vector containing classical variables

ψ_sc = semiclassical.State(ψ, u)
ρ_sc = semiclassical.State(ρ, u)

Note, that the vector of classical variables must be complex.

Furthermore, we can calculate the dynamic, semi-classical time evolution. This can be very practical when, for example, the system Hamiltonian depends on classical variables which themselves have to be calculated dynamically. Much like the time-dependent time evolution solvers timeevolution.schroedinger_dynamic and timeevolution.master_dynamic, the solver requires functional based parameter input. However, instead of just allowing variables to depend on time, these functions can be much more general: the Hamiltonian and jump operators can depend on variables that require simultaneous solution of their own dynamical equations which can in turn depend on the system density matrix. The time evolution can then be obtained by using semiclassical.schroedinger_dynamic or semiclassical.master_dynamic.

The functions, must have the form

# fquantum for a Schrödinger time evolution
function fquantum_schroedinger(t, ψ, u)
  # update H (Hamiltonian) according to dependencies on
  # classical variables u and quantum state ψ
  return H

# fquantum for a Master time evolution
function fquantum_master(t, ψ, u)
  # update H (Hamiltonian), J and Jdagger (jump operators)
  # according to dependencies on classical variables
  # u and quantum state ψ
  return H, J, Jdagger

# fclassical
function fclassical(t, ψ, u, du)
  # update du (vector of derivatives of classical variables)
  # according to dependencies on classical variables
  # u and quantum state ψ
  # -- no return statement!

Note, that the way these functions are defined has a great impact on the performance of the calculation since they are evaluated at every time step. To save time try to optimize them, e.g. by calculating required operator products outside of the functions.

We can then calculate the time evolution by passing the functions to the solver like so

tout, ψt = semiclassical.schroedinger_dynamic(tspan, ψ_sc, fquantum_schroedinger, fclassical)
tout, ρt = semiclassical.master_dynamic(tspan, ψ_sc, fquantum_master, fclassical)