Quantum simulators enable the exploration of complex quantum phenomena in condensed-matter systems by reproducing their dynamics on controllable quantum devices. However, experimental constraints often restrict the class of Hamiltonians that can be realized natively. Hamiltonian engineering addresses this limitation by expanding the set of accessible target Hamiltonians from a fixed system Hamiltonian defined by the hardware. We introduce a new framework for fermionic Hamiltonian engineering based on conjugating free evolution under the system Hamiltonian with sequences of experimentally feasible local fermionic unitaries. The required sequences and free-evolution times are obtained efficiently via a linear program. By interleaving system evolution with these local unitaries, our method realizes effective time evolution under a broad class of target Hamiltonians, with intrinsic robustness to finite-pulse-time errors. In particular, we demonstrate that arbitrary complex tunnelling coefficients can be realized, constrained only by the connectivity of the underlying system Hamiltonian. We illustrate this capability by engineering the dynamics of the non-interacting Harper-Hofstadter model on a 1088-mode lattice and an interacting Fermi-Hubbard chain with complex tunnelling coefficients. By construction, our approach avoids the continuous energy absorption inherent to Floquet engineering.
Quantum simulators enable the exploration of complex quantum phenomena in condensed-matter systems by reproducing their dynamics on controllable quantum devices. However, experimental constraints often restrict the class of Hamiltonians that can be realized natively. Hamiltonian engineering addresses this limitation by expanding the set of accessible target Hamiltonians from a fixed system Hamiltonian defined by the hardware. We introduce a new framework for fermionic Hamiltonian engineering based on conjugating free evolution under the system Hamiltonian with sequences of experimentally feasible local fermionic unitaries. The required sequences and free-evolution times are obtained efficiently via a linear program. By interleaving system evolution with these local unitaries, our method realizes effective time evolution under a broad class of target Hamiltonians, with intrinsic robustness to finite-pulse-time errors. In particular, we demonstrate that arbitrary complex tunnelling coefficients can be realized, constrained only by the connectivity of the underlying system Hamiltonian. We illustrate this capability by engineering the dynamics of the non-interacting Harper-Hofstadter model on a 1088-mode lattice and an interacting Fermi-Hubbard chain with complex tunnelling coefficients. By construction, our approach avoids the continuous energy absorption inherent to Floquet engineering.
Quantum simulators enable the exploration of complex quantum phenomena in condensed-matter systems by reproducing their dynamics on controllable quantum devices. However, experimental constraints often restrict the class of Hamiltonians that can be realized natively. Hamiltonian engineering addresses this limitation by expanding the set of accessible Hamiltonians beyond the fixed system Hamiltonian provided by the hardware. We introduce a framework for fermionic Hamiltonian engineering based on interleaving free evolution under the system Hamiltonian with a sequence of experimentally feasible local fermionic unitaries. The required sequences and free-evolution times are obtained efficiently through a linear program, enabling effective time evolution under a broad class of target Hamiltonians. In particular, our method realizes arbitrary complex tunnelling coefficients, constrained only by the connectivity of the underlying system Hamiltonian. We further show that symmetries of the Hamiltonian engineering task can be exploited to drastically reduce the cost of computing the control parameters, and that errors arising from finite pulse times can be compensated within the same framework through a simple modification of the linear program. To our knowledge, this is the first automated method to realize locally tunable artificial gauge fields in interacting fermionic lattice models using only site-resolved energy shifts as control. We demonstrate the method by engineering the dynamics of interacting Fermi-Hubbard chains with complex tunnelling coefficients and of the Harper-Hofstadter model on a 1088-mode lattice. For the latter, we exploit the translation symmetry of the model to determine a provably optimal six-pulse sequence that produces the desired dynamics independent of the lattice size.
Implementing the time evolution under a desired target Hamiltonian is critical for various applications in quantum science. Due to the exponential increase of parameters in the system size and due to experimental imperfections, this task can be challenging in quantum many-body settings. We introduce an efficient and robust scheme to engineer arbitrary local many-body Hamiltonians. To this end, our scheme applies single-qubit π or π / 2 pulses to an always-on system Hamiltonian, which we assume to be native to a given platform. These sequences are constructed by efficiently solving a linear program (LP) which minimizes the total evolution time. In this way, we can engineer target Hamiltonians that are only limited by the locality of the interactions in the system Hamiltonian. Based on average Hamiltonian theory and by using robust composite pulses, we make our schemes robust against errors, including finite-pulse-time errors and various control errors. To demonstrate the performance of our scheme, we provide numerical simulations. In particular, we solve the Hamiltonian-engineering problem on a laptop for arbitrary two-local Hamiltonians on a two-dimensional square lattice with 196 qubits in only 60 s. Moreover, we simulate the engineering of general Heisenberg Hamiltonians from Ising Hamiltonians with imperfect single-qubit pulses for smaller system sizes, and achieve a fidelity larger than 99.9%, which is orders of magnitude better than nonrobust implementations.
Hamiltonian engineering (HE) simulates quantum dynamics under a target Hamiltonian using a native entangling system Hamiltonian and restricted control, with applications from quantum gate design to analog quantum simulation. Observing that, for many systems, formulating HE as a linear program only requires specific commutation relations between the control pulses and the system Hamiltonian, we provide a unified framework for broad classes of 2-local qubit, qudit, and fermionic systems. We reformulate time-optimal HE as a complex k-cut polytope problem, with k the number of distinct phases in these commutation relations, and prove it NP-complete for every finite k. We therefore relax the polytope to the elliptope and apply a Krivine-type rounding, yielding pulses informed by the system and target Hamiltonians, which we then use in a linear program. We derive how many informed pulses are necessary and sufficient for the linear program to be reliably feasible. Additionally mixing in uniformly sampled pulses guarantees a solution close to the relaxation value from O(m) pulses, with m the number of interaction terms. Together with upper and lower bounds on the optimal quantum run time, this yields an O(m) approximation ratio. In benchmarks on a fully connected Ising model, a chiral clock model for qudits, and the fermionic Harper-Hofstadter model, our approach attains near-optimal quantum run times wherever the optimum is computable, reaches run times that saturate with the lattice size in the fermionic case, and outperforms state-of-the-art methods with O(m) pulses. This establishes a unified approach with provable guarantees to the automatic programming of analog quantum simulators.
Implementing the time evolution under a desired target Hamiltonian is critical for various applications in quantum science. Due to the exponential increase in the number of parameters with system size and experimental imperfections, this task can be challenging in quantum many-body settings. We introduce an efficient and robust scheme to engineer arbitrary local many-body Hamiltonians. To this end, our scheme applies single-qubit π or π/2 pulses to an always-on system Hamiltonian, which we assume to be native to a given platform. These sequences are constructed by efficiently solving a linear program (LP) which minimizes the total evolution time. In this way, we can engineer target Hamiltonians that are only limited by the locality of the interactions in the system Hamiltonian. Based on average Hamiltonian theory and using robust composite pulses, we make our schemes robust against errors, including finite pulse time errors and various control errors. To demonstrate the performance of our scheme, we provide numerical simulations. In particular, we solve the Hamiltonian engineering problem on a laptop for arbitrary two-local Hamiltonians on a 2D square lattice with 196 qubits in only 60 seconds. Moreover, we simulate the engineering of general Heisenberg Hamiltonians from Ising Hamiltonians using imperfect single-qubit pulses for smaller system sizes and achieve a fidelity exceeding 99.9%, which is orders of magnitude better than non-robust implementations.
Implementing the time evolution under a desired target Hamiltonian is critical for various applications in quantum science. Due to the exponential increase of parameters in the system size and due to experimental imperfections, this task can be challenging in quantum many-body settings. We introduce an efficient and robust scheme to engineer arbitrary local many-body Hamiltonians. To this end, our scheme applies single-qubit π or π/2 pulses to an always-on system Hamiltonian, which we assume to be native to a given platform. These sequences are constructed by efficiently solving a linear program (LP) which minimizes the total evolution time. In this way, we can engineer target Hamiltonians that are only limited by the locality of the interactions in the system Hamiltonian. Based on average Hamiltonian theory and by using robust composite pulses, we make our schemes robust against errors, including finite-pulse-time errors and various control errors. To demonstrate the performance of our scheme, we provide numerical simulations. In particular, we solve the Hamiltonian-engineering problem on a laptop for arbitrary two-local Hamiltonians on a two-dimensional square lattice with 196 qubits in only 60 s. Moreover, we simulate the engineering of general Heisenberg Hamiltonians from Ising Hamiltonians with imperfect single-qubit pulses for smaller system sizes, and achieve a fidelity larger than 99.9%, which is orders of magnitude better than nonrobust implementations.
Analog quantum simulators emulate complex many-body dynamics through native continuous-time evolution under hardware-defined interactions. Yet once a platform is specified, its interaction structure is largely fixed by the underlying hardware, restricting the Hamiltonians that can be realized and limiting programmability. Here we introduce universal analog quantum simulation (UAQS), a hybrid framework that systematically expands the range of accessible quantum evolutions within a given analog platform. UAQS employs optimized continuous-time control fields to engineer target dynamics directly, avoiding decomposition into discrete gate sequences. By preserving native analog evolution while extending the set of achievable Hamiltonians, UAQS transforms fixed-interaction analog devices into programmable simulators. Numerical studies on representative architectures, including superconducting circuits and Rydberg-atom arrays, show that UAQS accurately reproduces non-trivial many-body dynamics beyond the intrinsic interaction structure of the hardware. These results establish UAQS as a practical route toward programmable analog quantum simulation.