Stack Overflow
Stack Overflow

Questions tagged [qutip]

QuTiP stands for Quantum Toolbox in Python. It basically aims to be one-stop shop to carry out any quantum mechanics based calculations in python.

From qutip.org

QuTiP is open-source software for simulating the dynamics of open quantum systems. The QuTiP library depends on the excellent Numpy, Scipy, and Cython numerical packages. In addition, graphical output is provided by Matplotlib. QuTiP aims to provide user-friendly and efficient numerical simulations of a wide variety of Hamiltonians, including those with arbitrary time-dependence, commonly found in a wide range of physics applications such as quantum optics, trapped ions, superconducting circuits, and quantum nanomechanical resonators. QuTiP is freely available for use and/or modification on all major platforms such as Linux, Mac OSX, and Windows*. Being free of any licensing fees, QuTiP is ideal for exploring quantum mechanics and dynamics in the classroom.

63 questions
15
votes
1 answer

Open quantum system modelling

I have been working for a long time now on modelling an open quantum system using the Lindblad Equation. The Hamiltonian is the following: However, two other matrices are added to the Hamiltonian. One of them has all the diagonal terms equal to…
TanMath
  • 598
  • 1
  • 9
  • 27
5
votes
2 answers

python animation.FuncAnimation error : object is not iterable

I'm new to Python and now studying matplotlib to use animation function. I'm following qutip tutorial because of my study. But when I copied and pasted example code of qutip tutorial, it didn't work Error message was Axes3D object is not…
park yunseok
  • 53
  • 1
  • 1
  • 4
4
votes
2 answers

Fast eigenvectors finding on a HPC using Qutip and slepc4py

I am doing some numeric simulations of quantum computation, and I wish to find the eigenvectors of a big hermitian matrix (~2^14 rows/columns) I am running on a 24 core/48 threads XEON machine. The code was originally written with the help of the…
oyon
  • 81
  • 3
3
votes
1 answer

How do you remove the exponential text in axis values in matplotlib?

On the x-axis, my units are in ns as this is what the physics dictates. However, I don't want the exponential text highlighted to be seen. Is there anyway this can be removed while keeping the x-axis units the same? One further question is with…
user6061964
3
votes
1 answer

Python QuTip: Plot full lines in Bloch sphere

I've recently started using QuTip for Python. I've want to plot the evolution of a two-level system in a Bloch sphere using qutip.Bloch(), but I don't know how to make such a plot as a continuous line. This is the code I currently use: bola =…
Pablo Poggi
  • 31
  • 2
2
votes
0 answers

Using qutip.mesolve() to calculate the time evolved Hamiltonian with time dependent variable

I'm trying to find the time evolution of some initial state with a time dependent hailtonian where also the z variable (gradient of time evolution) is time dependent. Usually the qutip.mesolve() function can be used for this when the z variable is a…
Calvin
  • 21
  • 1
2
votes
1 answer

How to repeat a function for a list of inputs?

I'm trying to use the following function of QuTip library: coherent_dm(N=10, alpha = c, offset=0, method='operator') It's output is a matrix and it's input is a complex number "c", however I want to get the matrices for a list of multiple "c"…
Bruno Piveta
  • 25
  • 3
2
votes
1 answer

Qutip Raise Warning When Specify Bra Directly From Numpy Array Transpose

I try to define a bra of a ket using qutip. A bra is a row vector and a ket is a column vector. A bra can be defined as a ket conjugate transpose. However, if I specify a bra directly using this definition in qutip a warning always raised, and the…
Willy satrio nugroho
  • 908
  • 1
  • 16
  • 27
2
votes
1 answer

When plotting the Wigner function of a coherent state using QuTiP strange patterns appear

I noticed something strange this day when I plotted the Wigner function of a coherent state using the open source quantum toolbox QuTiP in python. When I do the plot I noticed these strange patterns just around the edge of the plot that are not…
Turbotanten
  • 386
  • 3
  • 14
2
votes
1 answer

Live plotting on bloch sphere

I am trying to plot live data on a bloch sphere using Qutip's function bloch(). So far, the code always interrupts, when I have a b.show() in there. I found a lot of solutions online to similar problems, but most of them make use of direct…
Mechanix
  • 95
  • 1
  • 8
2
votes
0 answers

Why is this qutip code not working?

I have been coding to calculate the following equation: Where rho(t) is calculated by solving a quantum mechanics equation. <3| is a row vector with 1 at the 3rd element. |3> is a column vector with 1 at the 3rd element. To solve the quantum…
TanMath
  • 598
  • 1
  • 9
  • 27
2
votes
1 answer

Integration not successful in Python QuTiP

I have been trying to use QuTiP to solve a quantum mechanics matrix differential equation (a Lindblad equation). Here is the code: from qutip import * from matplotlib import * import numpy as np hamiltonian = np.array([[215, -104.1, 5.1, -4.3 …
TanMath
  • 598
  • 1
  • 9
  • 27
1
vote
0 answers

how to use qutip in visual studio code?

I am a beginner of QuTip (a library of python?) user. I have already installed my qutip in my environment by conda. I can execute simple program in jupyter book but not in VScode. How can I do it in VScode? What I tried: import…
Hsu Bill
  • 111
  • 2
1
vote
1 answer

Calculate an Qobjevo operator of qutip

I define a time-dependent operator in the form of an Qobjevo and pass it to sesolve() to solve for Schrodinger's equation. Below is an MWE. from qutip import * import numpy as np tlist = np.linspace(0, np.pi / 2) H = [ sigmaz(), [sigmax(),…
Neo
  • 1,031
  • 2
  • 11
  • 27
1
vote
1 answer

Calculate the time dependence expectation value of a quantum operator with tensorflow

I'm trying to calculate the expectation value of a quantum operator in function of time. I obtained yes the states for every istance t_l = tf.constant( np.linspace(0, 4 / gam_q, 1000, dtype = np.complex128) ) dt = t_l[1] - t_l[0] H_t =…
Alberto Mercurio
  • 21
  • 2
1
2 3 4 5