We study the chiral Ising, the chiral XY, and the chiral Heisenberg models at four-loop order with the perturbative renormalization group in 4 − ϵ dimensions and compute critical exponents for the GrossNeveu-Yukawa fixed points to order Oðϵ 4 Þ. Further, we provide Padé estimates for the correlation length exponent, the boson and fermion anomalous dimension, as well as the leading correction to scaling exponent in 2 þ 1 dimensions. We also confirm the emergence of supersymmetric field theories at four loops for the chiral Ising and the chiral XY models with N ¼ 1=4 and N ¼ 1=2 fermions, respectively. Furthermore, applications of our results relevant to various quantum transitions in the context of Dirac and Weyl semimetals are discussed, including interaction-induced transitions in graphene and surface states of topological insulators.