In an embodiment, a method of time-domain simulation for simulating scattering spectra is described. The method may provide for computing spatial derivatives including computing spectral derivatives in some portion of the domain and finite difference derivatives in some other portion of the domain; forming an equation for non-reflecting boundary conditions; and computing a time-stepping scheme using a high-order unconditionally stable computation method.