Stable Evaluation of Lefschetz Thimble Intersection Numbers: Towards Real-Time Path Integrals

We introduce a robust numerical method for determining intersection numbers of Lefschetz thimbles in multivariable settings. Our approach employs the multiple shooting method to solve the upward flow equations from the saddle points to the original integration cycle, which also enables us to determine the signs of the intersection numbers. The method demonstrates stable and reliable performance, and has been tested for systems with up to $20$ variables, which can be further extended by adopting quadruple-precision arithmetic. We determine intersection numbers for several complex saddle points in a discretized path integral, providing new insights into the structure of real-time path integrals. The proposed method is broadly applicable to a wide range of problems involving oscillatory integrals in physics and mathematics.