Hamiltonian simulation

Hamiltonian simulation (also referred to as quantum simulation) is a problem in quantum information science that attempts to find the computational complexity and quantum algorithms needed for simulating quantum systems. Hamiltonian simulation is a problem that demands algorithms which implement the evolution of a quantum state efficiently. The Hamiltonian simulation problem was proposed by Richard Feynman in 1982, where he proposed a quantum computer as a possible solution since the simulation of general Hamiltonians seem to grow exponentially with respect to the system size.^[1]

Problem statement

In the Hamiltonian simulation problem, given a Hamiltonian $H$ ( $2^{n}\times 2^{n}$ hermitian matrix acting on $n$ qubits), a time $t$ and maximum simulation error $\epsilon$ , the goal is to find an algorithm that approximates $U$ such that $||U-e^{-iHt}||\leq \epsilon$ , where $e^{-iHt}$ is the ideal evolution and $||\cdot ||$ is the spectral norm. A special case of the Hamiltonian simulation problem is the local Hamiltonian simulation problem. This is when $H$ is a k-local Hamiltonian on $n$ qubits where $H=\sum _{j\mathop {=} 1}^{m}H_{j}$ and $H_{j}$ acts non-trivially on at most $k$ qubits instead of $n$ qubits.^[2] The local Hamiltonian simulation problem is important because most Hamiltonians that occur in nature are k-local.^[2]

Techniques

Product formulas

Also known as Trotter formulas or Trotter–Suzuki decompositions, Product formulas simulate the sum-of-terms of a Hamiltonian by simulating each one separately for a small time slice.^[3]^[4] If $H=A+B+C$ , then $U=e^{-i(A+B+C)t}=(e^{-iAt/r}e^{-iBt/r}e^{-iCt/r})^{r}$ for a large $r$ ; where $r$ is the number of time steps to simulate for. The larger the $r$ , the more accurate the simulation.

If the Hamiltonian is represented as a Sparse matrix, the distributed edge coloring algorithm can be used to decompose it into a sum of terms; which can then be simulated by a Trotter–Suzuki algorithm.^[5]

Taylor series

$e^{-iHt}=\sum _{n\mathop {=} 0}^{\infty }{\frac {(-iHt)^{n}}{n!}}=I-iHt-{\frac {H^{2}t^{2}}{2}}+{\frac {iH^{3}t^{3}}{6}}+\cdots$ by the Taylor series expansion.^[6] This says that during the evolution of a quantum state, the Hamiltonian is applied over and over again to the system with a various number of repetitions. The first term is the identity matrix so the system doesn't change when it is applied, but in the second term the Hamiltonian is applied once. For practical implementations, the series has to be truncated $\left(\sum _{n\mathop {=} 0}^{N}{\frac {(-iHt)^{n}}{n!}}\right)$ , where the bigger the $N$ , the more accurate the simulation.^[7] This truncated expansion is then implemented via the linear combination of unitaries (LCU) technique for Hamiltonian simulation.^[6] Namely, one decomposes the Hamiltonian $H=\sum _{\ell =1}^{L}\alpha _{\ell }H_{\ell }$ such that each $H_{\ell }$ is unitary (for instance, the Pauli operators always provide such a basis), and so each $H^{n}=\sum _{\ell _{1},\ldots ,\ell _{n}=1}^{L}\alpha _{\ell _{1}}\cdots \alpha _{\ell _{n}}H_{\ell _{1}}\cdots H_{\ell _{n}}$ is also a linear combination of unitaries.

Quantum walk

In the quantum walk, a unitary operation whose spectrum is related to the Hamiltonian is implemented then the Quantum phase estimation algorithm is used to adjust the eigenvalues. This makes it unnecessary to decompose the Hamiltonian into a sum-of-terms like the Trotter-Suzuki methods.^[6]

Quantum signal processing

The quantum signal processing algorithm works by transducing the eigenvalues of the Hamiltonian into an ancilla qubit, transforming the eigenvalues with single qubit rotations and finally projecting the ancilla.^[8] It has been proved to be optimal in query complexity when it comes to Hamiltonian simulation.^[8]

