In this article, molecular dynamics based simulations were carried out to study the fracture toughness of single crystals of niobium (Nb) and zirconium (Zr). Separate set of simulations were performed with different orientations of crack plane in Nb and Zr. In each case, efforts were made to align the crack front with principal planes of corresponding crystal structure, that was bcc for Nb and hcp for Zr. Deformation in single crystal of Nb and Zr was governed either by twinning or emergence of dislocations from the crack tip and surfaces. The deformation mechanism in single crystal of Nb and Zr also helps in governing the overall toughness of the material, while deforming via twinning leads to higher change in the crack length, whereas dislocations emerging from the crack tip blunts the opening and improves the toughness.