LAUSR

We analyze emergent quantum multicriticality for strongly interacting, massless Dirac fermions in two spatial dimensions ($d=2$) within the framework of Gross-Neveu-Yukawa models, by considering the competing order parameters that give rise to fully gapped (insulating or superconducting) ground states. We focus only on those competing orders which can be rotated into each other by generators of an exact or emergent chiral symmetry of massless Dirac fermions, and break $\text{O}({S}_{1})$ and $\text{O}({S}_{2})$ symmetries in the ordered phase. Performing a renormalization-group analysis by using the $\ensuremath{\epsilon}=(3\ensuremath{-}d)$ expansion scheme, we show that all the coupling constants in the critical hyperplane flow toward a new attractive fixed point, supporting an enlarged $\text{O}({S}_{1}+{S}_{2})$ chiral symmetry. Such a fixed point acts as an exotic quantum multicritical point (MCP), governing the continuous semimetal-insulator as well as insulator-insulator (for example, antiferromagnet to valence bond solid) quantum phase transitions. In comparison with the lower symmetric semimetal-insulator quantum critical points, possessing either $\text{O}({S}_{1})$ or $\text{O}({S}_{2})$ chiral symmetry, the MCP displays enhanced correlation length exponents, and anomalous scaling dimensions for both fermionic and bosonic fields. We discuss the scaling properties of the ratio of bosonic and fermionic masses, and the increased dc resistivity at the MCP. By computing the scaling dimensions of different local fermion bilinears in the particle-hole channel, we establish that most of the four fermion operators or generalized density-density correlation functions display faster power-law decays at the MCP compared to the free fermion and lower symmetric itinerant quantum critical points. Possible generalization of this scenario to higher-dimensional Dirac fermions is also outlined.