Solve Polarization under dynamic scenario

TDGL equation is the slow-limit version of Newton's 2nd law. We should revert it back to the original Newton's with P defined as a function of displacements: https://ieeexplore.ieee.org/stamp/stamp.jsp?tp=&arnumber=7990241