Backends

Backends allow execution of Qadence abstract quantum circuits. They could be chosen from a variety of simulators, emulators and hardware and can enable circuit differentiability. The primary way to interact and configure a backend is via the high-level API QuantumModel.

Not all backends are equivalent

Not all backends support the same set of operations, especially while executing analog blocks. Qadence will throw descriptive errors in such cases.

Execution backends

PyQTorch: An efficient, large-scale simulator designed for quantum machine learning, seamlessly integrated with the popular PyTorch deep learning framework for automatic differentiability. It also offers analog computing for time-(in)dependent pulses. See PyQTorchBackend.

Pulser: A Python library for pulse-level/analog control of neutral atom devices. Execution via QuTiP. See PulserBackend.

More: Proprietary Qadence extensions provide more high-performance backends based on tensor networks or differentiation engines. For more enquiries, please contact: info@pasqal.com.

Differentiation backend

The DifferentiableBackend class enables different differentiation modes for the given backend. This can be chosen from two types:

• Automatic differentiation (AD): available for PyTorch based backends (PyQTorch).
• Parameter Shift Rules (PSR): available for all backends. See this section for more information on differentiability and PSR.

In practice, only a diff_mode should be provided in the QuantumModel. Please note that diff_mode defaults to None:

import sympy
import torch
from qadence import Parameter, RX, RZ, Z, CNOT, QuantumCircuit, QuantumModel, chain, BackendName, DiffMode

x = Parameter("x", trainable=False)
y = Parameter("y", trainable=False)
fm = chain(
    RX(0, 3 * x),
    RX(0, x),
    RZ(1, sympy.exp(y)),
    RX(0, 3.14),
    RZ(1, "theta")
)

ansatz = CNOT(0, 1)
block = chain(fm, ansatz)

circuit = QuantumCircuit(2, block)

observable = Z(0)

# DiffMode.GPSR is available for any backend.
# DiffMode.AD is only available for natively differentiable backends.
model = QuantumModel(circuit, observable, backend=BackendName.PYQTORCH, diff_mode=DiffMode.GPSR)

# Get some values for the feature parameters.
values = {"x": (x := torch.tensor([0.5], requires_grad=True)), "y": torch.tensor([0.1])}

# Compute expectation.
exp = model.expectation(values)

# Differentiate the expectation wrt x.
dexp_dx = torch.autograd.grad(exp, x, torch.ones_like(exp))

dexp_dx = (tensor([3.6398]),)

Low-level backend_factory interface

Every backend in Qadence inherits from the abstract Backend class: Backend and implement the following methods:

• run: propagate the initial state according to the quantum circuit and return the final wavefunction object.
• sample: sample from a circuit.
• expectation: computes the expectation of a circuit given an observable.
• convert: convert the abstract QuantumCircuit object to its backend-native representation including a backend specific parameter embedding function.

Backends are purely functional objects which take as input the values for the circuit parameters and return the desired output from a call to a method. In order to use a backend directly, embedded parameters must be supplied as they are returned by the backend specific embedding function.

Here is a simple demonstration of the use of the PyQTorch backend to execute a circuit in non-differentiable mode:

from qadence import QuantumCircuit, FeatureParameter, RX, RZ, CNOT, hea, chain

# Construct a feature map.
x = FeatureParameter("x")
z = FeatureParameter("y")
fm = chain(RX(0, 3 * x), RZ(1, z), CNOT(0, 1))

# Construct a circuit with an hardware-efficient ansatz.
circuit = QuantumCircuit(3, fm, hea(3,1))

The abstract QuantumCircuit can now be converted to its native representation via the PyQTorch backend.

from qadence import backend_factory

# Use only PyQtorch in non-differentiable mode:
backend = backend_factory("pyqtorch")

# The `Converted` object
# (contains a `ConvertedCircuit` with the original and native representation)
conv = backend.convert(circuit)

conv.circuit.original = ChainBlock(0,1,2)
├── ChainBlock(0,1)
│   ├── RX(0) [params: ['3*x']]
│   ├── RZ(1) [params: ['y']]
│   └── CNOT(0, 1)
└── ChainBlock(0,1,2) [tag: HEA]
    ├── ChainBlock(0,1,2)
    │   ├── KronBlock(0,1,2)
    │   │   ├── RX(0) [params: ['theta_0']]
    │   │   ├── RX(1) [params: ['theta_1']]
    │   │   └── RX(2) [params: ['theta_2']]
    │   ├── KronBlock(0,1,2)
    │   │   ├── RY(0) [params: ['theta_3']]
    │   │   ├── RY(1) [params: ['theta_4']]
    │   │   └── RY(2) [params: ['theta_5']]
    │   └── KronBlock(0,1,2)
    │       ├── RX(0) [params: ['theta_6']]
    │       ├── RX(1) [params: ['theta_7']]
    │       └── RX(2) [params: ['theta_8']]
    └── ChainBlock(0,1,2)
        ├── KronBlock(0,1)
        │   └── CNOT(0, 1)
        └── KronBlock(1,2)
            └── CNOT(1, 2)
