Lefschetz Thimble Decomposition

• Lefschetz thimble decomposition is a framework that deforms complex oscillatory integrals into constant-phase steepest-descent manifolds, simplifying sign-problem evaluations.
• It isolates critical points of the holomorphic action and replaces the original contour with a sum of thimbles weighted by topological intersection numbers, managing residual Jacobian phases.
• The approach has practical implications in finite-density models, lattice simulations, and real-time quantum systems by enabling multi-thimble interference analysis and accurate phase cancellation.

Lefschetz thimble decomposition is a reformulation of oscillatory or complex path integrals in which the original real integration cycle is complexified and deformed into a sum of steepest-descent manifolds attached to critical points of the holomorphic action. In Picard–Lefschetz terms, the original contour is replaced by a homologically equivalent sum of Lefschetz thimbles weighted by integer intersection numbers with the corresponding dual cycles. Its practical importance lies in the sign problem: on each thimble the imaginary part of the action is constant, so the dominant phase oscillations of the original integrand are removed, while the remaining difficulty is shifted to a residual Jacobian phase and, in many physically relevant cases, to interference among multiple contributing thimbles (Cristoforetti et al., 2012, Cristoforetti et al., 2013).

1. Mathematical definition and homological structure

For a theory defined on a real domain $\mathcal D \subset \mathbb R^n$, expectation values have the form

$$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$

The Lefschetz-thimble construction assumes that the fields are complexified, that $S(\phi)$ is holomorphic, and that its critical points $\phi^\sigma$ satisfy

$$\frac{\partial S}{\partial \phi_i}=0, \qquad \det\!\left[\frac{\partial^2 S}{\partial \phi_i \partial \phi_j}\right]\neq 0.$$

Each critical point organizes a steepest-descent integration cycle.

A Lefschetz thimble $\mathcal J_\sigma$ is the union of all flow lines ending at $\phi^\sigma$ under the downward gradient flow. In one standard convention,

$$\frac{d\phi}{d\tau} = - \overline{\frac{\partial S}{\partial \phi}}.$$

The dual thimble $\mathcal K_\sigma$ is defined by the opposite asymptotic condition. The exact decomposition of the original contour is then

$$\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)} = \sum_\sigma m_\sigma \int_{\mathcal J_\sigma} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)},$$

where $$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$0 is the intersection number of the original cycle with the dual thimble. In the one-site fermion model this is written as

$$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$1

with $$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$2 given by the intersection of the original contour and $$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$3 (Cristoforetti et al., 2013, Tanizaki et al., 2015).

Geometrically, a thimble is an $$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$4-dimensional real manifold embedded in $$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$5. Along the flow, $$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$6 is constant while $$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$7 changes monotonically. This is the multidimensional analogue of deforming a contour onto steepest-descent paths in ordinary complex analysis. The decomposition is exact at the level of homology, not merely asymptotic, provided the standard holomorphy and nondegeneracy conditions hold (Cristoforetti et al., 2013).

2. Constant-phase manifolds and the sign problem

The sign problem arises when $$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$8 is complex, so the weight $$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$9 is not a positive measure. Standard importance sampling then fails because cancellations from rapidly oscillating phases render naive Monte Carlo exponentially inefficient. The basic appeal of a thimble is that $S(\phi)$0 is constant on each $S(\phi)$1, so

$S(\phi)$2

has only a global action phase on that thimble. The severe oscillations associated with $S(\phi)$3 are therefore absent from the local sampling weight (Cristoforetti et al., 2012, Cristoforetti et al., 2013).

This does not mean that the sign problem disappears completely. The change of variables from the original coordinates to coordinates on a curved thimble introduces a Jacobian, and its phase gives a residual sign problem. In the review formulation this appears through the measure factor

$S(\phi)$4

and the remaining fluctuation is associated with

$S(\phi)$5

The residual phase is expected to be much milder than the original phase oscillation because the action phase is constant and only the geometry of the embedding contributes (Cristoforetti et al., 2013).

