From c9ca9a58ad3db6abdc0bd23d83d8ad50bf5af2e1 Mon Sep 17 00:00:00 2001 From: Divyanshu Singh <55018955+divshacker@users.noreply.github.com> Date: Fri, 15 Jan 2021 17:43:52 +0530 Subject: [PATCH 1/2] Delete 01_algorithms_introduction.ipynb --- .../01_algorithms_introduction.ipynb | 355 ------------------ 1 file changed, 355 deletions(-) delete mode 100644 tutorials/algorithms/01_algorithms_introduction.ipynb diff --git a/tutorials/algorithms/01_algorithms_introduction.ipynb b/tutorials/algorithms/01_algorithms_introduction.ipynb deleted file mode 100644 index 808ccda3a..000000000 --- a/tutorials/algorithms/01_algorithms_introduction.ipynb +++ /dev/null @@ -1,355 +0,0 @@ -{ - "cells": [ - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "# An Introduction to Algorithms in Qiskit\n", - "\n", - "This is an introduction to algorithms in Qiskit and provides a high-level overview to help understand the various aspects of the functionality to get started. Other tutorials will provide more in-depth material, on given algorithms, and ways to use them etc." - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "## How is the algorithm library structured?\n", - "\n", - "Qiskit provides a number of [Algorithms](https://qiskit.org/documentation/apidoc/qiskit.aqua.algorithms.html) and they are grouped by category according to the task they can perform. For instance `Minimum Eigensolvers` to find the smallest eigen value of an operator, for example ground state energy of a chemistry Hamiltonian or a solution to an optimization problem when expressed as an Ising Hamiltonian. There are `Classifiers` for machine learning classification problems and `Amplitude Estimators` for value estimation that can be used say in financial applications. The full set of categories can be seen in the Algorithms documentation link above.\n", - "\n", - "Algorithms are configurable and often part of the configuration will be in the form of smaller building blocks, of which different instances of the building block type can be given. For instance with `VQE`, the Variational Quantum Eigensolver, it takes a trial wavefunction, in the form of a `QuantumCircuit` and a classical optimizer among other things. Such building blocks can be found as [Components](https://qiskit.org/documentation/apidoc/qiskit.aqua.components.html) and as circuits from the [Circuit Library](https://qiskit.org/documentation/apidoc/circuit_library.html).\n", - "\n", - "Let's take a look at an example to construct a VQE instance. Here `TwoLocal` is the variational form (trial wavefunction), a parameterized circuit which can be varied, and `SLSQP` a classical optimizer. These are created as separate instances and passed to VQE when it is constructed. Trying, for example, a different classical optimizer, or variational form is simply a case of creating an instance of the one you want and passing it into VQE." - ] - }, - { - "cell_type": "code", - "execution_count": 1, - "metadata": {}, - "outputs": [], - "source": [ - "from qiskit.aqua.algorithms import VQE\n", - "from qiskit.aqua.components.optimizers import SLSQP\n", - "from qiskit.circuit.library import TwoLocal\n", - "\n", - "num_qubits = 2\n", - "ansatz = TwoLocal(num_qubits, 'ry', 'cz')\n", - "opt = SLSQP(maxiter=1000)\n", - "vqe = VQE(var_form=ansatz, optimizer=opt)" - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "Let's draw the ansatz so we can see it's a QuantumCircuit where θ\\[0\\] through θ\\[7\\] will be the parameters that are varied as VQE optimizer finds the minimum eigenvalue. We'll come back to the parameters later in a working example below." - ] - }, - { - "cell_type": "code", - "execution_count": 2, - "metadata": {}, - "outputs": [ - { - "data": { - "text/html": [ - "
     ┌──────────┐   ┌──────────┐   ┌──────────┐   ┌──────────┐\n",
-       "q_0: ┤ RY(θ[0]) ├─■─┤ RY(θ[2]) ├─■─┤ RY(θ[4]) ├─■─┤ RY(θ[6]) ├\n",
-       "     ├──────────┤ │ ├──────────┤ │ ├──────────┤ │ ├──────────┤\n",
-       "q_1: ┤ RY(θ[1]) ├─■─┤ RY(θ[3]) ├─■─┤ RY(θ[5]) ├─■─┤ RY(θ[7]) ├\n",
-       "     └──────────┘   └──────────┘   └──────────┘   └──────────┘
" - ], - "text/plain": [ - " ┌──────────┐ ┌──────────┐ ┌──────────┐ ┌──────────┐\n", - "q_0: ┤ RY(θ[0]) ├─■─┤ RY(θ[2]) ├─■─┤ RY(θ[4]) ├─■─┤ RY(θ[6]) ├\n", - " ├──────────┤ │ ├──────────┤ │ ├──────────┤ │ ├──────────┤\n", - "q_1: ┤ RY(θ[1]) ├─■─┤ RY(θ[3]) ├─■─┤ RY(θ[5]) ├─■─┤ RY(θ[7]) ├\n", - " └──────────┘ └──────────┘ └──────────┘ └──────────┘" - ] - }, - "execution_count": 2, - "metadata": {}, - "output_type": "execute_result" - } - ], - "source": [ - "ansatz.draw()" - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "But more is needed before we can run the algorithm so let's get to that next." - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "## How to run an algorithm?\n", - "\n", - "In order to run an algorithm we need to have backend, a simulator or real device, on which the circuits that comprise the algorithm can be run. So for example we can use the `statevector_simulator` from the BasicAer provider." - ] - }, - { - "cell_type": "code", - "execution_count": 3, - "metadata": {}, - "outputs": [], - "source": [ - "from qiskit import BasicAer\n", - "\n", - "backend = BasicAer.get_backend('statevector_simulator')" - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "Now a backend on its own does not have information on how you might want to run the circuits etc. For example how many shots, do you want a noise model, even options around transpiling the circuits. For this Qiskit Aqua has a [QuantumInstance](https://qiskit.org/documentation/stubs/qiskit.aqua.QuantumInstance.html) which is provided both the backend as well as various settings around the circuit processing and execution so for instance:" - ] - }, - { - "cell_type": "code", - "execution_count": 4, - "metadata": {}, - "outputs": [], - "source": [ - "from qiskit.aqua import QuantumInstance\n", - "\n", - "backend = BasicAer.get_backend('qasm_simulator')\n", - "quantum_instance = QuantumInstance(backend=backend, shots=800, seed_simulator=99) " - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "Note: if you provide the backend directly then internally a QuantumInstance will be created from it, with default settings, so at all times the algorithms are working through a QuantumInstance.\n", - "\n", - "So now we would be able to call the [run()](https://qiskit.org/documentation/stubs/qiskit.aqua.algorithms.VQE.run.html) method, which is common to all algorithms and returns a result specific for the algorithm. In this case since VQE is a MinimumEigensolver we could use the [compute_mininum_eigenvalue()](https://qiskit.org/documentation/stubs/qiskit.aqua.algorithms.VQE.compute_minimum_eigenvalue.html) method. The latter is the interface of choice for the application modules, such as Chemistry and Optimization, in order that they can work interchangeably with any algorithm within the specific category." - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "## A complete working example\n", - "\n", - "Let's put what we have learned from above together and create a complete working example. VQE will find the minimum eigenvalue, i.e. minimum energy value of a Hamilitonian operator and hence we need such an operator for VQE to work with. Such an operator is given below. This was originally created by the Chemistry application module as the Hamiltonian for an H2 molecule at 0.735A interatomic distance. It's a sum of Pauli terms as below, but for now I am not going to say anything further about it since the goal is to run the algorithm, but further information on operators can be found in other tutorials." - ] - }, - { - "cell_type": "code", - "execution_count": 5, - "metadata": {}, - "outputs": [], - "source": [ - "from qiskit.aqua.operators import X, Z, I\n", - "\n", - "H2_op = (-1.052373245772859 * I ^ I) + \\\n", - " (0.39793742484318045 * I ^ Z) + \\\n", - " (-0.39793742484318045 * Z ^ I) + \\\n", - " (-0.01128010425623538 * Z ^ Z) + \\\n", - " (0.18093119978423156 * X ^ X)" - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "So let's build a VQE instance passing the operator, and a backend using a QuantumInstance, to the constructor. Now VQE does have setters so the operator and backend can also be passed later. Setting them later can be useful when creating an algorithm without this problem specific information and then later using it, say with the chemistry application module, which would create the operator for the specific chemistry problem being solved.\n", - "\n", - "Note: In order that you can run this notebook and see the exact same output the random number generator used throughout Aqua in aqua_globals, as well as the transpiler and simulator, via the QuantumInstance, are seeded. You do not have to do this but if want to be able to reproduce the exact same outcome each time then this is how it's done.\n", - "\n", - "So let's run VQE and print the result object it returns." - ] - }, - { - "cell_type": "code", - "execution_count": 6, - "metadata": { - "scrolled": true - }, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "{ 'cost_function_evals': 72,\n", - " 'eigenstate': array([-9.55448660e-05+2.12037105e-17j, 9.93766273e-01+2.25293943e-16j,\n", - " -1.11483565e-01+1.52657541e-16j, -1.77521351e-05+3.71607315e-17j]),\n", - " 'eigenvalue': (-1.857275017559769+0j),\n", - " 'optimal_parameters': { Parameter(θ[0]): 4.296520551468743,\n", - " Parameter(θ[1]): 4.426962086704216,\n", - " Parameter(θ[2]): 0.5470753710293924,\n", - " Parameter(θ[3]): 6.09294789784282,\n", - " Parameter(θ[4]): -2.598325857134344,\n", - " Parameter(θ[5]): 1.5683261371389359,\n", - " Parameter(θ[6]): -4.717618235040379,\n", - " Parameter(θ[7]): 0.3602072316165878},\n", - " 'optimal_point': array([ 4.29652055, 4.42696209, 0.54707537, 6.0929479 , -2.59832586,\n", - " 1.56832614, -4.71761824, 0.36020723]),\n", - " 'optimal_value': -1.857275017559769,\n", - " 'optimizer_evals': 72,\n", - " 'optimizer_time': 2.040440559387207}\n" - ] - } - ], - "source": [ - "from qiskit.aqua import aqua_globals\n", - "seed = 50\n", - "aqua_globals.random_seed = seed\n", - "qi = QuantumInstance(BasicAer.get_backend('statevector_simulator'), seed_transpiler=seed, seed_simulator=seed)\n", - "\n", - "ansatz = TwoLocal(rotation_blocks='ry', entanglement_blocks='cz')\n", - "slsqp = SLSQP(maxiter=1000)\n", - "vqe = VQE(operator=H2_op, var_form=ansatz, optimizer=slsqp, quantum_instance=qi)\n", - "result = vqe.run()\n", - "\n", - "import pprint\n", - "pp = pprint.PrettyPrinter(indent=4)\n", - "pp.pprint(result)" - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "From the above result we can see it took the optimizer `72` evaluations of parameter values until it found the minimum eigenvalue of `-1.85727` which is the electronic ground state energy of the given H2 molecule. The optimal parameters of the ansatz can also be seen which are the values that were in the ansatz at the minimum value." - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "## Using VQE as a MinimumEigensolver\n", - "\n", - "To close off let's use the MinimumEigensolver interface and also create a VQE instance without supplying either the operator or QuantumInstance. We later set the QuantumInstance and finally pass the operator on the `compute_minimum_eigenvalue` method (though we could have passed it in earlier via the setter instead, as long as by the time VQE runs it has an operator to work on)." - ] - }, - { - "cell_type": "code", - "execution_count": 7, - "metadata": {}, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "{ 'cost_function_evals': 72,\n", - " 'eigenstate': array([-9.55448660e-05+2.12037105e-17j, 9.93766273e-01+2.25293943e-16j,\n", - " -1.11483565e-01+1.52657541e-16j, -1.77521351e-05+3.71607315e-17j]),\n", - " 'eigenvalue': (-1.857275017559769+0j),\n", - " 'optimal_parameters': { Parameter(θ[2]): 0.5470753710293924,\n", - " Parameter(θ[7]): 0.3602072316165878,\n", - " Parameter(θ[6]): -4.717618235040379,\n", - " Parameter(θ[3]): 6.09294789784282,\n", - " Parameter(θ[5]): 1.5683261371389359,\n", - " Parameter(θ[0]): 4.296520551468743,\n", - " Parameter(θ[1]): 4.426962086704216,\n", - " Parameter(θ[4]): -2.598325857134344},\n", - " 'optimal_point': array([ 4.29652055, 4.42696209, 0.54707537, 6.0929479 , -2.59832586,\n", - " 1.56832614, -4.71761824, 0.36020723]),\n", - " 'optimal_value': -1.857275017559769,\n", - " 'optimizer_evals': 72,\n", - " 'optimizer_time': 1.5368030071258545}\n" - ] - } - ], - "source": [ - "aqua_globals.random_seed = seed\n", - "qi = QuantumInstance(BasicAer.get_backend('statevector_simulator'), seed_transpiler=seed, seed_simulator=seed)\n", - "\n", - "ansatz = TwoLocal(rotation_blocks='ry', entanglement_blocks='cz')\n", - "slsqp = SLSQP(maxiter=1000)\n", - "vqe = VQE(var_form=ansatz, optimizer=slsqp)\n", - "\n", - "vqe.quantum_instance = qi\n", - "result = vqe.compute_minimum_eigenvalue(operator=H2_op)\n", - "\n", - "pp.pprint(result)" - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "As the identical seeding was used as the prior example the result can be seen to be the same." - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "This concludes this introduction to algorithms in Qiskit. Please check out the other algorithm tutorials in this series for both broader as well as more in depth coverage of the algorithms." - ] - }, - { - "cell_type": "code", - "execution_count": 8, - "metadata": {}, - "outputs": [ - { - "data": { - "text/html": [ - "

Version Information

Qiskit SoftwareVersion
Qiskit0.23.1
Terra0.16.1
Aer0.7.1
Ignis0.5.1
Aqua0.8.1
IBM Q Provider0.11.1
System information
Python3.6.1 |Continuum Analytics, Inc.| (default, May 11 2017, 13:09:58) \n", - "[GCC 4.4.7 20120313 (Red Hat 4.4.7-1)]
OSLinux
CPUs1
Memory (Gb)5.827335357666016
Tue Nov 17 08:50:18 2020 EST
" - ], - "text/plain": [ - "" - ] - }, - "metadata": {}, - "output_type": "display_data" - }, - { - "data": { - "text/html": [ - "

This code is a part of Qiskit

© Copyright IBM 2017, 2020.

This code is licensed under the Apache License, Version 2.0. You may
obtain a copy of this license in the LICENSE.txt file in the root directory
of this source tree or at http://www.apache.org/licenses/LICENSE-2.0.

Any modifications or derivative works of this code must retain this
copyright notice, and modified files need to carry a notice indicating
that they have been altered from the originals.

" - ], - "text/plain": [ - "" - ] - }, - "metadata": {}, - "output_type": "display_data" - } - ], - "source": [ - "import qiskit.tools.jupyter\n", - "%qiskit_version_table\n", - "%qiskit_copyright" - ] - }, - { - "cell_type": "code", - "execution_count": null, - "metadata": {}, - "outputs": [], - "source": [] - } - ], - "metadata": { - "kernelspec": { - "display_name": "Python 3", - "language": "python", - "name": "python3" - }, - "language_info": { - "codemirror_mode": { - "name": "ipython", - "version": 3 - }, - "file_extension": ".py", - "mimetype": "text/x-python", - "name": "python", - "nbconvert_exporter": "python", - "pygments_lexer": "ipython3", - "version": "3.6.1" - } - }, - "nbformat": 4, - "nbformat_minor": 2 -} From 416e4ddcff26123fc98d486cff0b02a6680ae6b8 Mon Sep 17 00:00:00 2001 From: Divyanshu Singh <55018955+divshacker@users.noreply.github.com> Date: Fri, 15 Jan 2021 17:45:10 +0530 Subject: [PATCH 2/2] Fixed Some Broken Links in this tutorial --- .../01_algorithms_introduction.ipynb | 345 ++++++++++++++++++ 1 file changed, 345 insertions(+) create mode 100644 tutorials/algorithms/01_algorithms_introduction.ipynb diff --git a/tutorials/algorithms/01_algorithms_introduction.ipynb b/tutorials/algorithms/01_algorithms_introduction.ipynb new file mode 100644 index 000000000..fbe64f29f --- /dev/null +++ b/tutorials/algorithms/01_algorithms_introduction.ipynb @@ -0,0 +1,345 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# An Introduction to Algorithms in Qiskit\n", + "\n", + "This is an introduction to algorithms in Qiskit and provides a high-level overview to help understand the various aspects of the functionality to get started. Other tutorials will provide more in-depth material, on given algorithms, and ways to use them etc." + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## How is the algorithm library structured?\n", + "\n", + "Qiskit provides a number of [Algorithms](https://qiskit.org/documentation/apidoc/qiskit.aqua.algorithms.html) and they are grouped by category according to the task they can perform. For instance `Minimum Eigensolvers` to find the smallest eigen value of an operator, for example ground state energy of a chemistry Hamiltonian or a solution to an optimization problem when expressed as an Ising Hamiltonian. There are `Classifiers` for machine learning classification problems and `Amplitude Estimators` for value estimation that can be used say in financial applications. The full set of categories can be seen in the Algorithms documentation link above.\n", + "\n", + "Algorithms are configurable and often part of the configuration will be in the form of smaller building blocks, of which different instances of the building block type can be given. For instance with `VQE`, the Variational Quantum Eigensolver, it takes a trial wavefunction, in the form of a `QuantumCircuit` and a classical optimizer among other things. Such building blocks can be found as [Components](https://qiskit.org/documentation/apidoc/qiskit.aqua.components.html) and as circuits from the [Circuit Library](https://qiskit.org/documentation/apidoc/circuit_library.html).\n", + "\n", + "Let's take a look at an example to construct a VQE instance. Here `TwoLocal` is the variational form (trial wavefunction), a parameterized circuit which can be varied, and `SLSQP` a classical optimizer. These are created as separate instances and passed to VQE when it is constructed. Trying, for example, a different classical optimizer, or variational form is simply a case of creating an instance of the one you want and passing it into VQE." + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": {}, + "outputs": [], + "source": [ + "from qiskit.aqua.algorithms import VQE\n", + "from qiskit.aqua.components.optimizers import SLSQP\n", + "from qiskit.circuit.library import TwoLocal\n", + "\n", + "num_qubits = 2\n", + "ansatz = TwoLocal(num_qubits, 'ry', 'cz')\n", + "opt = SLSQP(maxiter=1000)\n", + "vqe = VQE(var_form=ansatz, optimizer=opt)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Let's draw the ansatz so we can see it's a QuantumCircuit where θ\\[0\\] through θ\\[7\\] will be the parameters that are varied as VQE optimizer finds the minimum eigenvalue. We'll come back to the parameters later in a working example below." + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAB7CAYAAACIG9xhAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAAV50lEQVR4nO3de1hUBf7H8ffMgAKCCqKSWAoilCiUl1K8oOWurD9/pmuWl9o0f1nolkputam1G6Vl2br7VFZb6rabukmmrmuWJeAFu2ipoBbeyUupqQmKIDC/PyZRAoajyJyZ6fN6nnmEc87DfPg4zHfOmcPBYrfb7YiIiBhgNTuAiIh4Dg0NERExTENDREQM09AQERHDNDRERMQwDQ0RETFMQ0NERAzT0BAREcM0NERExDANDRERMUxDQ0REDNPQEBERwzQ0RETEMA0NERExTENDREQM09AQERHDNDRERMQwDQ0RETFMQ0NERAzT0BAREcM0NERExDANDRERMUxDQ0REDNPQEBERwzQ0RETEMA0NERExzMfsAO7umzWQf9Sc+w5qBjG3mnPfZlDXrqOuXcesruuqZw2NGuQfhVMHzU7xy6CuXUddu463da3DUyIiYpiGhoiIGKahISIihmloiIiIYXoj/Cp5ZE5vdh7YiM3mi9VqIyw4ghG3TSExfqjZ0byKenYdde06ntS1hsZVNLLvNEb2nUppaQnLsl5mxoIRRIXfRHholNnRvIp6dh117Tqe0rUOT9UBm82H39xyP6VlJew5vMXsOF5LPbuOunYdd+9aQ6MOnC8pZkXWHABahkabnMZ7qWfXUdeu4+5da2hcRQs+eZZB0xoz4Al/5n04lZShbxLZIg6A6e+M4NMdK8q3fWr+IDZ985FZUatVVAL5hVBSanaS6nlDzwDnzju6Li0zO0n1vKXrwmIoOAdl6rrW3HpolJWV8eKLL9K2bVv8/PyIj48nMzOTmJgYxo4da3a8SkbcNoWlqadI+9Nxbr6+P1t3p5evS759NvM/nEZhUQHrspfQwK8RnWN+bWLaivYchb9nwOP/hmlL4I+L4d3P4YcCs5NV5sk9A+w4BK98DI+/6+h6Shos3QynC81OVpknd223w5f74aVVjsfz1PfgySWwciucLTY7XWWe0rVbD40xY8aQmprKAw88wAcffMCdd97J8OHD2bt3L506dTI7XrWCAoJJGfomn339X7JylgEQHNiMwT0m8Mqyh1nwyTM8OPAvJqe8aNM+eHm148nM/tOy86WwcRfM+gAOnzQ1XrU8rWeA9J3wRgbsvuRaROfOQ8bXjq5PuOGQBs/s+j9fwdsb4NsfLi4rKIKPcuCvH8KZIvOyOePuXbvt0Fi4cCHz589n+fLlTJ48mT59+jBlyhS6detGSUkJHTt2NDuiUw0DQhjSM4W5q56g7Kd94n5dRnHwWC6Duj9Mw4AQkxM6nDwDCzY6hoX9Z+vsQOF5mLsOyn6+0k14Ss/gePJa9qXjY3sVfZ4uhHc2ujbT5fCkrrcfgjU7HR9X9dA9ehre+8KlkS6LO3fttkNj+vTpJCUlkZiYWGF5VFQUvr6+xMU5jvXt37+fxMREoqOj6dChA+vWrTMjbpUG95zAidNHWL357fJlLZpEudUpdFm7nA8Eux2O58Ou71yX6XJ5Qs8A63LB4mS9HcdhwiOnXBToCnhM19+AxUnZdmBLnnseErzAXbt2y6Fx8OBBcnJyGDq08i+25OXlERsbS/369QF44IEHuOuuu8jNzeX1119n2LBhFBfXfMDSYrEYumVmZhjKPCs5g5F9p1ZY1sCvIUuePkG/LqMMfY2fy8zMMJzzSm9vL92IvaqXvZew2+2MnvBsnWcx0nVd9Oyqrj/5Iq/KV70/l3Tn79V1LW/b84qq3Ju7VJkdOiUO8dquL7dno9x2aACEhYVVWF5YWEhmZmb5oanjx4+zfv16xowZA0BCQgItWrQgPT0dMcZitRl4wNixWm0uyePNLAY7VNdXgdXYU5vR/xO5yC2HRmhoKAC5ubkVls+cOZMjR46Uvwmel5dH8+bNy/c6ACIiIjhw4ECN92G32w3dEhN7X71vDHh02HzaR/QwtG1iYm/DOa/0dkfSzTXmsFisvPz843We5Wp2fTk9g2u67tYh3OnhqQvee/uv6rqWt4hmvoa63rD6Xa/t+nJ7NsotLyMSGRlJXFwc06dPJyQkhPDwcNLS0li5ciWAW5855Wm6R8P6XdWvtwAN6kOHa10WyWv1iIZsJ3+MxwI0bwQRTV0WyWv1iHac4FEdiwXaNoemQa7L5C3cck/DarWyePFiYmNjSU5OZvTo0YSGhjJ+/HhsNlv5m+DXXXcd33//PUVFF8+d27dvH61atTIruse5pjEkdXB8/PNXZhYcP1x3dwebWz5SPEt0GCS0rXqdxQK+PjCym/M3cMWYzq2hQ8uq11mAgHowtOadbKmCW+5pAERHR1d6b+Kee+6hXbt2+Pv7A47DWN27d+ett95i3LhxZGVlcejQIfr06WNGZI+VFAfBDWB1Dhy/5PcEIpvB/8Q7/pXas1hgaBdo1hDSd8CPl5y5c8M1MOBGaBFsWjyvYrXCqJ6wervjTKoLv5NhtUDctY6uQ7WXcUXcdmhUZdOmTXTt2rXCstdee41Ro0Yxe/Zs6tWrx8KFC6lXr55JCT3XLW3g5kiYtMDx+dSB+qGqCxYL9L4eekVDykLHsj8NhsYB5ubyRjarYy+6bzuYvMix7M+DIcjf3FyezmOGRkFBAbm5uYwbN67C8sjISNauXWtSKoe3Vv6R7fs3ENu6Oy2bxrAofQYTh7xBfJtE3s14gazty2ge3Io/3DWf8yVFPPpGX8KbRPH4iH+ZmvvnLj0s4q4Do7qumwVfx8xFv8OChdBGLXls+D+xWW1MnTuAgsJTzB6/3uzoFVx6co+7Dgxnj2uAddlLmLNsAgumfkthUYHbPq59LjlByl0HhrOuN+euZtGaGZTZy3jgf2dxbdMYU7v2mCPVgYGBlJaW8tBDD5kdpYJ93+Vw5txpXhq3ltNnf+Bc8RmGJv6B+DaJnCw4ypY96cwev56Ia+LYkLMU//qBTBm5yOzYHslZ14F+jXlm9ApeGreWsJAIPv/acdLEM/etqOGrSlWcdX3Bum1pNG3sOENCj+sr56zrovOF/PfT13lu7GpmJWcQ3bKT6V17zNBwVzn71tM52nHhsI5tf1XhHPvcbzcRH9n7p3V92XnAja8R4QGcdR0UEEwD/0YA+Nh8sVp0/n1tOOsa4LOdK+nYti8Wi55CastZ1zsObMRisfLEm7/huYX3UFh8xqyY5fQ/Xkv5Z0/wj4+e4pE5vVnwybPknz1Rvu7MuVME+DUEoIFfIwrOnTIppXdw1vUFx388zObc1eU/hHJlaup69eZ/cFvHu01K512cdX0y/3tO5B9h+v99QGyrBP678XUTkzp4zHsa7iooIIR7+z1NQuxAPt2xgmM/XjwRv4FfI46dcnx+9txpAv0am5TSOzjrGqC4pIgX/n0vKUP/js2mh3ZtOOv6q91raNeqG74+OuHkaqjpOaR96x7YrDZujLqVxZkvmpjUQXsatdQ+ogfZex1vxG/dk0FZ2cW/XhR9bRe27c0E4MtdH3NDq65Vfg0xxlnXALPTxjIwYTytmrczI55Xcdb1/u9y2Lh9OX/8exIHvt/OvFVTq/syYoCzrmOu7ULeUcflevcc3kJYSIQpGS+loVFLEWHt8bH58sic3vjYfPGr16B8XXBgMzpE9mLiKz3Yc3gLCbGDzAvqBZx1vWP/RtbnLGHJutk8Mqc367PfNzGp53PW9eAeD/PCg2uYcf8qWjWPZXTSMyYm9XzOum4c2JS4yERSXu3Fh1/MY0C3B01M6qB9+KtgTP8Z5R+v3ZbGovTnCA9tS3ybRIb1eYxhfR4rX19YVMBzC+8m5touZkT1eM66Xv5MfqXtp84dQEjDa1wZ0Ws46/qCC6cy63FdO866HtJrEkN6TSpfb3bXFvvlXKnqF2jTIjjl5HpBdalxS+g8zLX3OfEdx7+zR7r2fkFdu5K6dh2zuq6rnrWnUYMgEy+hYeZ9m0Fdu466dh2zvt+6ul8NjRrE3Gp2gl8Ode066tp1vK1rvREuIiKGaWiIiIhhGhoiImKYhoaIiBimoSEiIoZpaIiIiGEaGiIiYpiGhoiIGKahISIihmloiIiIYRoaIiJimIaGiIgYpqEhIiKG6Sq3NfhmDeQfNee+g5p53xUynVHXrqOuXcesruuqZw2NGuQfNe+P1fzSqGvXUdeu421d6/CUiIgYpqEhIiKGaWiIiIhhek/jKnlkTm92HtiIzeaL1WojLDiCEbdNITF+qNnRvIp6dh117Tqe1LWGxlU0su80RvadSmlpCcuyXmbGghFEhd9EeGiU2dG8inp2HXXtOp7StQ5P1QGbzYff3HI/pWUl7Dm8xew4Xks9u466dh1371pDow6cLylmRdYcAFqGRpucxnupZ9dR167j7l1raFxFCz55lkHTGjPgCX/mfTiVlKFvEtkiDoDp74zg0x0ryrd9av4gNn3zkVlRPZp6dh117Tqe0rVbD42ysjJefPFF2rZti5+fH/Hx8WRmZhITE8PYsWPNjlfJiNumsDT1FGl/Os7N1/dn6+708nXJt89m/ofTKCwqYF32Ehr4NaJzzK9NTFuR3Q67v7/4+YotcPS0aXGc8uSeAcrssOPQxc9XZcPJM+blccbTuy4tg615Fz//eDvknzMvjzOe0rVbvxE+ZswYlixZwrRp0+jUqRNZWVkMHz6cY8eOkZKSYna8agUFBJMy9E3ufa4NWTnLSGh/O8GBzRjcYwKvLHuYPYe38PzYj82OWS7/HLyZAQd+uLjs4+2OW7couKML2Nzw5YWn9QzwQwG8nl5xIK/aBh9ug76x0D8eLBbz8lXHE7s+fBLeyIBTZy8uW7EFVm6FgR2h9/VmJXPO3bt2w6cCh4ULFzJ//nyWL1/O5MmT6dOnD1OmTKFbt26UlJTQsWNHsyM61TAghCE9U5i76gnKysoA6NdlFAeP5TKo+8M0DAgxOaFDaRm8tgbyfqh6/cbdsPRL12a6HJ7SM0BhMbz8MRyrYg/ODqzeDp/scHkswzyp61NnHV3/eLbyujI7LN0Mn+1xfS6j3Llrtx0a06dPJykpicTExArLo6Ki8PX1JS7OcazvySefJDo6GqvVSlpamhlRqzW45wROnD7C6s1vly9r0STKrU6hy/4WDp10PGlVZ31u1T987sITegb4fK/jMJSzrj/KgaISl0W6bJ7S9bpv4Gyx864/2AY/PR+7JXft2i0PTx08eJCcnBwmTZpUaV1eXh6xsbHUr18fgKSkJEaNGsV9993n6pgVzErOqLSsgV9Dljx9wvVhLsPne8GC8x8uux2+PAB9bnBVqup5as9grOviEsj5FjpFuCpV9Ty568/21rzNqbOw5yi0Dav7PDXxpK7dck/j4EHHJSHDwir+bxYWFpKZmVnh0FRCQgKRkZGXfR8Wi8XQLTMzo1bfS21kZmYYznmlt/QNm50+iQHYy8p4MvWFOs/i7V3n7jtSY9cA949PUde1vOWfLTWUZcBvh3tt15fbs1FuuacRGhoKQG5uLv379y9fPnPmTI4cOUKnTp3MilZrjw6bb3aECgpPH6WsrBSr1VbtNharlcL84y5MVXvu1jNA4eljBDRsjsXq/LWauq69c2dO4h8UWvN26vqyueXQiIyMJC4ujunTpxMSEkJ4eDhpaWmsXLkS4KoMDbvdyGs+2LTIvGvhJyb2xj7HWM4rtWkf/CvL+TYW4KOFzxMS+HzdZvHyrtfsgOVfOd/G1wZffvI2/vXedr5hLXl71+9vhsyvnW8T5Ad7t66u8zMDzeq6rnp2y8NTVquVxYsXExsbS3JyMqNHjyY0NJTx48djs9nK3wSX2rvxOmga5BgM1bk5EkICXRbJa93SxvFE5exIQJ8bwL+e6zJ5q14xUN/X+eO6Xwf3PJXc3bltZdHR0aSnp3PmzBny8vJITU0lOzubdu3a4e/vb3Y8r+Fjg/F9oXkjx+cWi+MH7cIP202tYOjNZqXzLg3qO7pu9NPD90LXF/SIhiS9HroqmgRC8q0Q4DhfxtH1JWUnxUH3tuZk83RueXiqOps2baJr164Vlk2bNo158+Zx7NgxsrOzmThxIpmZmbRp08aklJ6ncQA82h92HIYtB6DwPAQHOF4ZX9vE7HTeJawRTB0I276F7INwvtTxBNctyrFOrp7WofDUIMeZfzsPO34nKayRo+sm2nO+Yh4zNAoKCsjNzWXcuHEVlqemppKammpSKoe3Vv6R7fs3ENu6Oy2bxrAofQYTh7xBbOsEUl7txb7vsnlt0hbCQ6MoLCrg0Tf6Et4kisdH/MvU3JeyWqF9S8fNnVXXdeuwWJ6cNxCbzZcGfo2Yeve/KSsrdcuufWzQsbXj5s6q6zq+TSK3T2tEVIubAHjq3iU0DAhh6twBFBSeYvb49SYnv6ieD3Rt47i5s+q6buDXiDnLJwJw9OQBBvecwG97TjS1a7c9PPVzgYGBlJaW8tBDD5kdpYJ93+Vw5txpXhq3ltNnf+Bc8RmGJv6B+DaJ2Kw+/HnUUnp2uKN8e//6gUwZucjExJ7LWdeB/sH8Zdx6XkrOJDq8E5/uWKGua8FZ1wARYR2YlZzBrOSM8t9Ofua+Fc6+pFTDWddR4TeW9xxxTRy33DAAMLdrjxka7ipn33o6RzsuHNax7a8qnLpqsVgIDmpuVjSv46xrm9WG9adTWUvtpYSH6oB1bTjrGiDv6E4mvdqTN1c+bvhMRKlaTV0DFBaf4WT+d6b/Njh40OEpd5V/9gQrNr7Ge+v+QkHhKRLj76RxYDOzY3mlmrr+Ou9z/vb+OOr5+DG01yMmJvV8NXU9/7FdBPkH89f3HmTjjv+QEDvQxLSezchzyBdff0DnmCSTElakoVFLQQEh3NvvaRJiB/LpjhUc+9Gkk99/AWrq+vrrbubVCZtYnDmLVV/MZUivypehEWNq6vrCIamE9oPYfegrDY1aMPIcsiHnfe7s/agJ6SrT4alaah/Rg+y9awHYuieDsjJjly+Qy+es6/MlxeUfN/BrSD1fnZZdG866Liw+Q+lPn2/fv4EWTdz8XWY3V9NzSEnpefKO7qRNi3gz4lWioVFLEWHt8bH58sic3vjYfPGr16DC+tR/3snmXR8xc9G9ZOUsMymld3DW9Z7DW0iZk8jk1/rwxder+FWn35mY1PM56/rQsV38/m9dSHm1F8dOfUvPuDucfCWpSU3PIV/tXsONbW41KV1lOjx1FYzpP6P847Xb0liU/hzhoW2Jb5PItHverbBtYVEBzy28m5hru7g6pldw1vVLyZkVtlXXteOs6zkTK/+RlalzBxDS8BpXRvQazrruEtOPLjH9KmxvZtcWu059cMrMa/Q0bgmdh5lz32ZQ166jrl3HrK7rqmftadQgyMQTocy8bzOoa9dR165j1vdbV/erPQ0RETFMb4SLiIhhGhoiImKYhoaIiBimoSEiIoZpaIiIiGEaGiIiYpiGhoiIGKahISIihmloiIiIYRoaIiJimIaGiIgYpqEhIiKGaWiIiIhhGhoiImKYhoaIiBimoSEiIoZpaIiIiGEaGiIiYtj/AztR363q98R/AAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "execution_count": 2, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "ansatz.draw()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "But more is needed before we can run the algorithm so let's get to that next." + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## How to run an algorithm?\n", + "\n", + "In order to run an algorithm we need to have backend, a simulator or real device, on which the circuits that comprise the algorithm can be run. So for example we can use the `statevector_simulator` from the BasicAer provider." + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": {}, + "outputs": [], + "source": [ + "from qiskit import BasicAer\n", + "\n", + "backend = BasicAer.get_backend('statevector_simulator')" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Now a backend on its own does not have information on how you might want to run the circuits etc. For example how many shots, do you want a noise model, even options around transpiling the circuits. For this Qiskit Aqua has a [QuantumInstance](https://qiskit.org/documentation/stubs/qiskit.aqua.QuantumInstance.html) which is provided both the backend as well as various settings around the circuit processing and execution so for instance:" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": {}, + "outputs": [], + "source": [ + "from qiskit.aqua import QuantumInstance\n", + "\n", + "backend = BasicAer.get_backend('qasm_simulator')\n", + "quantum_instance = QuantumInstance(backend=backend, shots=800, seed_simulator=99) " + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Note: if you provide the backend directly then internally a QuantumInstance will be created from it, with default settings, so at all times the algorithms are working through a QuantumInstance.\n", + "\n", + "So now we would be able to call the [run()](https://qiskit.org/documentation/stubs/qiskit.aqua.algorithms.VQE.html#qiskit.aqua.algorithms.VQE.run) method, which is common to all algorithms and returns a result specific for the algorithm. In this case since VQE is a MinimumEigensolver we could use the [compute_mininum_eigenvalue()](https://qiskit.org/documentation/stubs/qiskit.aqua.algorithms.VQE.html#qiskit.aqua.algorithms.VQE.compute_minimum_eigenvalue) method. The latter is the interface of choice for the application modules, such as Chemistry and Optimization, in order that they can work interchangeably with any algorithm within the specific category." + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## A complete working example\n", + "\n", + "Let's put what we have learned from above together and create a complete working example. VQE will find the minimum eigenvalue, i.e. minimum energy value of a Hamilitonian operator and hence we need such an operator for VQE to work with. Such an operator is given below. This was originally created by the Chemistry application module as the Hamiltonian for an H2 molecule at 0.735A interatomic distance. It's a sum of Pauli terms as below, but for now I am not going to say anything further about it since the goal is to run the algorithm, but further information on operators can be found in other tutorials." + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": {}, + "outputs": [], + "source": [ + "from qiskit.aqua.operators import X, Z, I\n", + "\n", + "H2_op = (-1.052373245772859 * I ^ I) + \\\n", + " (0.39793742484318045 * I ^ Z) + \\\n", + " (-0.39793742484318045 * Z ^ I) + \\\n", + " (-0.01128010425623538 * Z ^ Z) + \\\n", + " (0.18093119978423156 * X ^ X)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "So let's build a VQE instance passing the operator, and a backend using a QuantumInstance, to the constructor. Now VQE does have setters so the operator and backend can also be passed later. Setting them later can be useful when creating an algorithm without this problem specific information and then later using it, say with the chemistry application module, which would create the operator for the specific chemistry problem being solved.\n", + "\n", + "Note: In order that you can run this notebook and see the exact same output the random number generator used throughout Aqua in aqua_globals, as well as the transpiler and simulator, via the QuantumInstance, are seeded. You do not have to do this but if want to be able to reproduce the exact same outcome each time then this is how it's done.\n", + "\n", + "So let's run VQE and print the result object it returns." + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": { + "scrolled": true + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "{ 'cost_function_evals': 72,\n", + " 'eigenstate': array([-9.55448660e-05+2.12037105e-17j, 9.93766273e-01+2.25293943e-16j,\n", + " -1.11483565e-01+1.52657541e-16j, -1.77521351e-05+3.71607315e-17j]),\n", + " 'eigenvalue': (-1.857275017559769+0j),\n", + " 'optimal_parameters': { Parameter(θ[0]): 4.296520551468743,\n", + " Parameter(θ[1]): 4.426962086704216,\n", + " Parameter(θ[2]): 0.5470753710293924,\n", + " Parameter(θ[7]): 0.3602072316165878,\n", + " Parameter(θ[6]): -4.717618235040379,\n", + " Parameter(θ[5]): 1.5683261371389359,\n", + " Parameter(θ[3]): 6.09294789784282,\n", + " Parameter(θ[4]): -2.598325857134344},\n", + " 'optimal_point': array([ 4.29652055, 4.42696209, 0.54707537, 6.0929479 , -2.59832586,\n", + " 1.56832614, -4.71761824, 0.36020723]),\n", + " 'optimal_value': -1.857275017559769,\n", + " 'optimizer_evals': 72,\n", + " 'optimizer_time': 1.310880184173584}\n" + ] + } + ], + "source": [ + "from qiskit.aqua import aqua_globals\n", + "seed = 50\n", + "aqua_globals.random_seed = seed\n", + "qi = QuantumInstance(BasicAer.get_backend('statevector_simulator'), seed_transpiler=seed, seed_simulator=seed)\n", + "\n", + "ansatz = TwoLocal(rotation_blocks='ry', entanglement_blocks='cz')\n", + "slsqp = SLSQP(maxiter=1000)\n", + "vqe = VQE(operator=H2_op, var_form=ansatz, optimizer=slsqp, quantum_instance=qi)\n", + "result = vqe.run()\n", + "\n", + "import pprint\n", + "pp = pprint.PrettyPrinter(indent=4)\n", + "pp.pprint(result)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "From the above result we can see it took the optimizer `72` evaluations of parameter values until it found the minimum eigenvalue of `-1.85727` which is the electronic ground state energy of the given H2 molecule. The optimal parameters of the ansatz can also be seen which are the values that were in the ansatz at the minimum value." + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Using VQE as a MinimumEigensolver\n", + "\n", + "To close off let's use the MinimumEigensolver interface and also create a VQE instance without supplying either the operator or QuantumInstance. We later set the QuantumInstance and finally pass the operator on the `compute_minimum_eigenvalue` method (though we could have passed it in earlier via the setter instead, as long as by the time VQE runs it has an operator to work on)." + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "{ 'cost_function_evals': 72,\n", + " 'eigenstate': array([-9.55448660e-05+2.12037105e-17j, 9.93766273e-01+2.25293943e-16j,\n", + " -1.11483565e-01+1.52657541e-16j, -1.77521351e-05+3.71607315e-17j]),\n", + " 'eigenvalue': (-1.857275017559769+0j),\n", + " 'optimal_parameters': { Parameter(θ[4]): -2.598325857134344,\n", + " Parameter(θ[5]): 1.5683261371389359,\n", + " Parameter(θ[6]): -4.717618235040379,\n", + " Parameter(θ[7]): 0.3602072316165878,\n", + " Parameter(θ[1]): 4.426962086704216,\n", + " Parameter(θ[0]): 4.296520551468743,\n", + " Parameter(θ[2]): 0.5470753710293924,\n", + " Parameter(θ[3]): 6.09294789784282},\n", + " 'optimal_point': array([ 4.29652055, 4.42696209, 0.54707537, 6.0929479 , -2.59832586,\n", + " 1.56832614, -4.71761824, 0.36020723]),\n", + " 'optimal_value': -1.857275017559769,\n", + " 'optimizer_evals': 72,\n", + " 'optimizer_time': 2.8010470867156982}\n" + ] + } + ], + "source": [ + "aqua_globals.random_seed = seed\n", + "qi = QuantumInstance(BasicAer.get_backend('statevector_simulator'), seed_transpiler=seed, seed_simulator=seed)\n", + "\n", + "ansatz = TwoLocal(rotation_blocks='ry', entanglement_blocks='cz')\n", + "slsqp = SLSQP(maxiter=1000)\n", + "vqe = VQE(var_form=ansatz, optimizer=slsqp)\n", + "\n", + "vqe.quantum_instance = qi\n", + "result = vqe.compute_minimum_eigenvalue(operator=H2_op)\n", + "\n", + "pp.pprint(result)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "As the identical seeding was used as the prior example the result can be seen to be the same." + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "This concludes this introduction to algorithms in Qiskit. Please check out the other algorithm tutorials in this series for both broader as well as more in depth coverage of the algorithms." + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "

Version Information

Qiskit SoftwareVersion
Qiskit0.23.1
Terra0.16.1
Aer0.7.1
Ignis0.5.1
Aqua0.8.1
IBM Q Provider0.11.1
System information
Python3.7.8 | packaged by conda-forge | (default, Jul 31 2020, 02:25:08) \n", + "[GCC 7.5.0]
OSLinux
CPUs8
Memory (Gb)31.40900421142578
Fri Jan 15 12:11:24 2021 UTC
" + ], + "text/plain": [ + "" + ] + }, + "metadata": {}, + "output_type": "display_data" + }, + { + "data": { + "text/html": [ + "

This code is a part of Qiskit

© Copyright IBM 2017, 2021.

This code is licensed under the Apache License, Version 2.0. You may
obtain a copy of this license in the LICENSE.txt file in the root directory
of this source tree or at http://www.apache.org/licenses/LICENSE-2.0.

Any modifications or derivative works of this code must retain this
copyright notice, and modified files need to carry a notice indicating
that they have been altered from the originals.

" + ], + "text/plain": [ + "" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "import qiskit.tools.jupyter\n", + "%qiskit_version_table\n", + "%qiskit_copyright" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3", + "language": "python", + "name": "python3" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 3 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython3", + "version": "3.7.8" + } + }, + "nbformat": 4, + "nbformat_minor": 2 +}