Predicting 2D THz Spectra Due to Nonlinear Phononics with First-Principles Calculations