A further subtlety is that the exact answer need not be representable by a single thimble. Some works motivate single-thimble formulations by universality and by dominance of the thimble attached to the global minimum of $S(\phi)$6. For a large class of theories, a single suitable thimble has the same degrees of freedom, symmetries, perturbative expansion, and naive continuum limit as the original theory. However, this is a contingent statement rather than a general theorem about all sign-problem systems. In models with Stokes jumps, determinant zeros, or competing saddle sectors, the physically relevant result can require a genuine multi-thimble sum (Cristoforetti et al., 2013, Cristoforetti et al., 2012).

3. Multi-thimble interference, Stokes phenomena, and Silver Blaze physics

The most explicit demonstrations of genuinely multi-thimble behavior appear in finite-density fermionic toy models. In the one-site repulsive Hubbard model with

$S(\phi)$7

the exact partition function is

$S(\phi)$8

At $S(\phi)$9, the density is the step function $\phi^\sigma$0 for $\phi^\sigma$1, $\phi^\sigma$2 for $\phi^\sigma$3, and $\phi^\sigma$4 for $\phi^\sigma$5. After a Hubbard–Stratonovich transformation, the path integral reduces to a one-dimensional oscillatory integral with effective action

$\phi^\sigma$6

and the saddle points $\phi^\sigma$7 form an infinite family in the low-temperature regime. Their actions satisfy

$\phi^\sigma$8

Thus many thimbles have comparable magnitudes while carrying different phases. The semiclassical sum is

$\phi^\sigma$9

which reproduces the zero-temperature nonanalytic jumps. By contrast, a phase-quenched or one-thimble approximation gives the wrong mean-field-like density

$$\frac{\partial S}{\partial \phi_i}=0, \qquad \det\!\left[\frac{\partial^2 S}{\partial \phi_i \partial \phi_j}\right]\neq 0.$$0

The analysis therefore identifies the physically relevant mechanism not as suppression of a single phase, but as interference among many thimbles (Tanizaki et al., 2015).

The same model makes the associated Stokes structure explicit. Around $$\frac{\partial S}{\partial \phi_i}=0, \qquad \det\!\left[\frac{\partial^2 S}{\partial \phi_i \partial \phi_j}\right]\neq 0.$$1 and $$\frac{\partial S}{\partial \phi_i}=0, \qquad \det\!\left[\frac{\partial^2 S}{\partial \phi_i \partial \phi_j}\right]\neq 0.$$2, the contributing thimble decomposition changes. When $$\frac{\partial S}{\partial \phi_i}=0, \qquad \det\!\left[\frac{\partial^2 S}{\partial \phi_i \partial \phi_j}\right]\neq 0.$$3 or $$\frac{\partial S}{\partial \phi_i}=0, \qquad \det\!\left[\frac{\partial^2 S}{\partial \phi_i \partial \phi_j}\right]\neq 0.$$4, only one thimble contributes and the sign problem is mild. In the intermediate regime $$\frac{\partial S}{\partial \phi_i}=0, \qquad \det\!\left[\frac{\partial^2 S}{\partial \phi_i \partial \phi_j}\right]\neq 0.$$5, many thimbles intersect the original contour, and destructive interference is essential. At half-filling, $$\frac{\partial S}{\partial \phi_i}=0, \qquad \det\!\left[\frac{\partial^2 S}{\partial \phi_i \partial \phi_j}\right]\neq 0.$$6, the phases become integer multiples of $$\frac{\partial S}{\partial \phi_i}=0, \qquad \det\!\left[\frac{\partial^2 S}{\partial \phi_i \partial \phi_j}\right]\neq 0.$$7, so the cancellation disappears and the sign problem is absent. Numerically, one-thimble and three-thimble truncations can produce thermodynamic instabilities, whereas five thimbles already reproduce the exact density well at $$\frac{\partial S}{\partial \phi_i}=0, \qquad \det\!\left[\frac{\partial^2 S}{\partial \phi_i \partial \phi_j}\right]\neq 0.$$8. A practical criterion in the most severe case $$\frac{\partial S}{\partial \phi_i}=0, \qquad \det\!\left[\frac{\partial^2 S}{\partial \phi_i \partial \phi_j}\right]\neq 0.$$9 gives roughly

