qiskit_bell_state_circuit_with_aer_estimator_simulation.py

1from qiskit import QuantumCircuit
2from qiskit.quantum_info import SparsePauliOp
3from qiskit_aer.primitives import Estimator
4from qiskit.visualization import plot_histogram
5
6# Create a new circuit with two qubits
7# by default, the qubits are set to the state |00>
8circuit = QuantumCircuit(2)
9
10# Add a Hadamard gate to qubit 0
11circuit.h(0)
12
13# Add a CX (CNOT) gate on control qubit 0 and target qubit 1
14# putting the qubits in a Bell state
15circuit.cx(0, 1)
16
17# Set up the observables to measure
18# In this case, we measure the ZZ expectation value
19observable = SparsePauliOp("ZZ")
20
21# Use the Estimator primitive to get the expectation value
22estimator = Estimator()
23job = estimator.run(circuit, observable)
24result = job.result()
25
26# Print the result
27print(f" > Expectation value: {result.values[0]}")
28print(f" > Circuit: \n{circuit.draw()}")