# 必要なモジュールをインポート
import numpy as np
from qiskit import QuantumCircuit, Aer
from qiskit.opflow import I, X, Z
from qiskit.circuit import Parameter
from qiskit.quantum_info import Statevector
from scipy.optimize import minimize

# ハミルトニアンを定義
hamiltonian = ((-3.8505 * I^I) + (-0.2288 * X^I) + (-0.2288 * I^X) 
            + (-1.0466 * Z^I) + (-1.0466 * I^Z) + (0.2613 * X^X) 
            + (0.2288 * X^Z) + (0.2288 * Z^X) + (0.2356 * Z^Z))/2