We introduce a novel construction of a contour deformation within the framework of Loop-Tree Duality for the numerical computation of loop integrals featuring threshold singularities in momentum space. The functional form of our contour deformation automatically satisfies all constraints without the need for fine-tuning. We demonstrate that our construction is systematic and efficient by applying it to more than 100 examples of finite scalar integrals featuring up to six loops. We also showcase a first step towards handling non-integrable singularities by applying our work to one-loop infrared divergent scalar integrals and to the one-loop amplitude for the ordered production of two and three photons. This requires the combination of our contour deformation with local counterterms that regulate soft, collinear and ultraviolet divergences. This work is an important step towards computing higher-order corrections to relevant scattering cross-sections in a fully numerical fashion. arXiv:1912.09291v1 [hep-ph] 19 Dec 2019 7 Results 53 7.1 Multi-loop finite integrals 54 7.2 Divergent one-loop four-and five-point scalar integrals 69 7.3 One-loop amplitude for qq → γ 1 γ 2 and qq → γ 1 γ 2 γ 3 69 -i -8 Conclusion 74 9 Acknowledgements 75 A Loop-Tree Duality example at two loops 76 B Expression for the qq → γ 1 γ 2 γ 3 amplitude and its counterterms 79 C Loop-Tree Duality with raised propagators 80 1 As discussed in ref.[83], our final expression in eq. (2.5) is also correct in the case of complex-valued external momenta, due to the fact that the right-most column of the matrix appearing in eq. (2.4) does not include the imaginary part Im[p 0 i ] of the external momenta. We note however, that the correct interpretation of the absence of this term in eq. (2.4) for complex-valued external kinematics is that the energy integrals are no longer performed along the real line but instead along a path including only one out of the two complex energy solutions of each propagator.3)The quantity ∇ k η( k), henceforth denoted as ∇η, is the outward pointing normal vector to the surface η( k) = 0. The contour deformation is defined in the (3n)-dimensional complex space and we parametrise it as k − i κ( k). It must satisfy constraints affecting two of its key characteristics, the direction and magnitude of the vector field κ( k):Direction: The deformation vector κ( k) must induce a sign of the imaginary part of the E-surface equation that matches the sign enforced by the causal prescription whenever k lies on a singular E-surfaces. This imposes conditions on the direction of the vector field κ( k). We derive these conditions by comparing the sign of the LTD prescription