Complexity

The table of the complexities of the Hamiltonian simulation algorithms mentioned above. The Hamiltonian simulation can be studied in two ways. This depends on how the Hamiltonian is given. If it is given explicitly, then gate complexity matters more than query complexity. If the Hamiltonian is described as an Oracle (black box) then the number of queries to the oracle is more important than the gate count of the circuit. The following table shows the gate and query complexity of the previously mentioned techniques.

Technique	Gate complexity	Query complexity
Product formula 1st order	$O\left({\frac {t^{2}}{\epsilon }}\right)$ ^[7]	$O\left(d^{3}t\left({\frac {dt}{\epsilon }}\right)^{{\frac {1}{2}}k}\right)$ ^[9]
Taylor series	$O\left({\frac {t\log ^{2}\left({\frac {t}{\epsilon }}\right)}{\log \log {\frac {t}{\epsilon }}}}\right)$ ^[7]	$O\left({\frac {d^{2}\|\|H\|\|_{\rm {max}}\log {\frac {d^{2}\|\|H\|\|_{max}}{\epsilon }}}{\log \log {\frac {d^{2}\|\|H\|\|_{\rm {max}}}{\epsilon }}}}\right)$ ^[6]
Quantum walk	$O\left({\frac {t}{\sqrt {\epsilon }}}\right)$ ^[7]	$O\left(d\|\|H\|\|_{\rm {max}}{\frac {t}{\sqrt {\epsilon }}}\right)$ ^[5]
Quantum signal processing	$O\left(t+\log {\frac {1}{\epsilon }}\right)$ ^[7]	$O\left(td\|\|H\|\|_{\rm {max}}+{\frac {\log {\frac {1}{\epsilon }}}{\log \log {\frac {1}{\epsilon }}}}\right)$ ^[8]

Where $||H||_{\rm {max}}$ is the largest entry of $H$ .

References

↑ Richard P Feynman (1982). "Simulating physics with computers". International Journal of Theoretical Physics. 21 (6): 467–488. Bibcode:1982IJTP...21..467F. doi:10.1007/BF02650179. S2CID 124545445. Retrieved 2019-05-04.
1 2 Lloyd, S. (1996). "Universal quantum simulators". Science. 273 (5278): 1073–8. Bibcode:1996Sci...273.1073L. doi:10.1126/science.273.5278.1073. PMID 8688088. S2CID 43496899.
↑ Suzuki, Masuo (1991). "General theory of fractal path integrals with applications to many‐body theories and statistical physics". Journal of Mathematical Physics. 32 (2): 400–407. Bibcode:1991JMP....32..400S. doi:10.1063/1.529425.
↑ Berry, Dominic; Ahokas, Graeme; Cleve, Richard; Sanders, Barry (2007). "Efficient Quantum Algorithms for Simulating Sparse Hamiltonians". Communications in Mathematical Physics. 270 (2): 359–371. arXiv:quant-ph/0508139. Bibcode:2007CMaPh.270..359B. doi:10.1007/s00220-006-0150-x. S2CID 37923044.
1 2 Berry, Dominic; Childs, Andrew; Kothari, Robin (2015). "Hamiltonian simulation with nearly optimal dependence on all parameters". 2015 IEEE 56th Annual Symposium on Foundations of Computer Science. pp. 792–809. arXiv:1501.01715. Bibcode:2015arXiv150101715B. doi:10.1109/FOCS.2015.54. ISBN 978-1-4673-8191-8. S2CID 929117.
1 2 3 4 Berry, Dominic; Childs, Andrew; Cleve, Richard; Kothari, Robin; Somma, Rolando (2015). "Simulating Hamiltonian dynamics with a truncated Taylor series". Physical Review Letters. 114 (9): 090502. arXiv:1412.4687. Bibcode:2015PhRvL.114i0502B. doi:10.1103/PhysRevLett.114.090502. PMID 25793789. S2CID 15682119.
1 2 3 4 5 Childs, Andrew; Maslov, Dmitri; Nam, Yunseong (2017). "Toward the first quantum simulation with quantum speedup". Proceedings of the National Academy of Sciences. 115 (38): 9456–9461. arXiv:1711.10980. Bibcode:2018PNAS..115.9456C. doi:10.1073/pnas.1801723115. PMC 6156649. PMID 30190433.
1 2 3 Low, Guang Hao; Chuang, Isaac (2017). "Optimal Hamiltonian Simulation by Quantum Signal Processing". Physical Review Letters. 118 (1): 010501. arXiv:1606.02685. Bibcode:2017PhRvL.118a0501L. doi:10.1103/PhysRevLett.118.010501. PMID 28106413. S2CID 1118993.
↑ Kothari, Robin (Dec 8, 2017). Quantum algorithms for Hamiltonian simulation: Recent results and open problems (Youtube). United States: IBM Research.

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.

