13

There are plenty of different variants, particularly with regards to the conditions on the Hamiltonian. It's a bit of a game, for example, to try and find the simplest possible class of Hamiltonians for which simulation is still BQP-complete. The statement will roughly be along the lines of: let $|\psi\rangle$ be a (normalised) product state, $H$ be a ...


9

One way order to perform Z rotations by arbitrary angles is to approximate them with a sequence of Hadamard and T gates. If you need the approximation to have maximum error $\epsilon$, there are known constructions that do this using roughly $3 \lg \frac{1}{\epsilon}$ T gates. See "Optimal ancilla-free Clifford+T approximation of z-rotations" by Ross et al. ...


7

A conventional Hamiltonian is Hermitian. Hence, if it contains a non-Hermitian term, it must either also contain its Hermitian conjuagte as another term, or have 0 weight. In this particular case, since $Z\otimes X\otimes Y$ is Hermitian itself, the coefficient would have to be 0. So, if you're talking about conventional Hamiltonians, you've probably made a ...


6

Your reasoning is correct if your two Hamiltonians commute. But, as you say, it doesn't work if they don't commute. In that case, the trick is to find something that approximates the the thing you want. So, what you should really be thinking about is taking terms in the opposite order: $$ e^{iH_1t/2}e^{iH_2t/2}e^{iH_1t/2}e^{iH_2t/2}\approx(e^{i(H_1+H_2)t/2})...


6

Your question remains very unclear as to what it actually is that you want to calculate. There is no direct correspondence between a system Hamiltonian and the quantum state of the system. No matter what the Hamiltonian, any quantum state is a valid state of the system. Where a Hamiltonian comes in useful is, if you know the state at some time (say, $t=0$),...


5

I don't know why/how the authors of that paper do what they do. However, here's how I'd go about it for this special case (and it is a very special case): You can write the Hamiltonian as a Pauli decomposition $$ A=15\mathbb{I}\otimes\mathbb{I}+9Z\otimes X+5X\otimes Z-3Y\otimes Y. $$ Update: It should be $+3Y\otimes Y$. But I don't want to redraw all my ...


5

The insight that suggests that sparse matrices are useful goes along the lines of: for any $H$, we can decompose it in terms of a set of $H_i$ whose individual components all commute (making diagonalisation straightforward), $$ H=\sum_{i=1}^mH_i. $$ If the matrix is sparse, then you shouldn't need too many distinct $H_i$. Then you can simulate the ...


5

Note: I'm deliberately leaving a few gaps here. Hopefully I'm saying enough to let you piece te rest together! Let's say that you want to implement $V$ on some state $$ \sum_{x\in\{0,1\}^n}\alpha_x|x\rangle $$ You can fairly easily write down what that state produces. Think about the Hamiltonian $Z^{\otimes n}$. What eigenvalues does it have? $\lambda=\pm 1$...


4

When using a simulator, it doesn't really matter what kind of qubit you refer to. You can even mix-and-match the types. The type of qubit only becomes relevant when you intend to run on a device, because devices have qubits at specific locations. For example, if you wanted to run on Bristlecone, you would limit yourself to GridQubit instances that actually ...


4

You want to start by being careful with the sizes of the operators. $\hat U$ acts on $q$ qubits, and $\hat H$ acts on $n<q$ qubits. I believe that $|G\rangle$ is a state of $q-n$ qubits. So, what we really need to talk about is two distinct sets of qubits. Let me call them sets $A$ and $B$. $A$ contains $n$ qubits, and $B$ contains $q-n$ qubits. I'll use ...


4

