qamomile.observable

Overview¶

Functions¶

X [source]¶

def X(index: int) -> Hamiltonian

Creates a Pauli X operator for a specified qubit.

Parameters:

Returns:

Hamiltonian — A Pauli X Hamiltonian operator acting on the specified qubit.

Example:

>>> X0 = X(0)
>>> print(X0)
Hamiltonian((X0,): 1.0)

Y [source]¶

def Y(index: int) -> Hamiltonian

Creates a Pauli Y operator for a specified qubit.

Parameters:

Returns:

Hamiltonian — A Pauli Y Hamiltonian operator acting on the specified qubit.

Example:

>>> Y1 = Y(1)
>>> print(Y1)
Hamiltonian((Y1,): 1.0)

Z [source]¶

def Z(index: int) -> Hamiltonian

Creates a Pauli Z operator for a specified qubit.

Parameters:

Returns:

Hamiltonian — A Pauli Z Hamiltonian operator acting on the specified qubit.

Example:

>>> Z2 = Z(2)
>>> print(Z2)
Hamiltonian((Z2,): 1.0)

Classes¶

Hamiltonian [source]¶

class Hamiltonian

Represents a quantum Hamiltonian as a sum of Pauli operator products.

The Hamiltonian is stored as a dictionary where keys are tuples of PauliOperators and values are their corresponding coefficients.

Example:

>>> H = Hamiltonian()
>>> H.add_term((PauliOperator(Pauli.X, 0), PauliOperator(Pauli.Y, 1)), 0.5)
>>> H.add_term((PauliOperator(Pauli.Z, 2),), 1.0)
>>> print(H.terms)
{(X0, Y1): 0.5, (Z2,): 1.0}

Constructor¶

def __init__(self, num_qubits: int | None = None) -> None

Attributes¶

• constant: float | complex

• num_qubits: int Calculates the number of qubits in the Hamiltonian.

• terms: dict[tuple[PauliOperator, ...], complex] Getter for the terms of the Hamiltonian.

Methods¶

add_term¶

def add_term(self, operators: tuple[PauliOperator, ...], coeff: float | complex)

Adds a term to the Hamiltonian.

This method adds a product of Pauli operators with a given coefficient to the Hamiltonian. If the term already exists, the coefficients are summed.

Parameters:

Example:

>>> H = Hamiltonian()
>>> H.add_term((PauliOperator(Pauli.X, 0), PauliOperator(Pauli.Y, 1)), 0.5)
>>> H.add_term((PauliOperator(Pauli.X, 0), PauliOperator(Pauli.Y, 1)), 0.5j)
>>> print(H.terms)
{(X0, Y1): (0.5+0.5j)}

identity¶

@classmethod
def identity(
    cls,
    coeff: float | complex = 1.0,
    num_qubits: int | None = None,
) -> Hamiltonian

Create a scalar times identity Hamiltonian.

remap_qubits¶

def remap_qubits(self, qubit_map: dict[int, int]) -> Hamiltonian

Remap qubit indices according to the given mapping.

This is used to translate Pauli indices (logical indices within an expval call) to physical qubit indices in the actual quantum circuit.

Parameters:

Returns:

Hamiltonian — New Hamiltonian with remapped qubit indices.

single_pauli¶

@classmethod
def single_pauli(cls, pauli: Pauli, index: int, coeff: float | complex = 1.0) -> Hamiltonian

Create a single Pauli term Hamiltonian.

to_latex¶

def to_latex(self) -> str

Converts the Hamiltonian to a LaTeX representation.

This function does not add constant term when we show the Hamiltonian. This function does not add $ symbols.

Returns:

str — A LaTeX representation of the Hamiltonian.

import qamomile.core.operator as qm_o
import IPython.display as ipd

h = qm_o.Hamiltonian()
h += -qm_o.X(0) * qm_o.Y(1) - 2.0 * qm_o.Z(0) * qm_o.Z(1)

# Show the Hamiltonian in LaTeX at Jupyter Notebook
ipd.display(ipd.Latex("$" + h.to_latex() + "$"))

zero¶

@classmethod
def zero(cls, num_qubits: int | None = None) -> Hamiltonian

