Herbert Eßl, Stefan Rohshap, Marcel Gievers, Markus Wallerberger, Alessandro Toschi, Anna Kauch
We systematically analyze the stability of the iterative solution of the parquet equations by studying the spectrum of the Jacobian associated with the commonly used damped fixed-point iteration procedure. In this context, we provide an explicit criterion that determines when the physical fixed point of the parquet iteration becomes unstable. Importantly, we demonstrate that misleading convergence issues, observed in parquet calculation at intermediate-to-high interaction values, are not restricted to parameter regions where the two-particle irreducible vertex diverges, but can also arise in absence of vertex divergences. Hence, the misleading convergence issues of parquet-based algorithms are not directly caused by the crossings of two solutions of the (multivalued) Luttinger-Ward functional, that are associated with vertex divergences. Building on these insights, we introduce a controlled stabilization strategy that allows the convergence to the physical solution in the instability regimes. We apply this procedure to the zero-point model and the Hubbard model in the atomic limit, where we successfully stabilize the physical solution deep in the non-perturbative regime, even across multiple divergence lines.
This paper provides a systematic analysis of why iterative parquet equation solvers fail to converge at intermediate-to-strong coupling, and offers a principled stabilization strategy. The core contributions are threefold:
First, the authors derive an explicit stability criterion for the physical fixed point of the damped fixed-point iteration by computing the full Jacobian of the parquet iteration map. They show that eigenvalues of the matrix Π with negative real parts render the physical solution unreachable by any damping parameter p ∈ (0,1].
Second, they demonstrate that convergence failure is *not* exclusively caused by divergences of the two-particle irreducible vertex Γ (which are connected to crossings of branches of the Luttinger-Ward functional). Instead, instabilities can arise through zero-crossings of real eigenvalues or through complex conjugate pairs acquiring negative real parts — mechanisms that occur even in absence of vertex divergences. This is a conceptually important clarification that disentangles two phenomena previously assumed to be linked.
Third, they introduce a stabilization method that flips the sign of the damping along unstable eigendirections of the Jacobian, making the physical fixed point locally stable. This is applied successfully to the zero-point model and the Hubbard atom deep into non-perturbative regimes.
The methodology is thorough and mathematically well-grounded. The paper builds on the general framework of Ref. [1] (by overlapping authors) and specializes it to the parquet context. Key strengths include:
One limitation is that the Jacobian is evaluated at the *known* physical fixed point for these proof-of-principle calculations. The authors acknowledge this and discuss strategies for realistic applications where the fixed point is unknown (Sec. 5), but this remains undemonstrated.
The parquet formalism is central to the study of strongly correlated electrons, particularly within the dynamical vertex approximation (DΓA) framework. Convergence problems at strong coupling have been a persistent practical bottleneck. This work:
The impact extends to adjacent self-consistent diagrammatic methods (GW, FLEX, etc.) that face analogous stability issues, and the companion paper [28] apparently explores this direction.
This work addresses a well-recognized bottleneck in computational many-body physics. The parquet formalism has seen renewed interest due to advances in vertex compression (quantics tensor trains, sparse intermediate representations) and the single-boson exchange formalism. Several groups are actively pushing parquet-based methods to realistic lattice problems. The convergence barrier at intermediate coupling has been a major obstacle, making this contribution timely.
The connection to the broader question of Luttinger-Ward functional multivaluedness — a topic of intense debate since the 2015 Kozik-Ferrero-Georges paper — adds fundamental relevance.
The paper is well-structured but very long (53 pages). The pedagogical examples in Sec. 3.1 effectively communicate the key ideas. The discussion in Sec. 5 about routes to realistic implementation is honest about remaining challenges, which strengthens credibility. The suggestion to combine stabilization with quantics tensor train compression or finite-difference parquet methods points to concrete future research directions.
Generated Jun 4, 2026
Paper 1 explores domain walls in altermagnets, a rapidly growing and highly topical area in condensed matter physics. Its findings on spin excitation localization and magnon propagation have direct implications for next-generation spintronics and magnonics devices. Paper 2, while methodologically rigorous and valuable for computational many-body theorists, focuses on algorithmic stability for solving parquet equations. Paper 1's combination of novel material physics with tangible device applications gives it a broader and more immediate potential scientific impact across both theoretical and experimental domains.
Paper 2 likely has higher scientific impact due to broad methodological relevance: it provides a general stability criterion and a controlled stabilization strategy for parquet solvers, addressing a widespread computational bottleneck in strongly correlated-electron theory. This can improve reliability of many-body calculations across models and materials, with applications spanning condensed matter, quantum chemistry, and numerical method development. Paper 1 is novel and experimentally compelling but is narrower in scope (specific perovskites and magnetic mechanism), so its cross-field and tool-building impact is more limited.
Paper 1 addresses a fundamental computational challenge in many-body physics—the stability of parquet equation solvers—with broad implications for correlated electron calculations across condensed matter and beyond. It provides both a diagnostic criterion and a practical stabilization strategy applicable to non-perturbative regimes, which could enable new physics results from parquet-based methods. Paper 2, while experimentally solid, corrects a previous structural determination of a specific kagome magnet, representing a more incremental and narrowly focused contribution with less methodological breadth.
Paper 2 addresses the highly active field of moiré materials and fractional Chern insulators. By exploring how anisotropy drives quantum phase transitions, it provides testable predictions for experimentalists using heterostrain. Its relevance to recent zero-field FCI discoveries gives it broader interest and higher potential impact in condensed matter physics compared to the narrower, methodological focus on algorithmic convergence in Paper 1.
Paper 1 addresses a high-profile topic in spintronics—p-wave altermagnetism—by providing the first direct experimental evidence that strong electronic correlations can quench predicted nonrelativistic spin splittings. This challenges a widely pursued theoretical paradigm (spin-space group symmetry predictions) and establishes fundamental constraints on when geometric symmetry classifications apply, impacting a broad community in condensed matter physics and spintronics. Paper 2, while methodologically rigorous and valuable for the computational many-body physics community, addresses a more technical numerical convergence issue in parquet equations with a narrower audience and scope of impact.
Paper 2 has higher likely impact: it addresses a timely, widely interesting topic (altermagnetism in heavy-fermion/f-electron systems) with direct experimental evidence (high-resolution ARPES) plus theory benchmarking (DFT vs DFT+U) that clarifies why a key predicted signature is absent. This has clear real-world relevance for spintronics/materials design and broad implications across condensed-matter subfields (magnetism, correlated electrons, spectroscopy, electronic-structure methods). Paper 1 is novel and rigorous but is more specialized and primarily methodological within parquet/diagrammatic many-body computations.
Paper 2 addresses a fundamental methodological challenge in computational many-body physics (convergence of the parquet equations) and introduces a stabilization strategy. This algorithmic advancement has broad applicability across various strongly correlated systems, offering high utility for researchers. Paper 1 presents highly novel theoretical findings on topological phases in a specific model, but its impact is narrower compared to the broad methodological improvements of Paper 2.
Paper 1 likely has higher scientific impact: it demonstrates an experimentally accessible, electrically read out probe (spin Seebeck effect) for fractionalized excitations (emergent monopoles) in a 3D magnetic insulator, advancing both quantum magnetism and spintronics with clear real-world measurement and device-interface relevance. Its novelty and cross-field breadth (condensed matter, materials, spin caloritronics) are strong and timely. Paper 2 is methodologically rigorous and valuable for computational many-body physics, but its impact is narrower (algorithmic stabilization of parquet solvers) and more specialized.
Paper 1 addresses a fundamental and notoriously difficult problem in condensed matter physics—understanding Non-Fermi liquids—by introducing a highly novel dimensional-separation framework. Taming severe logarithmic divergences to allow controlled perturbative analysis represents a major theoretical breakthrough. While Paper 2 offers a valuable algorithmic stabilization for the parquet equations, its impact is largely restricted to computational many-body methodology. Paper 1 has a significantly broader scope, higher theoretical novelty, and the potential to reshape our understanding of quantum phase transitions and strongly correlated systems.
Paper 2 likely has higher impact because it delivers a concrete methodological advance: an explicit stability criterion for parquet fixed-point iterations plus a controlled stabilization strategy, demonstrated on canonical models. This directly improves reliability of a widely used many-body computational framework and is broadly applicable across correlated-electron physics and materials modeling. Its rigor (Jacobian spectrum analysis, tested across divergence lines) and timeliness (addressing known convergence pathologies) suggest strong uptake. Paper 1 is a valuable perspective summarizing thin-film advances but is primarily review-like and narrower in immediate transferable methodology.