Irene Papaefstathiou,Daniel Malz,J. Ignacio Cirac,Mari Carmen Bañuls
Abstract:We present a numerical method to simulate a system of multiple emitters coupled to a non-interacting bath, in any parameter regime. Our method relies on a Block Lanczos transformation that maps the whole system onto a strip geometry, whose width is given by the number of emitters. Utilizing the spatial symmetries of the problem and identifying the relevant range of energies of the bath we achieve a more efficient simulation, which we perform using tensor network techniques. As a demonstration, we study the collective emission from multiple emitters coupled to a square lattice of bosons and observe how the departure from Markovianity as coupling strength and emitter number is increased prevents collective radiation. We also simulate the dynamic preparation of an excitation in a bound state from a multi-excitation initial state. Our work opens new possibilities for the systematic exploration of non-Markovian effects in the dynamics and equilibrium properties of multi-emitter systems. Furthermore, it can easily be extended to other setups, including finite bath temperature or impurities coupled to fermionic environments.
What problem does this paper attempt to address?
### Problems the paper attempts to solve
This paper aims to solve the problem of efficient numerical simulation when multi - emitter systems interact with non - Markovian environments. Specifically, the researchers proposed a method based on Block Lanczos transformation, which maps the entire system onto a strip - like geometric structure, whose width is determined by the number of emitters. By exploiting the spatial symmetry of the problem and identifying the relevant energy range in the bath, the researchers achieved more efficient simulations and used tensor network techniques for calculations.
### Key problems and methods
1. **Non - Markovian dynamics of multi - emitter systems**:
- The researchers focused on the situation where multiple emitters are coupled to non - interacting baths, especially the dynamic behavior in any parameter region.
- The traditional Markov approximation is valid under the conditions of weak coupling and fast - decaying baths, but under the conditions of strong coupling or slow - decaying baths, non - Markovian effects become significant, making numerical simulations very challenging.
2. **Block Lanczos transformation**:
- Through Block Lanczos transformation, the researchers mapped the system onto a strip - like geometric structure, where each emitter corresponds to a "rung".
- This transformation enables the system to be represented as a quasi - one - dimensional chain, so that efficient simulations can be carried out using tensor network techniques such as matrix product states (MPS).
3. **Symmetry enhancement**:
- By exploiting the spatial symmetry of the bath, the researchers further optimized the algorithm, reducing the width of the strip, thereby reducing the computational cost.
- For example, in a two - dimensional square lattice, by selecting a specific axis of symmetry, the bath modes can be grouped, thereby reducing the width of each strip.
4. **Energy truncation**:
- In order to improve the simulation efficiency, the researchers truncated the bath modes in the energy space, only retaining the modes close to the system transition frequencies.
- This truncation method limits the hopping rate to a smaller range, so that shorter chains and larger Trotter step sizes can be used, improving the simulation efficiency.
### Application examples
1. **Non - Markovian effects of collective radiation**:
- By simulating the situation where multiple emitters are coupled to a two - dimensional bosonic lattice, the researchers observed that as the coupling strength and the number of emitters increase, non - Markovian effects lead to the suppression of collective radiation.
2. **Formation of dynamically excited states**:
- The researchers also simulated the process of dynamically preparing a bound state from a multi - excited initial state, demonstrating the influence of non - Markovian effects in this process.
### Conclusion
The methods proposed in this paper provide new possibilities for systematically exploring non - Markovian effects in multi - emitter systems. These methods are not only applicable to current research, but can also be extended to other settings, such as baths at finite temperatures or impurity systems coupled to fermionic environments.