Effective light cone and digital quantum simulation of interacting bosons

Tomotaka Kuwahara,Tan Van Vu,Keiji Saito

DOI: https://doi.org/10.1038/s41467-024-46501-7

2024-03-22

Abstract:The speed limit of information propagation is one of the most fundamental features in non-equilibrium physics. The region of information propagation by finite-time dynamics is approximately restricted inside the effective light cone that is formulated by the Lieb-Robinson bound. To date, extensive studies have been conducted to identify the shape of effective light cones in most experimentally relevant many-body systems. However, the Lieb-Robinson bound in the interacting boson systems, one of the most ubiquitous quantum systems in nature, has remained a critical open problem for a long time. This study reveals a tight effective light cone to limit the information propagation in interacting bosons, where the shape of the effective light cone depends on the spatial dimension. To achieve it, we prove that the speed for bosons to clump together is finite, which in turn leads to the error guarantee of the boson number truncation at each site. Furthermore, we applied the method to provide a provably efficient algorithm for simulating the interacting boson systems. The results of this study settle the notoriously challenging problem and provide the foundation for elucidating the complexity of many-body boson systems.

Quantum Physics,Quantum Gases,Statistical Mechanics,High Energy Physics - Theory,Mathematical Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the shape of the effective light cone in the interacting boson system and its limitation on the speed of information propagation. Specifically:

1. **Shape of the effective light cone**: The paper explores the shape of the effective light cone of the interacting boson system in different spatial dimensions. This helps to understand the propagation range and speed of information in these systems.

2. **Speed limit of information propagation**: The paper proves that the speed of boson aggregation is limited, thus providing an error guarantee for the truncation of the number of bosons on each lattice site. This result is of great significance for simulating the interacting boson system because it provides an effective method to limit the speed of information propagation.

3. **Efficient algorithm for digital quantum simulation**: Based on the above - mentioned theoretical results, the paper proposes an efficient algorithm for simulating the time evolution of the interacting boson system. This algorithm can achieve high - precision quantum simulation with low gate complexity.

### Specific problem description

- **Background**: The Lieb - Robinson bound provides a fundamental upper limit for the speed of information propagation in non - equilibrium structures during real - time evolution. However, for the interacting boson system, the exact form of the Lieb - Robinson bound has been an unsolved problem.

- **Objective**: The paper aims to reveal a compact effective light cone to limit the information propagation in the interacting boson system. Specifically, the paper proves that the speed of boson aggregation is limited and proposes an efficient algorithm to simulate these systems.

### Main contributions

1. **Compact effective light cone**: The paper proves that in the interacting boson system, the shape of the effective light cone depends on the spatial dimension. This provides a theoretical basis for understanding the speed and range of information propagation.

2. **Finiteness of the speed of boson aggregation**: The paper proves that the speed of boson aggregation is limited, which provides an error guarantee for the truncation of the number of bosons on each lattice site.

3. **Efficient simulation algorithm**: Based on the above - mentioned theoretical results, the paper proposes an efficient algorithm for simulating the time evolution of the interacting boson system. The gate complexity of this algorithm is \( |\Lambda| t^{D + 1}\text{polylog}(|\Lambda| t / \epsilon) \), where \( |\Lambda| \) is the number of lattice sites in the system, \( t \) is the evolution time, and \( \epsilon \) is the error.

### Mathematical formulas

- **Boson number operator**:
\[
\hat{n}_i = b_i^\dagger b_i
\]
where \( b_i \) and \( b_i^\dagger \) are the boson annihilation and creation operators respectively.

- **Bose - Hubbard Hamiltonian**:
\[
\sum_{\langle i, j \rangle} J_{i,j} (b_i b_j^\dagger + \text{h.c.}) + V
\]
where \( J_{i,j} \) is the transition matrix element, and \( V = f(\{\hat{n}_i\}_{i \in \Lambda}) \) is the boson - boson interaction term.

- **Upper bound of the effective light cone**:
\[
R \propto t^D
\]
where \( D \) is the spatial dimension.

- **Error guarantee**:
\[
\| \rho_0(t) - \rho_0(\tilde{H}, t) \|_1 \leq |\Lambda| e^{-c_3 \left[ \frac{\bar{q}}{(t \log t)^D} \right]^{1/\kappa}}
\]
where \( \bar{q} = (t \log t)^D \log^\kappa \left( \frac{|\Lambda|}{\epsilon} \right) \).

Through these theories and algorithms, the paper provides important tools and methods for understanding and simulating the information propagation in the interacting boson system.