Owing to its ability to incorporate Schmid’s law at each integration point, crystal plasticity has proven a powerful tool to simulate and predict the slip behavior at the grain level and the ensuing heterogeneous stress/strain localization and texture evolution at the macroscopic level. Unfortunately, notwithstanding substantial efforts during the last three decades, this remarkable capability has not been replicated for materials where twinning becomes a noticeable deformation mechanism, namely in the case of low-stacking fault energy cubic, orthorhombic, and hexagonal close-packed structures. The culprit lies in the widely adopted unphysical pseudo-slip approach for capturing twin formation. While the slip is diffuse, twinning is a localized event that occurs as a drastic burst of a confined number of partial twinning dislocations establishing an interface that pursues growth through a thread of perfect twinning dislocations in the sense of bicrystallography. Moreover, at earlier stages, twin nucleation may require atomic diffusion (Shuffling) and faceting, generally demanding higher stress levels not necessarily on the twin shear plane, while triaxiality at adequate sites might be needed or preferred such as lower grain boundary misorientations or other twin boundaries. Identifying a mathematical framework in the constitutive equations for capturing these twin formation sensitivities has been a daunting challenge for crystal plasticity modelers, which has stalled ameliorating the design of key hexagonal materials for futuristic climate change-related industries. This paper reviews existing approaches to incorporating twinning in crystal plasticity models, discusses their capabilities, addresses their limitations, and suggests prospective views to fill gaps. The incorporation of a new physics-based twin nucleation criterion in crystal plasticity models holds groundbreaking potential for substantial progress in the field of computational material science.