[Feynman1982-1] Richard P Feynman (1982). "Simulating physics with computers". International Journal of Theoretical Physics. 21 (6): 467–488. Bibcode:1982IJTP...21..467F. doi:10.1007/BF02650179. S2CID 124545445. Retrieved 2019-05-04.

[Lloyd1996-2] 1 2 Lloyd, S. (1996). "Universal quantum simulators". Science. 273 (5278): 1073–8. Bibcode:1996Sci...273.1073L. doi:10.1126/science.273.5278.1073. PMID 8688088. S2CID 43496899.

[3] Suzuki, Masuo (1991). "General theory of fractal path integrals with applications to many‐body theories and statistical physics". Journal of Mathematical Physics. 32 (2): 400–407. Bibcode:1991JMP....32..400S. doi:10.1063/1.529425.

[4] Berry, Dominic; Ahokas, Graeme; Cleve, Richard; Sanders, Barry (2007). "Efficient Quantum Algorithms for Simulating Sparse Hamiltonians". Communications in Mathematical Physics. 270 (2): 359–371. arXiv:quant-ph/0508139. Bibcode:2007CMaPh.270..359B. doi:10.1007/s00220-006-0150-x. S2CID 37923044.

[walk-5] 1 2 Berry, Dominic; Childs, Andrew; Kothari, Robin (2015). "Hamiltonian simulation with nearly optimal dependence on all parameters". 2015 IEEE 56th Annual Symposium on Foundations of Computer Science. pp. 792–809. arXiv:1501.01715. Bibcode:2015arXiv150101715B. doi:10.1109/FOCS.2015.54. ISBN 978-1-4673-8191-8. S2CID 929117.

[taylor-6] 1 2 3 4 Berry, Dominic; Childs, Andrew; Cleve, Richard; Kothari, Robin; Somma, Rolando (2015). "Simulating Hamiltonian dynamics with a truncated Taylor series". Physical Review Letters. 114 (9): 090502. arXiv:1412.4687. Bibcode:2015PhRvL.114i0502B. doi:10.1103/PhysRevLett.114.090502. PMID 25793789. S2CID 15682119.

[large-7] 1 2 3 4 5 Childs, Andrew; Maslov, Dmitri; Nam, Yunseong (2017). "Toward the first quantum simulation with quantum speedup". Proceedings of the National Academy of Sciences. 115 (38): 9456–9461. arXiv:1711.10980. Bibcode:2018PNAS..115.9456C. doi:10.1073/pnas.1801723115. PMC 6156649. PMID 30190433.

[LC2016-8] 1 2 3 Low, Guang Hao; Chuang, Isaac (2017). "Optimal Hamiltonian Simulation by Quantum Signal Processing". Physical Review Letters. 118 (1): 010501. arXiv:1606.02685. Bibcode:2017PhRvL.118a0501L. doi:10.1103/PhysRevLett.118.010501. PMID 28106413. S2CID 1118993.

[9] Kothari, Robin (Dec 8, 2017). Quantum algorithms for Hamiltonian simulation: Recent results and open problems (Youtube). United States: IBM Research.

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]