Photonic Crystal Fiber Custom Wave Propagation

[kuroshfrz]

Hello All,

I am trying to simulate wave propagation of a Gaussian wave (and a custom input light) through a Photonic Crystal Fiber.
I am struggling to find a way to plot the beam profile (Intensity in the xy plane). I keep getting an error while trying to extract the output beam profile, stating that the corresponding functions from the Processing.jl gives modes only in multimode simulations.
I would appreciate it if you could help me do these tasks for my simulations.

[chrisbrahms]

Hi,

From what you're saying, it sounds like you are running a single-mode (mode-averaged) simulation. Is that correct? If so, the error you're getting is expected. There is little point in plotting the beam profile of a single-mode simulation, because it's trivial--it will always look like the single mode you're propagating. If you want to do it anyway, there is probably a way to hack it, though.

To help you further it would be much easier if you could post an example of the code you're running and which is giving you the error. Otherwise I'm just guessing :)

[kuroshfrz]

Hello Chris,

Thanks for your reply. That's true. I'm using the prop_gnlse to propagate light through the PCF. What I want to do is to extract the time-domain electric field in the x and y directions at some zslice values. Here is a sample code of what I'm currently working on. It's technically the same as the example for redoing the Fig.3 of Dudley et. al, RMP 78 1135 (2006).
I was wondering if Luna is also capable of simulating the nonlinear polarization rotation.

γ = 0.11
flength = 15e-2
βs = [0.0, 0.0, -1.1830e-26, 8.1038e-41, -9.5205e-56, 2.0737e-70, -5.3943e-85, 1.3486e-99, -2.5495e-114, 3.0524e-129, -1.7140e-144]
τfwhm = 50e-15
λ0 = 835e-9
power = 10000.0
output = prop_gnlse(γ, flength, βs; λ0, τfwhm, power, pulseshape=:sech, λlims=(400e-9, 2400e-9), trange=12.5e-12, ramanmodel=:sdo, τ1=12.2e-15, τ2=32e-15)

Best regards