Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1806.02793 · Physical Review A

Quantum simulation of the Sachdev-Ye-Kitaev model by asymmetric qubitization

Ryan Babbush, Dominic William Berry, Hartmut Neven

Published Apr 4, 2019• 70 citations
arXivHTMLPDF

Abstract

We show that one can quantum simulate the dynamics of a Sachdev-Ye-Kitaev model with $N$ Majorana modes for time $t$ to precision $\ensuremath{\epsilon}$ with gate complexity $O({N}^{7/2}t+{N}^{5/2}t\phantom{\rule{0.16em}{0ex}}\mathrm{polylog}(N/\ensuremath{\epsilon}))$. In addition to scaling sublinearly in the number of Hamiltonian terms, this gate complexity represents an exponential improvement in $1/\ensuremath{\epsilon}$ and large polynomial improvement in $N$ and $t$ over prior state-of-the-art algorithms which scale as $O({N}^{10}{t}^{2}/\ensuremath{\epsilon})$. Our approach involves a variant of the qubitization technique in which we encode the Hamiltonian $H$ as an asymmetric projection of a signal oracle $U$ onto two different signal states prepared by state oracles, $A|0\ensuremath{\rangle}\ensuremath{\mapsto}|A\ensuremath{\rangle}$ and $B|0\ensuremath{\rangle}\ensuremath{\mapsto}|B\ensuremath{\rangle}$, such that $H=\ensuremath{\langle}B|U|A\ensuremath{\rangle}$. Our strategy for applying this method to the Sachdev-Ye-Kitaev model involves realizing $B$ using only Hadamard gates and realizing $A$ as a random quantum circuit.

DOI
10.1103/physreva.99.040301

Related papers

Sachdev-Ye-Kitaev model on a noisy quantum computer

Muhammad Asaduzzaman, Raghav G. Jha, Bharath Sambasivam

We study the Sachdev-Ye-Kitaev (SYK) model---an important toy model for quantum gravity on IBM's superconducting qubit quantum computers. By using a graph-coloring algorithm to minimize the number of commuting clusters of terms in the qubitized Hamiltonian, we find the gate complexity of the time evolution using the first-order product formula for $N$ Majorana fermions is $\mathcal{O}({N}^{5}{J}^{2}{t}^{2}/\ensuremath{\epsilon})$ where $J$ is the dimensionful coupling parameter, $t$ is the evolution time, and $\ensuremath{\epsilon}$ is the desired precision. With this improved resource requirement, we perform the time evolution for $N=6$, 8 with maximum two-qubit circuit depth and gate count of 343. We perform different error mitigation schemes on the noisy hardware results and find good agreement with the exact diagonalization results on classical computers and noiseless simulators. In particular, we compute vacuum return probability after time $t$ and out-of-time order correlators which is a standard observable of quantifying the chaotic nature of quantum systems.

arXiv:2311.17991•Physical review. D/Physical review. D.•May 1, 2024•26 citations

Quantum Resource Estimation for Simulating the SYK Model with Trotterization, qDRIFT, and Asymmetric Qubitization

Brian Goldsmith, Larissa Kroell, Nishna Aerabati

The Sachdev-Ye-Kitaev (SYK) model has been identified as a promising candidate to run on early fault-tolerant quantum computers due to the relatively modest resources required to probe non-trivial physics (namely holographic duality and AdS/CFT correspondence). As such, it is crucial that the details of how to run such a simulation are well understood. Using PsiQuantum's Construct platform, we implement and analyze three different approaches to simulate the SYK model: Trotterization, qDRIFT, and asymmetric qubitization with Quantum Signal Processing. We provide an open-source library containing implementations for SYK simulation using all three methods, which we use to obtain quantum resource estimates for qubit and T gate count as functions of the number of Majorana modes and precision. We find that while qDRIFT and Trotterization benefit from a lower qubit count, the large number of T gates required lead to asymmetric qubitization being advantageous in most cases. This reinforces previous theoretical considerations. We intend both the implementations and the estimates to be useful for researchers to continue to study the SYK model and understand how the techniques and resources vary.

quant-pharXiv:2608.27573•Aug 31, 2026•0 citations

