qutip

QuTiP - Quantum Dynamics

QuTiP is designed to be a flexible and efficient framework for quantum mechanics. It allows for easy creation of quantum objects and provides powerful solvers to track their evolution in both closed and open environments.

When to Use

• Simulating the time evolution of a quantum system (Schrödinger or Master equation)
• Calculating steady states of open quantum systems (e.g., a cavity under drive and dissipation)
• Analyzing entanglement, Wigner functions, and quantum correlations
• Studying light-matter interactions (Jaynes-Cummings model, Rabi oscillations)
• Pulse sequence optimization and quantum control
• Calculating spectrums and multi-time correlation functions
• Parallelizing quantum simulations across multiple CPUs

Reference Documentation

Official docs: https://qutip.org/docs/latest/
Tutorials: https://qutip.org/tutorials.html
Search patterns: qutip.Qobj, qutip.mesolve, qutip.wigner, qutip.expect, qutip.basis

Core Principles

The Qobj (Quantum Object)

The central data structure. A Qobj can represent a state (ket or bra), an operator, or a superoperator. It automatically handles dimensions and validates operations (e.g., preventing the addition of a ket and an operator).

Composite Systems

QuTiP uses the Kronecker product (tensor) to represent systems composed of multiple sub-systems (e.g., a qubit coupled to a resonator).

Solver Workflow

1. Define the Hamiltonian (H)
2. Define collapse operators (C_n) for dissipation
3. Set initial state (ρ₀)
4. Define time sequence (t)
5. Run the solver (mesolve, sesolve, etc.)

Quick Reference

Installation

pip install qutip

Standard Imports

import qutip as qt
import numpy as np
import matplotlib.pyplot as plt

Basic Pattern - Rabi Oscillations

import qutip as qt
import numpy as np

# 1. Define Hamiltonian: H = h_bar * omega * sigma_x / 2
omega = 1.0 * 2 * np.pi
H = 0.5 * omega * qt.sigmax()

# 2. Initial state: Spin down (|1>)
psi0 = qt.basis(2, 1)

# 3. Time points
times = np.linspace(0.0, 10.0, 100)

# 4. Solve Schrödinger Equation
result = qt.sesolve(H, psi0, times, [qt.sigmaz()])

# 5. Get expectation values
sz_expt = result.expect[0]