In the simulations, a train of 3.5 ps long Gaussian pulses approximating Fig. 1(a) is launched through 5mm of Te in the phase matched orientation. The most intense pulse is set at 20 GW/cm\(^2\) , and noise is modeled in the time domain with magnitude proportional to the pulse intensity in each discretization bin [1]}. The Raman resopnse \(h_R(t)\)  was constructed by fitting a response function [2]} to the known Raman spectrum [3]} in the frequency domain.
