def test_uccsd_hf_aer_qasm(self):
""" uccsd hf test with Aer qasm_simulator. """
try:
# pylint: disable=import-outside-toplevel
from qiskit import Aer
except Exception as ex: # pylint: disable=broad-except
self.skipTest(
"Aer doesn't appear to be installed. Error: '{}'".format(
str(ex)))
return
backend = Aer.get_backend('qasm_simulator')
optimizer = SPSA(maxiter=200, last_avg=5)
solver = VQE(var_form=self.var_form,
optimizer=optimizer,
expectation=PauliExpectation(),
quantum_instance=QuantumInstance(
backend=backend,
seed_simulator=aqua_globals.random_seed,
seed_transpiler=aqua_globals.random_seed))
gsc = GroundStateEigensolver(self.fermionic_transformation, solver)
result = gsc.solve(self.driver)
self.assertAlmostEqual(result.energy, -1.138, places=2)
def test_excitation_preserving(self):
"""Test the excitation preserving wavefunction on a chemistry example."""
driver = HDF5Driver(self.get_resource_path('test_driver_hdf5.hdf5'))
fermionic_transformation = FermionicTransformation(
qubit_mapping=QubitMappingType.PARITY, two_qubit_reduction=False)
qubit_op, _ = fermionic_transformation.transform(driver)
optimizer = SLSQP(maxiter=100)
initial_state = HartreeFock(
fermionic_transformation.molecule_info['num_orbitals'],
fermionic_transformation.molecule_info['num_particles'],
qubit_mapping=fermionic_transformation._qubit_mapping,
two_qubit_reduction=fermionic_transformation._two_qubit_reduction)
wavefunction = ExcitationPreserving(qubit_op.num_qubits,
initial_state=initial_state)
solver = VQE(var_form=wavefunction,
optimizer=optimizer,
quantum_instance=QuantumInstance(
BasicAer.get_backend('statevector_simulator'),
seed_simulator=aqua_globals.random_seed,
seed_transpiler=aqua_globals.random_seed))
gsc = GroundStateEigensolver(fermionic_transformation, solver)
result = gsc.solve(driver)
self.assertAlmostEqual(result.energy, self.reference_energy, places=4)
def test_uccsd_hf_qpUCCD(self):
""" paired uccd test """
optimizer = SLSQP(maxiter=100)
initial_state = HartreeFock(
self.fermionic_transformation.molecule_info['num_orbitals'],
self.fermionic_transformation.molecule_info['num_particles'],
qubit_mapping=self.fermionic_transformation._qubit_mapping,
two_qubit_reduction=self.fermionic_transformation._two_qubit_reduction)
var_form = UCCSD(
num_orbitals=self.fermionic_transformation.molecule_info['num_orbitals'],
num_particles=self.fermionic_transformation.molecule_info['num_particles'],
active_occupied=None, active_unoccupied=None,
initial_state=initial_state,
qubit_mapping=self.fermionic_transformation._qubit_mapping,
two_qubit_reduction=self.fermionic_transformation._two_qubit_reduction,
num_time_slices=1,
shallow_circuit_concat=False,
method_doubles='pucc',
excitation_type='d'
)
solver = VQE(var_form=var_form, optimizer=optimizer,
quantum_instance=QuantumInstance(
backend=BasicAer.get_backend('statevector_simulator')))
gsc = GroundStateEigensolver(self.fermionic_transformation, solver)
result = gsc.solve(self.driver)
self.assertAlmostEqual(result.total_energies[0], self.reference_energy_pUCCD, places=6)
def test_vqe_uccsd(self):
""" Test VQE UCCSD case """
solver = VQEUCCSDFactory(
QuantumInstance(BasicAer.get_backend('statevector_simulator')))
calc = GroundStateEigensolver(self.transformation, solver)
res = calc.solve(self.driver)
self.assertAlmostEqual(res.energy, self.reference_energy, places=6)
def test_tapered_op(self):
""" tapered op test """
optimizer = SLSQP(maxiter=1000)
init_state = HartreeFock(
num_orbitals=self.fermionic_transformation.molecule_info['num_orbitals'],
qubit_mapping=self.fermionic_transformation._qubit_mapping,
two_qubit_reduction=self.fermionic_transformation._two_qubit_reduction,
num_particles=self.fermionic_transformation.molecule_info['num_particles'],
sq_list=self.z2_symmetries.sq_list)
var_form = UCCSD(
num_orbitals=self.fermionic_transformation.molecule_info['num_orbitals'],
num_particles=self.fermionic_transformation.molecule_info['num_particles'],
active_occupied=None,
active_unoccupied=None,
initial_state=init_state,
qubit_mapping=self.fermionic_transformation._qubit_mapping,
two_qubit_reduction=self.fermionic_transformation._two_qubit_reduction,
num_time_slices=1,
z2_symmetries=self.z2_symmetries)
solver = VQE(var_form=var_form, optimizer=optimizer,
quantum_instance=QuantumInstance(
backend=BasicAer.get_backend('statevector_simulator')))
gsc = GroundStateEigensolver(self.fermionic_transformation, solver)
result = gsc.solve(self.driver)
self.assertAlmostEqual(result.total_energies[0], self.reference_energy, places=6)
def test_npme(self):
""" Test NumPyMinimumEigensolver """
solver = NumPyMinimumEigensolverFactory()
calc = GroundStateEigensolver(self.transformation, solver)
res = calc.solve(self.driver)
self.assertAlmostEqual(res.total_energies[0],
self.reference_energy,
places=6)
def test_npme_with_default_filter(self):
""" Test NumPyMinimumEigensolver with default filter """
solver = NumPyMinimumEigensolverFactory(
use_default_filter_criterion=True)
calc = GroundStateEigensolver(self.transformation, solver)
res = calc.solve(self.driver)
self.assertAlmostEqual(res.total_energies[0],
self.reference_energy,
places=6)
def _setup_evaluation_operators(self):
# first we run a ground state calculation
solver = VQEUCCSDFactory(QuantumInstance(BasicAer.get_backend('statevector_simulator')))
calc = GroundStateEigensolver(self.transformation, solver)
res = calc.solve(self.driver)
# now we decide that we want to evaluate another operator
# for testing simplicity, we just use some pre-constructed auxiliary operators
_, aux_ops = self.transformation.transform(self.driver)
return calc, res, aux_ops
def test_uccsd_hf(self):
""" uccsd hf test """
backend = BasicAer.get_backend('statevector_simulator')
solver = VQE(var_form=self.var_form, optimizer=self.optimizer,
quantum_instance=QuantumInstance(backend=backend))
gsc = GroundStateEigensolver(self.fermionic_transformation, solver)
result = gsc.solve(self.driver)
self.assertAlmostEqual(result.total_energies[0], self.reference_energy, places=6)
def test_with_numpy_minimum_eigensolver(self):
""" Test with NumPyMinimumEigensolver """
bosonic_transformation = BosonicTransformation(
qubit_mapping=BosonicQubitMappingType.DIRECT,
transformation_type=BosonicTransformationType.HARMONIC,
basis_size=2,
truncation=2)
solver = NumPyMinimumEigensolverFactory(use_default_filter_criterion=True)
gsc = GroundStateEigensolver(bosonic_transformation, solver)
result = gsc.solve(self.driver)
self.assertAlmostEqual(result.computed_vibronic_energies[0],
self.reference_energy, places=4)
def test_uccsd_hf_qasm(self):
""" uccsd hf test with qasm_simulator. """
backend = BasicAer.get_backend('qasm_simulator')
optimizer = SPSA(maxiter=200, last_avg=5)
solver = VQE(var_form=self.var_form, optimizer=optimizer,
expectation=PauliExpectation(),
quantum_instance=QuantumInstance(backend=backend,
seed_simulator=aqua_globals.random_seed,
seed_transpiler=aqua_globals.random_seed))
gsc = GroundStateEigensolver(self.fermionic_transformation, solver)
result = gsc.solve(self.driver)
self.assertAlmostEqual(result.total_energies[0], -1.138, places=2)
def test_aux_ops_reusability(self):
""" Test that the auxiliary operators can be reused """
# Regression test against #1475
solver = NumPyMinimumEigensolverFactory()
calc = GroundStateEigensolver(self.transformation, solver)
modes = 4
h_1 = np.eye(modes, dtype=np.complex)
h_2 = np.zeros((modes, modes, modes, modes))
aux_ops = [FermionicOperator(h_1, h_2)]
aux_ops_copy = copy.deepcopy(aux_ops)
_ = calc.solve(self.driver, aux_ops)
assert all([a == b for a, b in zip(aux_ops, aux_ops_copy)])
def test_vqe_uvccsd_factory(self):
""" Test with VQE plus UCCSD """
bosonic_transformation = BosonicTransformation(
qubit_mapping=BosonicQubitMappingType.DIRECT,
transformation_type=BosonicTransformationType.HARMONIC,
basis_size=2,
truncation=2)
optimizer = COBYLA(maxiter=5000)
solver = VQEUVCCSDFactory(QuantumInstance(BasicAer.get_backend('statevector_simulator')),
optimizer=optimizer)
gsc = GroundStateEigensolver(bosonic_transformation, solver)
result = gsc.solve(self.driver)
self.assertAlmostEqual(result.computed_vibronic_energies[0], self.reference_energy,
places=1)
def test_uccsd_hf_aer_statevector(self):
""" uccsd hf test with Aer statevector """
try:
# pylint: disable=import-outside-toplevel
from qiskit import Aer
except Exception as ex: # pylint: disable=broad-except
self.skipTest("Aer doesn't appear to be installed. Error: '{}'".format(str(ex)))
return
backend = Aer.get_backend('statevector_simulator')
solver = VQE(var_form=self.var_form, optimizer=self.optimizer,
quantum_instance=QuantumInstance(backend=backend))
gsc = GroundStateEigensolver(self.fermionic_transformation, solver)
result = gsc.solve(self.driver)
self.assertAlmostEqual(result.total_energies[0], self.reference_energy, places=6)
def _run_driver(driver,
transformation=TransformationType.FULL,
qubit_mapping=QubitMappingType.JORDAN_WIGNER,
two_qubit_reduction=False,
freeze_core=True):
fermionic_transformation = \
FermionicTransformation(transformation=transformation,
qubit_mapping=qubit_mapping,
two_qubit_reduction=two_qubit_reduction,
freeze_core=freeze_core,
orbital_reduction=[])
solver = NumPyMinimumEigensolver()
gsc = GroundStateEigensolver(fermionic_transformation, solver)
result = gsc.solve(driver)
return result
def test_uccsd_hf_aer_qasm_snapshot(self):
""" uccsd hf test with Aer qasm_simulator snapshot. """
try:
# pylint: disable=import-outside-toplevel
from qiskit import Aer
backend = Aer.get_backend('qasm_simulator')
except ImportError as ex: # pylint: disable=broad-except
self.skipTest(
"Aer doesn't appear to be installed. Error: '{}'".format(
str(ex)))
return
optimizer = SPSA(maxiter=200, last_avg=5)
solver = VQE(var_form=self.var_form,
optimizer=optimizer,
expectation=AerPauliExpectation(),
quantum_instance=QuantumInstance(backend=backend))
gsc = GroundStateEigensolver(self.fermionic_transformation, solver)
result = gsc.solve(self.driver)
self.assertAlmostEqual(result.total_energies[0],
self.reference_energy,
places=3)