Inchworm quasi Monte Carlo for quantum impurities

The inchworm expansion is a promising approach to solving strongly correlated quantum impurity models due to its reduction of the sign problem in real and imaginary time. However, inchworm Monte Carlo is computationally expensive, converging as $1/\sqrt{N}$ where $N$ is the number of samples. We show that the imaginary time integration is amenable to quasi Monte Carlo, with enhanced $1/N$ convergence, by mapping the Sobol low-discrepancy sequence from the hypercube to the simplex with the so-called Root transform. This extends the applicability of the inchworm method to, e.g., multi-orbital Anderson impurity models with off-diagonal hybridization, relevant for materials simulation, where continuous time hybridization expansion Monte Carlo has a severe sign problem.