The functionalization of nanoparticles (NPs) with DNA has proven to be an effective strategy for self-assembly of NPs into superlattices with a broad range of lattice symmetries. By combining this strategy with the DNA origami approach, the possible lattice structures have been expanded to include the cubic diamond lattice. This symmetry is of particular interest, both due to the inherent synthesis challenges, as well as the potential valuable optical properties, including a complete band-gap. Using these lattices in functional devices requires a robust and stable lattice. Here, we use molecular simulations to investigate how NP size and DNA stiffness affect the structure, stability, and crystallite shape of NP superlattices with diamond symmetry. We use the Wulff construction method to predict the equilibrium crystallite shape of the cubic diamond lattice. We find that, due to reorientation of surface particles, it is possible to create bonds at the surface with dangling DNA links on the interior, thereby reducing surface energy. Consequently, the crystallite shape depends on the degree to which such surface reorientation is possible, which is sensitive to DNA stiffness. Further, we determine dependence of the lattice stability on NP size and DNA stiffness by evaluating relative Gibbs free energy. We find that the free energy is dominated by the entropic component. Increasing NP size or DNA stiffness increases free energy, and thus decreases the relative stability of lattices. On the other hand, increasing DNA stiffness results in a more precisely defined lattice structure. Thus, there is a trade off between structure and stability of the lattice. Our findings should assist experimental design for controlling lattice stability and crystallite shape.