Dispersion Plot

[Manjushree Maity Res. Scholar, Physics, IIT(BHU)]

  I am trying to reproduce the dispersion results from a paper (attached below), but I consistently get  extra modes even though I’ve used the same parameters as reported. This issue isn’t limited to just this case—I’ve also tried reproducing results from one or two other papers, and the same problem occurs: one additional mode always appears. Could anyone suggest what might be causing this and how I can obtain results that match the papers?

Best Regards,

Manju

[Антон Луценко]

Hi,

It would probably be helpful if you linked the paper and attached your simulation script. But one option to isolate a given mode would be to plot mode profiles and then change your excitation profile to maximize excitation of the chosen mode. 

[Manjushree Maity Res. Scholar, Physics, IIT(BHU)]

Hello Sir, 

  Thank you for your response. I am currently attempting to reproduce the results of the attached paper. However, I am encountering two issues: first, an extra mode appears in the dispersion, and second, I am unable to generate the transmission spectra. Could you please advise on how to properly obtain the transmission spectra plot and provide guidance on what might be causing the unexpected mode in the dispersion?

  

OutputFormat = OVF2_TEXT
Nx := 600
Ny := 10
Nz := 1
setgridsize(Nx, Ny, Nz)
setcellsize(1e-9, 1e-9, 1e-9)
Msat = 1.6e6
alpha = 0.001
Aex = 6e-12
Ku1 = 2e6
anisU = vector(0, 0, 1)
Dind = 3.0e-3
m = uniform(0, 0, 1)
relax()
save(m)
B0 := 10e-3
fcut := 100e9
excite := rect(30e-9, 10e-9)
defRegion(15, excite)
B_ext.SetRegion(15, vector(B0*sinc(2*pi*fcut*t), 0, 0))
tableautosave(5e-12)
autosave(m, 5e-12)
run(30e-9)