$\mathcal J_\sigma$0

so the number of required thimbles grows linearly with $\mathcal J_\sigma$1 (Tanizaki et al., 2015, Tanizaki et al., 2016).

The one-dimensional lattice Thirring model at finite density shows the same qualitative pattern in a lattice setting. There the complexified link fields $\mathcal J_\sigma$2 give an effective action

$\mathcal J_\sigma$3

and each critical point is paired one-to-one with a zero of the fermion determinant. At small and large chemical potential the original contour is effectively equivalent to a single thimble, while in the crossover region multiple thimbles are necessary. In the low-temperature limit, the rapid crossover behavior is recovered only after adding multi-thimble contributions with alternating signs and partial cancellations (Fujii et al., 2015).

A common oversimplification is that the thimble method reduces every sign problem to a benign single-thimble computation. The finite-density Hubbard and Thirring analyses show that the dominant obstruction can instead be transferred to the relative phases among several saddle sectors. This suggests that in theories with Silver Blaze behavior, the sign problem is tied to exact phase interference rather than to a single global complex prefactor (Tanizaki et al., 2015, Fujii et al., 2015).

4. Numerical algorithms and approximate manifolds

The formal decomposition has motivated several algorithmic realizations. An early proposal is the Aurora algorithm, which samples a thimble attached to the global minimum of $\mathcal J_\sigma$4 by Langevin-like dynamics constrained to the manifold. The drift follows steepest descent of $\mathcal J_\sigma$5, while the noise is projected onto the tangent space at the saddle and then transported along the flow. The tangent transport equation is

$\mathcal J_\sigma$6

and an Iwasawa-type decomposition of the Hessian is used to construct a more stable orthogonal transport (Cristoforetti et al., 2012).

A distinct strategy maps the curved thimble of the full action to the flat manifold of the corresponding quadratic action near the saddle. In the Metropolis approach to the $\mathcal J_\sigma$7 one-plaquette model,

$\mathcal J_\sigma$8

the algorithm samples Gaussian coordinates and deterministically flows them onto the nonlinear thimble. The numerical results converge to the exact analytic value

$\mathcal J_\sigma$9

and the residual phase is reported not to represent a sign problem in that model (Mukherjee et al., 2013).

The first practical Monte Carlo calculations on a thimble approximation were carried out for the relativistic Bose gas. There the exact thimble is approximated by its tangent space $\phi^\sigma$0, and the antiholomorphic flow is used to move systematically toward the true thimble. Simulations were performed in $\phi^\sigma$1 on lattices up to $\phi^\sigma$2 with $\phi^\sigma$3, and the results for $\phi^\sigma$4 and $\phi^\sigma$5 agreed excellently with worm-algorithm benchmarks and earlier determinations. A transition around $\phi^\sigma$6 was observed. Moving closer to the thimble reduced fluctuations of the imaginary part of the action by factors of roughly $\phi^\sigma$7, $\phi^\sigma$8, and $\phi^\sigma$9 on $$\frac{d\phi}{d\tau} = - \overline{\frac{\partial S}{\partial \phi}}.$$0, $$\frac{d\phi}{d\tau} = - \overline{\frac{\partial S}{\partial \phi}}.$$1, and $$\frac{d\phi}{d\tau} = - \overline{\frac{\partial S}{\partial \phi}}.$$2 at $$\frac{d\phi}{d\tau} = - \overline{\frac{\partial S}{\partial \phi}}.$$3, and the residual phase on $$\frac{d\phi}{d\tau} = - \overline{\frac{\partial S}{\partial \phi}}.$$4 was statistically negligible in the cases studied (Cristoforetti et al., 2013).

