jrjohansson / qutip-lectures Goto Github PK

Hi.

I'm very new to Qutip and still in the process of refreshing my quantum physics knowledge so please forgive me if this is a silly question.

In the notebook on the Jaynes Cumming model I cannot understand why the hamiltonian is written with the sm.dag()*sm for the atomic part in the following:

H = wc * a.dag() * a + wa * sm.dag() * sm + g * (a.dag() * sm + a * sm.dag())

When I look at the description of Jaynes Cumming model that you give and also on Wikipedia it seems like we should be using something like sigmaz instead, i.e.

H = wc * a.dag() * a + wa *tensor(qeye(N),sigmaz())+ g * (a.dag() * sm + a * sm.dag())

I wondered whether they might be the same, but they do not appear to be when I calculate them.

Thanks for your time.

Matt

When trying to evolve the code:
mesolve(h_t, psi0, taulist, [], process_rho, args)

I get the following error. It looks like it is expecting parameters to be passed in, but I'm not sure which ones to use.
`TypeError Traceback (most recent call last)
in ()
1 # Evolve the system, request the solver to call process_rho at each time step.
2
----> 3 mesolve(h_t, psi0, taulist, [], process_rho, args)

/home/jhellmers/anaconda3/lib/python3.6/site-packages/qutip/mesolve.py in mesolve(H, rho0, tlist, c_ops, e_ops, args, options, progress_bar, _safe_mode)
224
225 if _safe_mode:
--> 226 _solver_safety_check(H, rho0, c_ops, e_ops, args)
227
228 if progress_bar is None:

/home/jhellmers/anaconda3/lib/python3.6/site-packages/qutip/solver.py in _solver_safety_check(H, state, c_ops, e_ops, args)
834 _structure_check(Hdims, Htype, _temp_state)
835
--> 836 for ii in range(len(e_ops)):
837 if isinstance(e_ops[ii], Qobj):
838 _temp_state = e_ops[ii]

TypeError: object of type 'function' has no len()`

With Python 3.8.8, the notebook throws the following:

DeprecationWarning: scipy.rand is deprecated and will be removed in SciPy 2.0.0, use numpy.random.rand instead

DeprecationWarning: scipy.real is deprecated and will be removed in SciPy 2.0.0, use numpy.real instead

DeprecationWarning: scipy.array is deprecated and will be removed in SciPy 2.0.0, use numpy.array instead