Search papers, labs, and topics across Lattice.
This paper introduces diagonal-budgeted Trotterization, a novel method for classical simulation of quantum Hamiltonian dynamics that maintains diagonal sparsity and controls fidelity loss. The authors implemented HamSim, which utilizes a specialized data layout and optimized C++/CUDA kernels to significantly enhance performance and efficiency compared to existing frameworks like Qiskit-Aer. Benchmark results reveal that HamSim achieves remarkable speedups of up to 1,269x on CPUs and 178x on GPUs, while preserving near-perfect fidelity without the need for exponential steps in Trotterization.
HamSim achieves unprecedented speedups in classical quantum Hamiltonian simulation while maintaining near-perfect fidelity, challenging the limits of traditional methods.
Efficient classical simulation of quantum Hamiltonian dynamics is often bottlenecked by exponential state growth and the overhead of generic sparse linear algebra. We introduce diagonal-budgeted Trotterization, a structure-aware strategy that decomposes Hamiltonians into factors preserving diagonal sparsity while tightly controlling fidelity loss. Our implementation, HamSim, utilizes a compact diagonal-sparse data layout and specialized C++/CUDA kernels to bypass the overheads of generic formats like CSR. By leveraging SIMD vectorization, multithreading, and GPU acceleration, HamSim achieves high performance across heterogeneous architectures. Benchmarks on the HamLib suite show that HamSim significantly outperforms Qiskit-Aer. On CPUs, HamSim attains speedups of $182$--$1,269\times$ on optimization instances (TSP, MaxCut) and $4.8$--$841\times$ on physical models (TFIM, Heisenberg). On GPUs, it achieves up to $178\times$ speedup for $12$--$16$ qubit problems. Unlike traditional Trotterization, HamSim maintains near-perfect fidelity without requiring exponential steps. This demonstrates that diagonal-aware numerical kernels provide a scalable foundation for high-fidelity classical Hamiltonian simulation.
