Seth K. Asante,Sebastian Steinhaus
Abstract:Numerical computations and methods have become increasingly crucial in the study of spin foam models across various regimes. This paper adds to this field by introducing new algorithms based on tensor network methods for computing amplitudes, focusing on topological SU(2) BF and Lorentzian EPRL spin foam models. By reorganizing the sums and tensors involved, vertex amplitudes are recast as a sequence of matrix contractions. This reorganization significantly reduces computational complexity and memory usage, allowing for scalable and efficient computations of the amplitudes for larger representation labels on standard consumer hardware--previously infeasible due to the computational demands of high-valent tensors.
We apply these tensor network algorithms to analyze the characteristics of various vertex configurations, including Regge and vector geometries for the SU(2) BF theory, demonstrating consistent scaling behavior and differing oscillation patterns. Our benchmarks reveal substantial improvements in computational time and memory allocations, especially for large representation labels. Additionally, these tensor network methods are applicable to generic 2-complexes with multiple vertices, where we introduce partial-coherent vertex amplitudes to streamline the computations. The implementation of these algorithms is available on GitHub for further exploration and use.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the efficient calculation of amplitudes in the spin - foam model. Specifically, the paper introduces a new algorithm based on the tensor network method for calculating vertex amplitudes, especially for the topological \( \text{SU}(2) \) BF model and the Lorentz EPRL spin - foam model. By reorganizing the sums and tensors involved, the vertex amplitudes are rewritten as a series of matrix contractions. This reorganization significantly reduces the computational complexity and memory usage, enabling the calculation of amplitudes with larger representation labels on standard consumer - level hardware, which was difficult to achieve in the past due to the computational requirements of high - order tensors.
### Main problems
1. **Computational complexity and memory usage**: Traditional spin - foam model calculation methods face serious computational complexity and memory usage problems when dealing with large representation labels, especially in the case of high - order tensors.
2. **Scalability**: Existing calculation methods are inefficient when dealing with multiple vertices or large - scale calculations, limiting the exploration of the properties of the spin - foam model and the semiclassical regime.
3. **Limitations of existing tools**: Although there are some numerical libraries (such as sl2cfoam and sl2cfoam - next) that can be used to calculate the coherent amplitudes of the Lorentz EPRL and \( \text{SU}(2) \) BF models, these tools still have performance bottlenecks when dealing with large representation labels.
### Solutions
The paper proposes a method based on tensor networks to optimize the calculation through the following steps:
1. **Tensor network representation**: Represent the vertex amplitudes as matrix contractions in the tensor network instead of directly calculating high - order tensors.
2. **Matrix contraction**: Reduce the computational complexity and memory usage by decomposing high - order tensors into low - order tensors (especially matrices).
3. **Partially coherent vertex amplitudes**: Introduce partially coherent vertex amplitudes to further simplify the calculation process, making it suitable for general 2 - complexes containing multiple vertices.
4. **Benchmarking**: Verify the significant improvement of the new method in terms of calculation time and memory allocation through benchmarking, especially in the case of large representation labels.
### Specific technical details
- **Matrix representation of the {15j}-symbol**: Represent the {15j}-symbol as a matrix product of {6j}-symbols, and each {6j}-symbol can be represented as a matrix, with two tensor indices as open indices.
- **Matrix contraction**: Calculate the vertex amplitudes through matrix contraction instead of directly calculating the contraction of high - order tensors.
- **Partially coherent vertex amplitudes**: For 2 - complexes containing multiple vertices, introduce partially coherent vertex amplitudes to simplify the calculation process.
### Applications and results
- **Numerical results**: The paper shows the numerical results of the new method under different vertex configurations, including Regge geometry and vector geometry, verifying the effectiveness and accuracy of the method.
- **Performance improvement**: Benchmarking shows that the new method has significant improvements in terms of calculation time and memory usage, especially in the case of large representation labels.
In conclusion, this paper solves the problems of high computational complexity and large memory usage in the spin - foam model by introducing a new algorithm based on tensor networks, improving the computational efficiency and scalability.