  We perform a numerical bootstrap study of scalar operators in the critical 3d
Gross-Neveu-Yukawa models, a family of conformal field theories containing N
Majorana fermions in the fundamental representation of an O(N) global symmetry.
We compute rigorous bounds on the scaling dimensions of the next-to-lowest
parity-even and parity-odd singlet scalars at N = 2, 4, and 8. All of these
dimensions have lower bounds greater than 3, implying that there are only two
relevant singlet scalars and placing constraints on the RG flow structure of
these theories.