To address multimodality at large flow time, the tempered Lefschetz thimble method introduces parallel tempering in flow time. A ladder of flowed manifolds $$\frac{d\phi}{d\tau} = - \overline{\frac{\partial S}{\partial \phi}}.$$5 is simulated jointly, so small-$$\frac{d\phi}{d\tau} = - \overline{\frac{\partial S}{\partial \phi}}.$$6 replicas preserve ergodicity while large-$$\frac{d\phi}{d\tau} = - \overline{\frac{\partial S}{\partial \phi}}.$$7 replicas reduce the sign problem. An HMC implementation on each flowed surface was developed for the Hubbard model on a small lattice; it reproduces the exact Trotterized density $$\frac{d\phi}{d\tau} = - \overline{\frac{\partial S}{\partial \phi}}.$$8, with HMC yielding $$\frac{d\phi}{d\tau} = - \overline{\frac{\partial S}{\partial \phi}}.$$9 and a lower effective cost than Metropolis, about $\mathcal K_\sigma$0 with swaps and less than $\mathcal K_\sigma$1 without swaps. The algorithm includes a momentum-flip mechanism when the molecular-dynamics step approaches zeros of fermion determinants (Fukuma et al., 2019).

The worldvolume tempered Lefschetz thimble method replaces discrete replicas by a continuous union of flowed surfaces,

$\mathcal K_\sigma$2

and performs HMC on this worldvolume. It preserves the ergodicity advantage of tempering while eliminating the need to compute the full flow Jacobian during configuration generation. In the Stephanov model, where complex Langevin is known to suffer from wrong convergence, WV-TLTM agrees well with exact values with controlled statistical errors (Fukuma et al., 2021).

5. Real-time path integrals, quantum cosmology, and intersection-number computation

Lefschetz thimble methods extend beyond Euclidean finite-density problems. In real-time quantum mechanics, the discretized path integral has oscillatory phase $\mathcal K_\sigma$3, and the generalized Lefschetz thimble method uses the anti-holomorphic gradient flow

$\mathcal K_\sigma$4

to deform the integration domain to a flowed manifold $\mathcal K_\sigma$5. In the $\mathcal K_\sigma$6 limit, the manifold decomposes into thimbles attached to critical points. For the double-well potential $\mathcal K_\sigma$7, the relevant saddles are complex classical trajectories rather than real ones, and tunneling is carried by these complex solutions. The weak value

$\mathcal K_\sigma$8

computed from the generalized thimble method matches direct Schrödinger evolution and becomes complex in the tunneling regime, giving an observable diagnostic of the complex saddle dynamics (Nishimura, 2023).

In Lorentzian quantum cosmology, the generalized thimble method has been used as a numerical realization of Picard–Lefschetz theory for the oscillatory gravitational path integral. In a mini-superspace model with metric

$\mathcal K_\sigma$9

the flow deforms the lapse and scale-factor integration variables onto a manifold where phase fluctuations are much smaller. The nonperturbative Monte Carlo analysis confirms that Dirichlet boundary conditions select Vilenkin-type saddles, while Robin boundary conditions with

$$\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)} = \sum_\sigma m_\sigma \int_{\mathcal J_\sigma} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)},$$0

can switch the relevant saddle to a Hartle–Hawking-like one. The same study isolates an “arc problem” associated with the lapse integration domain: for Robin conditions, the contour deformation can generate an additional arc contribution from $$\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)} = \sum_\sigma m_\sigma \int_{\mathcal J_\sigma} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)},$$1, which may dominate over the thimble contribution in some parameter regions (Chou et al., 2024).

A central unresolved technical problem in multivariable settings is the stable determination of intersection numbers. A recent development addresses this by solving the upward flow boundary-value problem with a multiple shooting method. For oscillatory integrals

$$\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)} = \sum_\sigma m_\sigma \int_{\mathcal J_\sigma} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)},$$2

the coefficients

$$\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)} = \sum_\sigma m_\sigma \int_{\mathcal J_\sigma} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)},$$3

