In this project we develop a quantum algorithm to realize finite temperature simulation on a quantum computer. As quantum computers use real-time evolution we did not use the imaginary time methods popular on classical algorithms. Instead, we implemented a real-time therom field dynamics formalism, which has the added benefit of being able to compute quantities that are both time-and temperature-dependent. To implement thermo field dynamics we apply a unitary transformation [1] to discrete quantum mechanical operators to make new Hamiltoni-ans with encoded temperature dependence. The method works normally for fermions, which have a finite representation, but needs some modification to work with bosons. These Hamilto-nians are then processed into a Pauli matrix representation in order to be used as input for IBM's Qiskit package [2]. We then use IBM's quantum simulator to calculate an approximation to the Hamiltonaian's ground state energy via the variational quantum eigensolver (VQE) algorithm [3]. This approximation is then compared to a classically calculated value for the exact energy. The thermo field dynamics quantum algorithm has general applications to material science, high-energy physics and nuclear physics, particularly in those situations involving realtime evolution at high temperature.
In this project we develop a quantum algorithm to realize finite temperature simulation on a quantum computer. As quantum computers use real-time evolution we did not use the imaginary time methods popular on classical algorithms. Instead, we implemented a real-time therom field dynamics formalism, which has the added benefit of being able to compute quantities that are both time-and temperature-dependent. To implement thermo field dynamics we apply a unitary transformation [1] to discrete quantum mechanical operators to make new Hamiltoni-ans with encoded temperature dependence. The method works normally for fermions, which have a finite representation, but needs some modification to work with bosons. These Hamilto-nians are then processed into a Pauli matrix representation in order to be used as input for IBM's Qiskit package [2]. We then use IBM's quantum simulator to calculate an approximation to the Hamiltonaian's ground state energy via the variational quantum eigensolver (VQE) algorithm [3]. This approximation is then compared to a classically calculated value for the exact energy. The thermo field dynamics quantum algorithm has general applications to material science, high-energy physics and nuclear physics, particularly in those situations involving realtime evolution at high temperature.
Quantum computers have the potential to transform the ways in which we tackle some important problems. The efforts by companies like Google, IBM and Microsoft to construct quantum computers have been making headlines for years. Equally important is the challenge of translating problems into a state that can be fed to these machines. Because quantum computers are in essence controllable quantum systems, the problems that most naturally map to them are those of quantum mechanics. Quantum chemistry has seen particular success in the form of the variational quantum eigensolver (VQE) algorithm, which is used to determine the ground state energy of molecular systems. The goal of our work has been to use the matrix formulation of quantum mechanics to translate other systems so that they can be run through this same algorithm. We describe two ways of accomplishing this using a position basis approach and a Gaussian basis approach. We use this translation to compute finite temperature quantities such as the average energy as a function of temperature. We do this by constructing a 0+1 dimensional thermal field theory, constructing a thermal double of the system and using tensor products to construct operators in the system coupled to a heat bath. We also discuss how to include a chemical potential in the quantum computations. We then connect the 0+1 formalism to gauge theory by using an effective matrix model description used in nuclear theory. We study effective potentials for components of the gauge field and use the VQE algorithm to calculate the ground state energies. We also visualize the wave functions from the eigensolver and make comparisons to theoretical results obtained with continuous operators.
In this project we studied a simplified potential from [1] describing nuclear physics from quark and gluon fields, where spatial degrees of freedom have been integrated over. This potential deals with the phase structure of the strong interaction. We created a matrix representation of the Hamiltonian for this theory using discrete quantum mechanics [2]. We then feed this Hamiltonian to IBM's Qiskit package [3]. We run on IBM's 32-qubit quantum simulator, which approximates the ground state energy using the variational quantum eigensolver (VQE) algorithm [4]. This result is then compared with the exact ground state energy, which we can calculate classically since the matrices we use are still small enough. The bosonic Hamiltonian for a 6-qubit representation took over 1000 Pauli factors to represent fully. To improve our performance, we can look for a more efficient representation of this model's Hamiltonian. We will develop a better variational trial wave function and initial state, as well as include the effect of fermions and heavy quarks on the potential.
In this thesis we have implemented three many-body methods in quantum mechanics on quantum computers, namely the Quantum Phase Estimation (QPE) algorithm, the Variational Quantum Eigensolver (VQE) algorithm and the Quantum Adiabatic Time Evolution (QATE) algorithm.We utilized these methods to approximate the ground state energy of the pairing model Hamiltonian.We provided a comparison between the ideal simulation of a quantum computer and a simulation including the noise model from some of IBM's actual quantum devices.For the VQE algortihm, we also run the method on the IBM Q London five-qubit device.As the time complexity (circuit depth) of a quantum algorithm (quantum circuit) is often the bottleneck on current quantum computers, we also introduced a circuit optimization scheme, called Recursive Circuit Optimization, which can be utilized to approximate the action of a quantum circuit by learning the parameters of a parametrized unitary operation and only evaluating small parts of the circuit at a time.We tested out the method by trying to optimize a random circuit, the QATE algorithm and the unitary coupled cluster doubles (UCCD) ansatz.Finally, we implemented a quantum circuit that allows us to encode a real data set into the amplitudes of a quantum state, and we utilized this circuit to implement several realizations of dense and recurrent neural networks on quantum computers.We proposed neural network layers that are consisting of hardware-and parameter-efficient parametrized quantum circuits, and we have made a Python framework to construct and train deep quantum neural networks.We then tested the neural networks by trying to learn a simple non-linear function as well as finding the smallest eigenvalue of the pairing Hamiltonian Full Configuration Interaction (FCI) matrix.We found that all the quantum many-body methods showed promise in providing an upper bound on the ground state energy for the ideal simulation of a quantum computer.When performing the simulation with noise and the execution on the real quantum device however, we found that only the VQE algorithm was able to yield reasonable results due to short circuit depth.Still, the success of the method was found to be dependent on the interaction strength g of the quantum system, as the approximation of the ground state energy started to deviate from the actual value as we increased g.For the recursive circuit optimization of a random circuit, we found that we were able to approximate the state given by the random circuit, and we learned that we should expect some error in the approximation due to finite measurements and the flexibility of the parametrized unitary.When utilizing the recursive optimization scheme to approximate the QATE algorithm, we found that although we were not able to exactly reproduce the action of the operator, we did provide an upper bound to the ground state energy close to the one given by the actual QATE algorithm.In the same time, the circuit depth of the algorithm was drastically decreased.When utilizing the scheme for the VQE algorithm, we got promising results while reducing the circuit depth of the UCCD ansatz.Over to the neural networks, we managed to approximate a non-linear function utilizing a hardware-efficient three layered network, with only three qubits in the intermediate layers.We showed that our proposed neural network are suitable for near-term quantum computers by successfully approximating the lowest eigenvalue of the FCI matrix when running a simulation with the noise model from the IBM Q London five-qubit computer.i I would first like to thank my supervisor, Morten Hjorth-Jensen, for being helpful and motivational throughout my time at the University of Oslo.Your door has always been open when help is needed, which has been a great relief in stressful times.I would also like to thank my friends and family for believing in me and providing emotional support.Your continuous interest in my studies and well-being has meant a lot to me.My next thanks goes to Heidi K. J.I highly appreciate your never-ending support.Finally I would like to thank the computational physics study group.You have been invaluable to me in both a social and educational manner.
The developments of quantum computing algorithms and experiments for atomic scale simulations have largely focused on quantum chemistry for molecules, while their application in condensed matter systems is scarcely explored. Here we present a quantum algorithm to perform dynamical mean field theory (DMFT) calculations for condensed matter systems on currently available quantum computers, and demonstrate it on two quantum hardware platforms. DMFT is required to properly describe the large class of materials with strongly correlated electrons. The computation-ally challenging part arises from solving the effective problem of an interacting impurity coupled to a bath, which scales exponentially with system size on conventional computers. An exponential speedup is expected on quantum computers, but the algorithms proposed so far are based on real time evolution of the wavefunction, which requires high-depth circuits and hence very low noise levels in the quantum hardware. Here we propose an alternative approach, which uses the variational quantum eigensolver (VQE) method for ground and excited states to obtain the needed quantities as part of an exact diagonalization impurity solver. We present the algorithm for a two site DMFT system, which we benchmark using simulations on conventional computers as well as experiments on superconducting and trapped ion qubits, demonstrating that this method is suitable for running DMFT calculations on currently available quantum hardware.
We ran quantum simulations on lattices of interacting fermions in a background magnetic field [1]. We created discretized Hamiltonian matrices for each system and pre-processed them so that they could be read by the quantum simulator implemented in IBM's Qiskit package [2]. The simulator uses the variational quantum eigensolver (VQE) algorithm [3] to approximate the ground state energy of the input Hamiltonian, which is then compared to the classically-calculated exact value. This classical evaluation is possible because the systems studied are still within the grasp of classical computers. We also investigated ways to simplify the system before the pre-processing step. We found that in almost all cases these simplifications made our quantum calculations converge to better answers using simpler starting conditions. We are exploring other ways to interface with IBM's quantum computing software so that we no longer need to explicitly define the Hamiltonian matrices for our systems, a problem which will itself soon become intractable.
The developments of quantum computing algorithms and experiments for atomic scale simulations have largely focused on quantum chemistry for molecules, while their application in condensed matter systems is scarcely explored. Here we present a quantum algorithm to perform dynamical mean field theory (DMFT) calculations for condensed matter systems on currently available quantum computers, and demonstrate it on two quantum hardware platforms. DMFT is required to properly describe the large class of materials with strongly correlated electrons. The computationally challenging part arises from solving the effective problem of an interacting impurity coupled to a bath, which scales exponentially with system size on conventional computers. An exponential speedup is expected on quantum computers, but the algorithms proposed so far are based on real time evolution of the wavefunction, which requires high-depth circuits and hence very low noise levels in the quantum hardware. Here we propose an alternative approach, which uses the variational quantum eigensolver (VQE) method for ground and excited states to obtain the needed quantities as part of an exact diagonalization impurity solver. We present the algorithm for a two site DMFT system, which we benchmark using simulations on conventional computers as well as experiments on superconducting and trapped ion qubits, demonstrating that this method is suitable for running DMFT calculations on currently available quantum hardware.