conv.circuit.native = QuantumCircuit(
  (operations): ModuleList(
    (0): Sequence(
      (operations): ModuleList(
        (0): Sequence(
          (operations): ModuleList(
            (0): RX(target: (0,), param: 830af221-5b50-411c-aed1-1ba91a919a85)
            (1): RZ(target: (1,), param: aaf2f95b-cb3e-4c1c-a5c3-66fe442e4b77)
            (2): CNOT(control: (0,), target: (1,))
          )
        )
        (1): Sequence(
          (operations): ModuleList(
            (0): Sequence(
              (operations): ModuleList(
                (0): Merge(
                  (operations): ModuleList(
                    (0): RX(target: (0,), param: b40d1b0a-3570-4025-980f-eae1ead879a4)
                    (1): RY(target: (0,), param: 178c307c-f0d0-46e4-ac7a-d0bcc907b148)
                    (2): RX(target: (0,), param: e3e6d0e6-0f5e-4555-a531-b674622b3bbd)
                  )
                )
                (1): Merge(
                  (operations): ModuleList(
                    (0): RX(target: (1,), param: 206305a8-fab3-4950-bb1e-3785703b8091)
                    (1): RY(target: (1,), param: 38384b4f-f67f-46dc-9809-fee6d0959745)
                    (2): RX(target: (1,), param: ce9f0234-ebd7-46fd-a681-c72e7d8e7cb9)
                  )
                )
                (2): Merge(
                  (operations): ModuleList(
                    (0): RX(target: (2,), param: b5b0d188-d247-4b5b-9aa9-311b888dfb1d)
                    (1): RY(target: (2,), param: 622c21b7-511f-45a7-908b-192ea53ea078)
                    (2): RX(target: (2,), param: 68f78757-3a9d-4306-8a9f-28ec3c0b1e55)
                  )
                )
              )
            )
            (1): Sequence(
              (operations): ModuleList(
                (0): Sequence(
                  (operations): ModuleList(
                    (0): CNOT(control: (0,), target: (1,))
                  )
                )
                (1): Sequence(
                  (operations): ModuleList(
                    (0): CNOT(control: (1,), target: (2,))
                  )
                )
              )
            )
          )
        )
      )
    )
  )
)

Additionally, Converted contains all fixed and variational parameters, as well as an embedding function which accepts feature parameters to construct a dictionary of circuit native parameters. These are needed as each backend uses a different representation of the circuit parameters:

import torch

# Contains fixed parameters and variational (from the HEA)
conv.params

inputs = {"x": torch.tensor([1., 1.]), "y":torch.tensor([2., 2.])}

# get all circuit parameters (including feature params)
embedded = conv.embedding_fn(conv.params, inputs)

conv.params = {
  theta_5: tensor([0.5045], requires_grad=True)
  theta_6: tensor([0.7569], requires_grad=True)
  theta_8: tensor([0.5915], requires_grad=True)
  theta_1: tensor([0.2272], requires_grad=True)
  theta_4: tensor([0.1270], requires_grad=True)
  theta_0: tensor([0.6965], requires_grad=True)
  theta_7: tensor([0.7917], requires_grad=True)
  theta_2: tensor([0.4572], requires_grad=True)
  theta_3: tensor([0.3075], requires_grad=True)
}
embedded = {
  830af221-5b50-411c-aed1-1ba91a919a85: tensor([3., 3.], grad_fn=<ViewBackward0>)
  aaf2f95b-cb3e-4c1c-a5c3-66fe442e4b77: tensor([2., 2.])
  b40d1b0a-3570-4025-980f-eae1ead879a4: tensor([0.6965], grad_fn=<ViewBackward0>)
  178c307c-f0d0-46e4-ac7a-d0bcc907b148: tensor([0.3075], grad_fn=<ViewBackward0>)
  e3e6d0e6-0f5e-4555-a531-b674622b3bbd: tensor([0.7569], grad_fn=<ViewBackward0>)
  206305a8-fab3-4950-bb1e-3785703b8091: tensor([0.2272], grad_fn=<ViewBackward0>)
  38384b4f-f67f-46dc-9809-fee6d0959745: tensor([0.1270], grad_fn=<ViewBackward0>)
  ce9f0234-ebd7-46fd-a681-c72e7d8e7cb9: tensor([0.7917], grad_fn=<ViewBackward0>)
  b5b0d188-d247-4b5b-9aa9-311b888dfb1d: tensor([0.4572], grad_fn=<ViewBackward0>)
  622c21b7-511f-45a7-908b-192ea53ea078: tensor([0.5045], grad_fn=<ViewBackward0>)
  68f78757-3a9d-4306-8a9f-28ec3c0b1e55: tensor([0.5915], grad_fn=<ViewBackward0>)
}

With the embedded parameters, QuantumModel methods are accessible:

output = backend.run(conv.circuit, embedded)
print(f"{output = }")

output = tensor([[ 0.2903+0.0211j,  0.1016-0.1596j,  0.2012+0.2426j, -0.1817+0.4495j,
         -0.4665-0.2817j, -0.3051+0.1800j,  0.0967+0.1646j, -0.1523+0.2511j],
        [ 0.2903+0.0211j,  0.1016-0.1596j,  0.2012+0.2426j, -0.1817+0.4495j,
         -0.4665-0.2817j, -0.3051+0.1800j,  0.0967+0.1646j, -0.1523+0.2511j]],
       grad_fn=<TBackward0>)

Lower-level: the Backend representation

If there is a requirement to work with a specific backend, it is possible to access directly the native circuit. For example, should one wish to use PyQtorch noise features directly instead of using the NoiseHandler interface from Qadence:

from pyqtorch.noise import Depolarizing

inputs = {"x": torch.rand(1), "y":torch.rand(1)}
embedded = conv.embedding_fn(conv.params, inputs)

# Define a noise channel on qubit 0
noise = Depolarizing(0, error_probability=0.1)

# Add noise to circuit
conv.circuit.native.operations.append(noise)

When running With noise, one can see that the output is a density matrix:

density_result = backend.run(conv.circuit, embedded)
print(density_result.shape)

torch.Size([1, 8, 8])