1. Definitions Names and symbols used in this answer follow the ones defined in Quantum linear systems algorithms: a primer (Dervovic, Herbster, Mountney, Severini, Usher & Wossnig, 2018). A recall is done below. 1.1 Register names Register names are defined in Figure 5. of Quantum linear systems algorithms: a primer (Dervovic, Herbster, Mountney, ...


4

It's not an area I personally know much about, but I know that many of my physicist friends are excited about being able to investigate the Hubbard Model on larger lattices than we can simulate today. There are known and published algorithms for finding the ground state energy, computing Green's functions, and other important characteristics of the model. ...


4

Several graph theory problems such as Graph Coloring (which is NP-complete) can be cleverly mapped to finding ground states of some classes of Hamiltonians. Graph Partitioning using Quantum Annealing on the D-Wave System Quantum annealing of the graph coloring problem


4

Generally speaking, a realization of a quantum gate involves coherent manipulation of a two-level system (but this is nothing new to you, maybe). For example, you can use two long-lived electronic states in a trapped atom (neutral or ionized in vacuo) and use an applied electric field to implement single-qubit operations (see trapped ions or optical lattices,...


4

Reformulating your question: How to perform Hamiltonian Simulation for a generic square matrix $A$? Quick answer: it is not possible. The goal of Hamiltonian Simulation (HS) is to find a quantum circuit (i.e. a succession of gates) that acts like $U(t) = e^{-iAt}$ on a quantum state. Here $U(t)$ needs to be unitary (because of the properties of quantum ...


3

$\newcommand{\bra}[1]{\left\langle#1\right|}\newcommand{\ket}[1]{\left|#1\right\rangle}\newcommand{\proj}[1]{|#1\rangle\langle#1|}\newcommand{\half}{\frac12}$In answer to your first question, I wrote myself some notes some time ago about my understanding of how it worked. The notation is probably a bit different (I've tried to bring it more into line, but it'...


3

In this article the authors stated that they used this Group Leader's algorithm in order to obtain the circuit implementing the hamiltonian simulation used as a subroutine in an instance of HHL algorithm. Unfortunately though, I did not understand quite well how they actually managed to find the circuit with that method.


3

The simplest method to implement $e^{iA\theta}$ for a small, Hermitian matrix $A$ is to: Find the eigenvectors $|\lambda\rangle$ and eigenvalues $\lambda$ of $A$. Construct the unitary $U=\sum_i|i\rangle\langle\lambda_i|$. Implement the gate sequence: $U$ $e^{i\theta\sum_i\lambda_i|i\rangle\langle i|}$ $U^\dagger$ Now, for one qubit, you have the middle ...


3

Yes, in this special case the circuit will simplify as you suggest. The advantage of the circuit that was given is that it generalises more easily, and works for any $H$ which has $\pm 1$ eigenvalues. Here's a general form of the circuit for your reference: This essentially comes down to an issue of how you can reversibly compute a one-bit function $f:\{0,1\...


3

Let's say you have a Hamiltonian of the form $$ H=\sigma_1\otimes\sigma_2\otimes\sigma_2\otimes\ldots\otimes\sigma_n $$ There's a straightforward circuit construction that lets you implement its time evolution $e^{-iHt}$. The trick is basically to decompose the state that you're evolving into the components that are in the $\pm 1$ eigenspaces of $H$. Then, ...


3

If you are looking for a more complete implementation of a quantum variational algorithm in the context of Cirq, I would recommend looking at the second example in the OpenFermion-Cirq notebook found here. It uses a custom ansatz for hydrogen in a minimal basis, but makes a bit more explicit all the required pieces. Another good example, perhaps without ...


2

Hamiltonian simulation is usually talking about simulating a quantum system on another quantum system. In particular, a digital simulation would involve implementing the simulation on a quantum computer using, e.g. the gate model. In this setting, there are several different strategies, i.e. algorithms, which one can use to perform the simulation. These ...


2

Typically the second option: you map the fermionic occupation number for each single-particle state to a qubit. (Also, the usual convention is that $0$ denotes unoccupied, $1$ denotes occupied ;) This mapping is usually accomplished via the Jordan-Wigner or Bravyi-Kitaev transformation, or some hybrid of the two: see https://arxiv.org/abs/1208.5986 for a ...


2

I don't have a completely general method for doing what you ask. However, there are a few of the steps that I might take: The $4\times 4$ matrix $H$ can always be written in the form $$ H=a\mathbb{I}\otimes\mathbb{I}+\underline{n}_1\cdot\underline{\sigma}\otimes\mathbb{I}+\mathbb{I}\otimes\underline{n}_2\cdot\underline{\sigma}+\underline{\sigma}\cdot M\cdot\...


2

I consider the leading candidate (considering you explicitly excluded quantum chemistry and hence all of biochemistry) the calculation of nuclear properties: This will, with suitable quantum computers, one day allow to compare experimental, atomic physics with theoretical (quantum computer calculated) values from the standard model. It's not entirely science ...


2

Update on the subject: there are several implementations in the wild. I don't know if you still need them, but even if you don't it will hopefully be useful to other people. I chose to list the implementations by "provenance" rather than by the algorithm used because there are not that much implementations. This may change in the future. Qiskit-aqua: ...


Only top voted, non community-wiki answers of a minimum length are eligible