We investigate the dynamics of photonic crystal (PC) molecules with Kerr nonlinearity by numerical simulations in a pump-probe scheme based on the finite-difference time-domain technique. The constitutional PC atoms of the nonlinear PC molecules are intentionally chosen to be multimode so that the p