In this work, molecular dynamics (MD) is combined with continuum non-local crystal plasticity (CP) model[1] to investigate the deformation behaviour of the Inconel 718 (IN718) superalloy, which is mainly strengthened by coherent Ni3Nb γ" precipitates. The MD simulations reveal that dislocation shearing of the γ" precipitate occurs for the precipitate size and volume fraction considered. It was also found that interface dislocations develop and interact with antiphase boundaries (APBs) created by dislocations that move through the precipitates, and that γ" shearing was found to happen in four well-defined octahedral slip systems. A novel non-local CP modeling framework is proposed to account for the mechanical behaviour of individual IN718 grains as well as that of IN718’s individual phases. The CP model for the IN718 grains incorporates explicitly the strengthening effect of the mean γ" precipitate size and volume fraction. A multiscale approach is used to calibrate the constitutive parameters of the CP model for the different phases. Results from the CP simulations at the level of the individual γ" precipitates reveal that the three γ" variants behave differently under uniaxial loading, leading to strong strain and stress heterogeneities in the matrix phase. The size and volume fraction of the γ" precipitates were also found to affect the subsequent strain hardening behaviour at this level due to the development of geometrically necessary dislocations (GNDs) near the precipitate interfaces. An investigation into the effect of the γ" precipitate size on the overall mechanical behaviour of the superalloy is also performed at the level of IN718 grains.