  We investigate bulk structural properties of tetravalent associating
particles within the framework of classical density functional theory, building
upon Wertheim's thermodynamic perturbation theory. To this end, we calculate
density profiles within an effective test-particle geometry and compare to
radial distribution functions obtained from computer simulations. We
demonstrate that a modified version of the functional proposed by Yu and Wu [J.
Chem. Phys. 116, 7094 (2002)] based on fundamental measure theory for hard
spheres produces accurate results, although the functional does not satisfy the
exactly known low-density limit. However, at low temperatures where particles
start to form an amorphous tetrahedral network, quantitative differences
between simulations and theory emerge due to the absence of geometrical
informations regarding the patch arrangement in the latter. Indeed, here we
find that the theory fits better to simulations of the floating-bond model [J.
Chem. Phys. 127, 174501 (2007)], which exhibits a weaker tetrahedral order due
to more flexible bonds between particles. We also demonstrate that another
common density functional approach by Segura \textit{et al.} [Mol. Phys. 90,
759 (1997)] fails to capture fundamental structural properties.