are computed by propagating the unstable manifold from each saddle to the original integration cycle, including the orientation sign. The method uses normalized upward flow and has been tested for systems with up to $$\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)} = \sum_\sigma m_\sigma \int_{\mathcal J_\sigma} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)},$$4 variables; in discretized real-time double-well path integrals it identifies nontrivial complex saddles with $$\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)} = \sum_\sigma m_\sigma \int_{\mathcal J_\sigma} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)},$$5, $$\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)} = \sum_\sigma m_\sigma \int_{\mathcal J_\sigma} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)},$$6, or $$\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)} = \sum_\sigma m_\sigma \int_{\mathcal J_\sigma} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)},$$7, thereby exposing which complex trajectories actually contribute (Shoji et al., 7 Oct 2025).

6. Relations to complex Langevin, symmetry constraints, and resurgent viewpoints

The relation between thimble methods and complex Langevin is close but nontrivial. Both complexify the dynamical variables, but the generalized Lefschetz-thimble method deforms the contour geometrically, whereas complex Langevin samples a complexified distribution stochastically. A combined formulation applies complex Langevin to the real coordinates that parametrize a flowed contour $$\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)} = \sum_\sigma m_\sigma \int_{\mathcal J_\sigma} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)},$$8. In that framework, the partially phase-quenched version interpolates continuously between ordinary complex Langevin at $$\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)} = \sum_\sigma m_\sigma \int_{\mathcal J_\sigma} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)},$$9 and the original Lefschetz-thimble method as $$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$00. The construction clarifies that treating the Jacobian phase by reweighting is what makes the large-flow-time limit behave like genuine thimble integration (Nishimura et al., 2017).

The multi-thimble analyses of finite-density fermion models also explain why complex Langevin may fail. In the difficult Silver Blaze regime of the one-site model, the exact answer is a coherent sum of saddle contributions with complex coefficients, whereas the semiclassical complex-Langevin picture yields a positive-weight mixture of saddle values. That mismatch prevents correct reproduction of the required phase interference among saddles (Tanizaki et al., 2016).

Symmetry can constrain the thimble decomposition in a decisive way. For integrals satisfying

$$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$01

with an involutive linear map $$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$02, the corresponding anti-linear map

$$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$03

pairs saddles and thimbles so that the partition function remains manifestly real even when the saddles themselves are complex. In dense QCD mean-field models this is the $$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$04 symmetry, and the thimble decomposition can be organized into $$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$05-invariant saddles and conjugate thimble pairs. The resulting steepest-descent expansion preserves the reality of physical quantities order by order (Tanizaki et al., 2015).

Lefschetz-thimble decomposition also intersects with Borel resummation and Dyson–Schwinger truncation. In the zero-dimensional quartic model,

$$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$06

the perturbative saddle at $$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$07 and the nonperturbative saddles $$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$08 generate distinct asymptotic sectors. For $$\langle \mathcal O \rangle = \frac{\int_{\mathcal D} d\phi\, \mathcal O(\phi)\, e^{-S(\phi)}} {\int_{\mathcal D} d\phi\, e^{-S(\phi)}}.$$09, the perturbative saddle series is Borel summable, but it does not reconstruct the full integral; the full answer requires the nonperturbative thimbles as well. The corresponding conclusion for Dyson–Schwinger truncation is that one must use the large-order structure implied by the full thimble decomposition rather than forcing the highest correlator to vanish (Peng et al., 2024).

A related contemporary development uses additive weight regularizations to deform thimble structures in models with compact domains. The empirical criterion emerging from these studies is that complex Langevin converges correctly when the regularized system exhibits a single relevant compact thimble. Bias correction based on Dyson–Schwinger identities is then used to recover observables of the original theory. This does not provide a universal constructive prescription for lattice field theory, but it strengthens the view that thimble structure is a diagnostic for the success or failure of stochastic complexification methods (Boguslavski et al., 2024).

Lefschetz thimble decomposition is therefore best understood not as a single algorithm, but as a geometric framework for reorganizing complex integrals. Its exact content is homological, its computational leverage comes from constant-phase steepest-descent manifolds, and its physical significance is clearest in problems where interference among saddles encodes nonperturbative structure: Silver Blaze behavior at finite density, real-time tunneling, Lorentzian quantum cosmology, and trans-series completions of perturbation theory.