Relativistic Dirac-Coulomb ͑DC͒ Hartree-Fock calculations are employed to obtain the analytic electric field gradient ͑EFG͒ on the antimony nucleus in the SbN, SbP, SbF, and SbCl molecules. The electronic correlation contribution to the EFGs is included with the DC-CCSD͑T͒ and DC-CCSD-T approaches, also in the four-component framework, using a finite-difference method. The total EFG results, along with the experimental nuclear quadrupole coupling constants from microwave spectroscopy, allow to derive the nuclear quadrupole moments of 121 Sb and 123 Sb, respectively, as −543͑11͒ and −692͑14͒ mb.