def cost_func(theta):
    # 量子状態シミュレータを利用
    backend = Aer.get_backend('statevector_simulator')
    # パラメータ付き量子回路をパラメータをθとして計算
    job = backend.run(PQC(theta), shots = 1024)
    # 計算結果を取得
    result = job.result()
    # 計算結果の状態ベクトルを取得
    q_state = Statevector(result.get_statevector())
    # 期待値を返す
    return np.real(q_state.expectation_value(hamiltonian))