We derive the modified Young's equation for the equilibrium dihedral angle at the triple junction of the grain boundary groove by taking into account the discrete structure of the low angle grain boundary. For low angle grain boundaries, the geometric relation that the misorientation of the bicrystal is inversely proportional to the dislocation spacing naturally gives rise to the variation in the misorientation when the grain boundary length changes (holding the number of dislocations constant). The fact that the grain boundary energy increases as the grain boundary length decreases due to a smaller dislocation spacing leads to a larger dihedral angle compared to that of the classical theory. Two atomistic continuum modelling tools, namely the phase field crystal model and the amplitude equations, are used to simulate the equilibrium dihedral angle. The numerical results are in quantitatively good agreement with the derived modified Young's equation. Furthermore, the amplitude equations are employed to investigate the kinetics of the grain boundary grooving. The time-independent groove shape as predicted by Mullins is observed for bicrystal with high misorientation, and the groove width scales with time approximately as a power law t1/4. For bicrystals with small misorientations, the groove root exhibits a stick-slip motion due to volatile dislocation motion toward the liquid phase as it is close to the triple junction.