Create a zero Hamiltonian.

Pauli [source]¶

class Pauli(enum.Enum)

Enum class for Pauli operators.

Attributes¶

• I

• X

• Y

• Z

PauliOperator [source]¶

class PauliOperator

Represents a single Pauli operator acting on a specific qubit.

Example:

>>> X0 = PauliOperator(Pauli.X, 0)
>>> print(X0)
X0

Constructor¶

def __init__(self, pauli: Pauli, index: int) -> None

Attributes¶

• index: int

• pauli: Pauli

qamomile.observable.hamiltonian¶

This module provides the intermediate representation of Hamiltonian for quantum systems.

It defines classes and functions to create and manipulate Pauli operators and Hamiltonians, which are fundamental in quantum mechanics and quantum computing.

Key Components:

• Pauli: An enumeration of Pauli operators (X, Y, Z, I).

• PauliOperator: A class representing a single Pauli operator acting on a specific qubit.

• Hamiltonian: A class representing a quantum Hamiltonian as a sum of Pauli operator products.

Usage:

from qamomile.operator.hamiltonian import X, Y, Z, Hamiltonian

Create Pauli operators¶

X0 = X(0) # Pauli X operator on qubit 0 Y1 = Y(1) # Pauli Y operator on qubit 1

Create a Hamiltonian¶

H = Hamiltonian() H.add_term((X0, Y1), 0.5) # Add term 0.5 * X0 * Y1 to the Hamiltonian H.add_term((Z(2),), 1.0) # Add term 1.0 * Z2 to the Hamiltonian

Access Hamiltonian properties¶

print(H.terms) print(H.num_qubits)

Features:

• Factory methods for common Hamiltonians (zero, identity, single_pauli)

• Iterator protocol for looping over terms

• Qubit index remapping for circuit integration

• Modern type annotations (PEP 604)

Overview¶

Functions¶

X [source]¶

def X(index: int) -> Hamiltonian

Creates a Pauli X operator for a specified qubit.

Parameters:

Returns:

Hamiltonian — A Pauli X Hamiltonian operator acting on the specified qubit.

Example:

>>> X0 = X(0)
>>> print(X0)
Hamiltonian((X0,): 1.0)

Y [source]¶

def Y(index: int) -> Hamiltonian

Creates a Pauli Y operator for a specified qubit.

Parameters:

Returns:

Hamiltonian — A Pauli Y Hamiltonian operator acting on the specified qubit.

Example:

>>> Y1 = Y(1)
>>> print(Y1)
Hamiltonian((Y1,): 1.0)

Z [source]¶

def Z(index: int) -> Hamiltonian

Creates a Pauli Z operator for a specified qubit.

Parameters:

Returns:

Hamiltonian — A Pauli Z Hamiltonian operator acting on the specified qubit.

Example:

>>> Z2 = Z(2)
>>> print(Z2)
Hamiltonian((Z2,): 1.0)

multiply_pauli_same_qubit [source]¶

def multiply_pauli_same_qubit(pauli1: PauliOperator, pauli2: PauliOperator) -> tuple[PauliOperator, complex]

Multiplies two Pauli operators acting on the same qubit.

Parameters:

Returns:

tuple[PauliOperator, complex] — tuple[PauliOperator, complex]: A tuple containing the resulting Pauli operator and the complex coefficient.

Raises:

• ValueError — If the Pauli operators act on different qubits.

Example:

>>> X0 = PauliOperator(Pauli.X, 0)
>>> Y0 = PauliOperator(Pauli.Y, 0)
>>> Z0 = PauliOperator(Pauli.Z, 0)
>>> multiply_pauli_same_qubit(X0, Y0)
(Z0, 1j)

simplify_pauliop_terms [source]¶

def simplify_pauliop_terms(term: tuple[PauliOperator, ...]) -> tuple[tuple[PauliOperator, ...], complex]

Simplifies a tuple of Pauli operators by combining operators acting on the same qubit.

This function performs canonicalization by:

1. Multiplying Pauli operators on the same qubit using the Pauli algebra

2. Removing identity operators from the result

