A symmetry-preserving treatment of a vector $$\times $$
×
vector contact interaction is used to compute spectra of ground-state $$J^P = 0^\pm , 1^\pm $$
J
P
=
0
±
,
1
±
$$(f{\bar{g}})$$
(
f
g
¯
)
mesons, their partner diquark correlations, and $$J^P=1/2^\pm , 3/2^\pm $$
J
P
=
1
/
2
±
,
3
/
2
±
(fgh) baryons, where $$f,g,h \in \{u,d,s,c,b\}$$
f
,
g
,
h
∈
{
u
,
d
,
s
,
c
,
b
}
. Results for the leptonic decay constants of all mesons are also obtained, including scalar and pseudovector states involving heavy quarks. The spectrum of baryons produced by this chiefly algebraic approach reproduces the 64 masses known empirically or computed using lattice-regularised quantum chromodynamics with an accuracy of 1.4(1.2)%. It also has the richness of states typical of constituent-quark models and predicts many baryon states that have not yet been observed. The study indicates that dynamical, nonpointlike diquark correlations play an important role in all baryons; and, typically, the lightest allowed diquark is the most important component of a baryon’s Faddeev amplitude.