We propose a novel scheme for combining efficiently the truncated-unity functional renormalization group (TUFRG) and the mean-field theory. It follows the method of Wang, Eberlein, and Metzner that uses only the two-particle channel-irreducible part of the vertex as an input for the mean-field treatment. In the TUFRG, the neglect of fluctuation effects from other channels in the symmetry-broken regime is represented by applying the random phase approximation (RPA) in each individual channel, below the divergence scale. Then the Bethe-Salpeter equation for the four-point vertex is translated into the RPA matrix equations for the bosonic propagators that relates the singular and irreducible singular modes of the propagators. The universal symmetries for the irreducible singular modes are obtained from the antisymmetry of Grassmann variables. The mean-field equation based on these modes is derived by the saddle-point approximation in the framework of the path-integral formalism. By using our scheme, the power of the TUFRG, as a diagrammatically unbiased tool for identifying the many-body instabilities, could be elevated to a quantitatively reasonable level, and its application would be extended to a quantitatively reasonable analysis of the coexisting orders. As an illustration, we employ this scheme to study the coexistence phase of the chiral superconductivity and the chiral spin-density wave, predicted near van Hove filling of the honeycomb lattice.

• Received 11 December 2023
• Revised 8 April 2024
• Accepted 19 April 2024

DOI:https://doi.org/10.1103/PhysRevB.109.205118