Thank you for your continuing help! According to your line of thought, I checked the simulation detail as follow:
1. I calculate the pulse energy by integrating the poynting vector from 2D time monitor, which I set at 0.1 and 25 micrometers away from the TE mode source as input and output. Here’s the script I use to calculate the pulse energy:
p_in = getdata("0.1μm 2D Y-normal time monitor","power");
t_in = getdata("0.1μm 2D Y-normal time monitor","t");
power_in = integrate(p_in,2,t_in);
2. When I change the point for field profile from 10-20, the frequency range doesn’t change: for the mode source, the 200fs pulse centered at 1550nm has FWHM from 1540 to 1560nm. And for the monitor, I set the frequency range from 1400 to 1700nm in “set global monitor settings” in the DFT monitor(I didn’t find that in time monitor).
3. The source amplitude is set to 3e9, and the index change is about 0.291.(I calculated by the following formula in ”FDTD Analysis of Two-Photon Absorption and Free-Carrier Absorption in Si High-Index-Contrast Waveguides”)
This index change is pretty larger than 0.01, but according to the theory model of “kerrnonlinear”, all I need is to entry only one permittivity value(which is the square of index, 3.48), and this model seems to have nothing to do with wavelength?
4. I also tried to set lower source amplitude or lower Chi3, unsurprisingly the loss is gradully disappearing, but the broadened spectrum is also returning to its original source spectrum, making this simualtion seems to be a linear one.
Thank you for your patient guidance, here are my entire script(plugin parameters are listed before), could you please take a glance in case I made some stupid mistake? I ‘m really eager to know why this happened.