Hamiltonian Simulation by Qubitization

Guang Hao Low, Isaac L. Chuang

We present the problem of approximating the time-evolution operator e−iH^t to error ϵ , where the Hamiltonian H^=(⟨G|⊗I^)U^(|G⟩⊗I^) is the projection of a unitary oracle U^ onto the state |G⟩ created by another unitary oracle. Our algorithm solves this with a query complexity O(t+log⁡(1/ϵ)) to both oracles that is optimal with respect to all parameters in both the asymptotic and non-asymptotic regime, and also with low overhead, using at most two additional ancilla qubits. This approach to Hamiltonian simulation subsumes important prior art considering Hamiltonians which are d -sparse or a linear combination of unitaries, leading to significant improvements in space and gate complexity, such as a quadratic speed-up for precision simulations. It also motivates useful new instances, such as where H^ is a density matrix. A key technical result is `qubitization', which uses the controlled version of these oracles to embed any H^ in an invariant SU(2) subspace. A large class of operator functions of H^ can then be computed with optimal query complexity, of which e−iH^t is a special case.

arXiv:1610.06546•Quantum•Jul 12, 2019•854 citations

Hamiltonian Simulation by Qubitization

Guang Hao Low, Isaac L. Chuang

We present the problem of approximating the time-evolution operator e-iHt to error ϵ, where the Hamiltonian H = ((G|⊗I)U (|Gi⊗I) is the projection of a unitary oracle U onto the state |Gi created by another unitary oracle. Our algorithm solves this with a query complexity O ( t + log(1/ϵ) ) to both oracles that is optimal with respect to all parameters in both the asymptotic and non-asymptotic regime, and also with low overhead, using at most two additional ancilla qubits. This approach to Hamiltonian simulation subsumes important prior art considering Hamiltonians which are d-sparse or a linear combination of unitaries, leading to significant improvements in space and gate complexity, such as a quadratic speed-up for precision simulations. It also motivates useful new instances, such as where H is a density matrix. A key technical result is 'qubitization', which uses the controlled version of these oracles to embed any H in an invariant SU(2) subspace. A large class of operator functions of H can then be computed with optimal query complexity, of which e-iHt is a special case.

•Jul 1, 2019•652 citations

Parallel Quantum Algorithm for Hamiltonian Simulation

Zhicheng Zhang, Qisheng Wang, Mingsheng Ying

We study how parallelism can speed up quantum simulation. A parallel quantum algorithm is proposed for simulating the dynamics of a large class of Hamiltonians with good sparse structures, called uniform-structured Hamiltonians, including various Hamiltonians of practical interest like local Hamiltonians and Pauli sums. Given the oracle access to the target sparse Hamiltonian, in both query and gate complexity, the running time of our parallel quantum simulation algorithm measured by the quantum circuit depth has a doubly (poly-)logarithmic dependence polylog⁡log⁡(1/ϵ) on the simulation precision ϵ . This presents an exponential improvement over the dependence polylog⁡(1/ϵ) of previous optimal sparse Hamiltonian simulation algorithm without parallelism. To obtain this result, we introduce a novel notion of parallel quantum walk, based on Childs' quantum walk. The target evolution unitary is approximated by a truncated Taylor series, which is obtained by combining these quantum walks in a parallel way. A lower bound Ω(log⁡log⁡(1/ϵ)) is established, showing that the ϵ -dependence of the gate depth achieved in this work cannot be significantly improved.Our algorithm is applied to simulating three physical models: the Heisenberg model, the Sachdev-Ye-Kitaev model and a quantum chemistry model in second quantization. By explicitly calculating the gate complexity for implementing the oracles, we show that on all these models, the total gate depth of our algorithm has a polylog⁡log⁡(1/ϵ) dependence in the parallel setting.

arXiv:2105.11889•Quantum•Jan 15, 2024•12 citations

Digital Quantum Simulation of Minimal AdS/CFT

Laura García-Álvarez, I. L. Egusquiza, Lucas Lamata +3 more

We propose the digital quantum simulation of a minimal AdS/CFT model in controllable quantum platforms. We consider the Sachdev-Ye-Kitaev model describing interacting Majorana fermions with randomly distributed all-to-all couplings, encoding nonlocal fermionic operators onto qubits to efficiently implement their dynamics via digital techniques. Moreover, we also give a method for probing nonequilibrium dynamics and the scrambling of information. Finally, our approach serves as a protocol for reproducing a simplified low-dimensional model of quantum gravity in advanced quantum platforms as trapped ions and superconducting circuits.

arXiv:1607.08560•Physical Review Letters•Jul 25, 2017•152 citations
We show that one can quantum simulate the dynamics of a Sachdev-Ye-Kitaev model with NNN Majorana modes for time ttt to precision \ensuremathϵ\ensuremath{\epsilon}\ensuremathϵ with gate complexity O(N7/2t+N5/2tpolylog(N/\ensuremathϵ))O({N}^{7/2}t+{N}^{5/2}t\phantom{\rule{0.16em}{0ex}}\mathrm{polylog}(N/\ensuremath{\epsilon}))O(N7/2t+N5/2tpolylog(N/\ensuremathϵ)). In addition to scaling sublinearly in the number of Hamiltonian terms, this gate complexity represents an exponential improvement in 1/\ensuremathϵ1/\ensuremath{\epsilon}1/\ensuremathϵ and large polynomial improvement in NNN and ttt over prior state-of-the-art algorithms which scale as O(N10t2/\ensuremathϵ)O({N}^{10}{t}^{2}/\ensuremath{\epsilon})O(N10t2/\ensuremathϵ). Our approach involves a variant of the qubitization technique in which we encode the Hamiltonian HHH as an asymmetric projection of a signal oracle UUU onto two different signal states prepared by state oracles, A∣0\ensuremath⟩\ensuremath↦∣A\ensuremath⟩A|0\ensuremath{\rangle}\ensuremath{\mapsto}|A\ensuremath{\rangle}A∣0\ensuremath⟩\ensuremath↦∣A\ensuremath⟩ and B∣0\ensuremath⟩\ensuremath↦∣B\ensuremath⟩B|0\ensuremath{\rangle}\ensuremath{\mapsto}|B\ensuremath{\rangle}B∣0\ensuremath⟩\ensuremath↦∣B\ensuremath⟩, such that H=\ensuremath⟨B∣U∣A\ensuremath⟩H=\ensuremath{\langle}B|U|A\ensuremath{\rangle}H=\ensuremath⟨B∣U∣A\ensuremath⟩. Our strategy for applying this method to the Sachdev-Ye-Kitaev model involves realizing BBB using only Hadamard gates and realizing AAA as a random quantum circuit.
We study the Sachdev-Ye-Kitaev (SYK) model---an important toy model for quantum gravity on IBM's superconducting qubit quantum computers. By using a graph-coloring algorithm to minimize the number of commuting clusters of terms in the qubitized Hamiltonian, we find the gate complexity of the time evolution using the first-order product formula for NNN Majorana fermions is $\mathcal{O}({N}^{5}{J}^{2}{t}^{2}/\ensuremath{\epsilon})$ where JJJ is the dimensionful coupling parameter, ttt is the evolution time, and $\ensuremath{\epsilon}$ is the desired precision. With this improved resource requirement, we perform the time evolution for N=6N=6N=6, 8 with maximum two-qubit circuit depth and gate count of 343. We perform different error mitigation schemes on the noisy hardware results and find good agreement with the exact diagonalization results on classical computers and noiseless simulators. In particular, we compute vacuum return probability after time ttt and out-of-time order correlators which is a standard observable of quantifying the chaotic nature of quantum systems.
The Sachdev-Ye-Kitaev (SYK) model has been identified as a promising candidate to run on early fault-tolerant quantum computers due to the relatively modest resources required to probe non-trivial physics (namely holographic duality and AdS/CFT correspondence). As such, it is crucial that the details of how to run such a simulation are well understood. Using PsiQuantum's Construct platform, we implement and analyze three different approaches to simulate the SYK model: Trotterization, qDRIFT, and asymmetric qubitization with Quantum Signal Processing. We provide an open-source library containing implementations for SYK simulation using all three methods, which we use to obtain quantum resource estimates for qubit and T gate count as functions of the number of Majorana modes and precision. We find that while qDRIFT and Trotterization benefit from a lower qubit count, the large number of T gates required lead to asymmetric qubitization being advantageous in most cases. This reinforces previous theoretical considerations. We intend both the implementations and the estimates to be useful for researchers to continue to study the SYK model and understand how the techniques and resources vary.
We present the problem of approximating the time-evolution operator e−iH^t to error ϵ , where the Hamiltonian H^=(⟨G|⊗I^)U^(|G⟩⊗I^) is the projection of a unitary oracle U^ onto the state |G⟩ created by another unitary oracle. Our algorithm solves this with a query complexity O(t+log⁡(1/ϵ)) to both oracles that is optimal with respect to all parameters in both the asymptotic and non-asymptotic regime, and also with low overhead, using at most two additional ancilla qubits. This approach to Hamiltonian simulation subsumes important prior art considering Hamiltonians which are d -sparse or a linear combination of unitaries, leading to significant improvements in space and gate complexity, such as a quadratic speed-up for precision simulations. It also motivates useful new instances, such as where H^ is a density matrix. A key technical result is `qubitization', which uses the controlled version of these oracles to embed any H^ in an invariant SU(2) subspace. A large class of operator functions of H^ can then be computed with optimal query complexity, of which e−iH^t is a special case.
We present the problem of approximating the time-evolution operator e-iHt to error ϵ, where the Hamiltonian H = ((G|⊗I)U (|Gi⊗I) is the projection of a unitary oracle U onto the state |Gi created by another unitary oracle. Our algorithm solves this with a query complexity O ( t + log(1/ϵ) ) to both oracles that is optimal with respect to all parameters in both the asymptotic and non-asymptotic regime, and also with low overhead, using at most two additional ancilla qubits. This approach to Hamiltonian simulation subsumes important prior art considering Hamiltonians which are d-sparse or a linear combination of unitaries, leading to significant improvements in space and gate complexity, such as a quadratic speed-up for precision simulations. It also motivates useful new instances, such as where H is a density matrix. A key technical result is 'qubitization', which uses the controlled version of these oracles to embed any H in an invariant SU(2) subspace. A large class of operator functions of H can then be computed with optimal query complexity, of which e-iHt is a special case.
We study how parallelism can speed up quantum simulation. A parallel quantum algorithm is proposed for simulating the dynamics of a large class of Hamiltonians with good sparse structures, called uniform-structured Hamiltonians, including various Hamiltonians of practical interest like local Hamiltonians and Pauli sums. Given the oracle access to the target sparse Hamiltonian, in both query and gate complexity, the running time of our parallel quantum simulation algorithm measured by the quantum circuit depth has a doubly (poly-)logarithmic dependence polylog⁡log⁡(1/ϵ) on the simulation precision ϵ . This presents an exponential improvement over the dependence polylog⁡(1/ϵ) of previous optimal sparse Hamiltonian simulation algorithm without parallelism. To obtain this result, we introduce a novel notion of parallel quantum walk, based on Childs' quantum walk. The target evolution unitary is approximated by a truncated Taylor series, which is obtained by combining these quantum walks in a parallel way. A lower bound Ω(log⁡log⁡(1/ϵ)) is established, showing that the ϵ -dependence of the gate depth achieved in this work cannot be significantly improved.Our algorithm is applied to simulating three physical models: the Heisenberg model, the Sachdev-Ye-Kitaev model and a quantum chemistry model in second quantization. By explicitly calculating the gate complexity for implementing the oracles, we show that on all these models, the total gate depth of our algorithm has a polylog⁡log⁡(1/ϵ) dependence in the parallel setting.
We propose the digital quantum simulation of a minimal AdS/CFT model in controllable quantum platforms. We consider the Sachdev-Ye-Kitaev model describing interacting Majorana fermions with randomly distributed all-to-all couplings, encoding nonlocal fermionic operators onto qubits to efficiently implement their dynamics via digital techniques. Moreover, we also give a method for probing nonequilibrium dynamics and the scrambling of information. Finally, our approach serves as a protocol for reproducing a simplified low-dimensional model of quantum gravity in advanced quantum platforms as trapped ions and superconducting circuits.