We compute numerically the Lorentzian Engle-Pereira-Rovelli-Livine (EPRL) spinfoam propagator on a 4-simplex, by adapting the methods of Lefschetz thimble and Markov Chain Monte-Carlo to oscillatory spinfoam integrals. Our method can compute any spinfoam observables at relatively large spins. We obtain the numerical results of the propagators at different spins and demonstrate their consistency with the expected spinfoam semi-classical behavior in the large spin limit. Our results exhibit significant quantum corrections at smaller spins. Our method is reliable and thus can be employed to discover the semi-classical and quantum behaviors of the spinfoam model.
Comment: 26 pages, 8 figures, 11 tables