3 qubit resonator

In [ ]:

import numpy as np
import qutip as qt

from qubex.simulator import (
    Control,
    Coupling,
    QuantumSimulator,
    QuantumSystem,
    Qubit,
    Resonator,
)

import numpy as np import qutip as qt from qubex.simulator import ( Control, Coupling, QuantumSimulator, QuantumSystem, Qubit, Resonator, )

In [ ]:

# Create the quantum system with a qubit and a resonator (Jaynes-Cummings model)

qubit = Qubit(
    label="Q01",
    frequency=7.648,
)

resonator = Resonator(
    label="R01",
    dimension=10,
    frequency=10.123,
)

system = QuantumSystem(
    objects=[qubit, resonator],
    couplings=[
        Coupling(
            pair=(qubit, resonator),
            strength=0.01,
        ),
    ],
)

simulator = QuantumSimulator(system)

# Create the quantum system with a qubit and a resonator (Jaynes-Cummings model) qubit = Qubit( label="Q01", frequency=7.648, ) resonator = Resonator( label="R01", dimension=10, frequency=10.123, ) system = QuantumSystem( objects=[qubit, resonator], couplings=[ Coupling( pair=(qubit, resonator), strength=0.01, ), ], ) simulator = QuantumSimulator(system)

In [ ]:

system.hamiltonian

system.hamiltonian

In [ ]:

control = Control(
    target=resonator.label,
    frequency=qubit.frequency,
    waveform=[2 * np.pi * 0.5] * 1000,
    durations=[1.0] * 1000,
)
control.plot()

control = Control( target=resonator.label, frequency=qubit.frequency, waveform=[2 * np.pi * 0.5] * 1000, durations=[1.0] * 1000, ) control.plot()

In [ ]:

result = simulator.mesolve(
    controls=[control],
    initial_state={
        "Q01": qt.basis(2, 0),
        "R01": qt.basis(10, 0),
    },
)

result = simulator.mesolve( controls=[control], initial_state={ "Q01": qt.basis(2, 0), "R01": qt.basis(10, 0), }, )

In [ ]:

result.plot_population_dynamics(qubit.label)
result.plot_population_dynamics(resonator.label)

result.plot_population_dynamics(qubit.label) result.plot_population_dynamics(resonator.label)

In [ ]:

result.display_bloch_sphere(qubit.label)

result.display_bloch_sphere(qubit.label)