3. Tracking the accumulated phase factor from the multiplications

Parameters:

Returns:

tuple[tuple[PauliOperator, ...], complex] — tuple[tuple[PauliOperator,...],complex]: A tuple containing the simplified Pauli operators and the phase factor.

Example:

>>> X0 = PauliOperator(Pauli.X, 0)
>>> Y0 = PauliOperator(Pauli.Y, 0)
>>> Z1 = PauliOperator(Pauli.Z, 1)
>>> simplify_pauliop_terms((X0, Y0, Z1))
((Z0, Z1), 1j)

Classes¶

Hamiltonian [source]¶

class Hamiltonian

Represents a quantum Hamiltonian as a sum of Pauli operator products.

The Hamiltonian is stored as a dictionary where keys are tuples of PauliOperators and values are their corresponding coefficients.

Example:

>>> H = Hamiltonian()
>>> H.add_term((PauliOperator(Pauli.X, 0), PauliOperator(Pauli.Y, 1)), 0.5)
>>> H.add_term((PauliOperator(Pauli.Z, 2),), 1.0)
>>> print(H.terms)
{(X0, Y1): 0.5, (Z2,): 1.0}

Constructor¶

def __init__(self, num_qubits: int | None = None) -> None

Attributes¶

• constant: float | complex

• num_qubits: int Calculates the number of qubits in the Hamiltonian.

• terms: dict[tuple[PauliOperator, ...], complex] Getter for the terms of the Hamiltonian.

Methods¶

add_term¶

def add_term(self, operators: tuple[PauliOperator, ...], coeff: float | complex)

Adds a term to the Hamiltonian.

This method adds a product of Pauli operators with a given coefficient to the Hamiltonian. If the term already exists, the coefficients are summed.

Parameters:

Example:

>>> H = Hamiltonian()
>>> H.add_term((PauliOperator(Pauli.X, 0), PauliOperator(Pauli.Y, 1)), 0.5)
>>> H.add_term((PauliOperator(Pauli.X, 0), PauliOperator(Pauli.Y, 1)), 0.5j)
>>> print(H.terms)
{(X0, Y1): (0.5+0.5j)}

identity¶

@classmethod
def identity(
    cls,
    coeff: float | complex = 1.0,
    num_qubits: int | None = None,
) -> Hamiltonian

Create a scalar times identity Hamiltonian.

remap_qubits¶

def remap_qubits(self, qubit_map: dict[int, int]) -> Hamiltonian

Remap qubit indices according to the given mapping.

This is used to translate Pauli indices (logical indices within an expval call) to physical qubit indices in the actual quantum circuit.

Parameters:

Returns:

Hamiltonian — New Hamiltonian with remapped qubit indices.

single_pauli¶

@classmethod
def single_pauli(cls, pauli: Pauli, index: int, coeff: float | complex = 1.0) -> Hamiltonian

Create a single Pauli term Hamiltonian.

to_latex¶

def to_latex(self) -> str

Converts the Hamiltonian to a LaTeX representation.

This function does not add constant term when we show the Hamiltonian. This function does not add $ symbols.

Returns:

str — A LaTeX representation of the Hamiltonian.

import qamomile.core.operator as qm_o
import IPython.display as ipd

h = qm_o.Hamiltonian()
h += -qm_o.X(0) * qm_o.Y(1) - 2.0 * qm_o.Z(0) * qm_o.Z(1)

# Show the Hamiltonian in LaTeX at Jupyter Notebook
ipd.display(ipd.Latex("$" + h.to_latex() + "$"))

zero¶

@classmethod
def zero(cls, num_qubits: int | None = None) -> Hamiltonian

Create a zero Hamiltonian.

Pauli [source]¶

class Pauli(enum.Enum)

Enum class for Pauli operators.

Attributes¶

• I

• X

• Y

• Z

PauliOperator [source]¶

class PauliOperator

Represents a single Pauli operator acting on a specific qubit.

Example:

>>> X0 = PauliOperator(Pauli.X, 0)
>>> print(X0)
X0

Constructor¶

def __init__(self, pauli: Pauli, index: int) -> None

Attributes¶

• index: int

• pauli: Pauli