Hi

To perform such an act ...

The following example can be modeled.

Which is a Lumerical example.

https://support.lumerical.com/hc/en-us/articles/360034395254-Calculating-absorbed-optical-power-Higher-accuracy-method-with-multiple-materials

Motahari

Dear Mohamad Reza Motahari,

Great job in figuring it out yourself. There is no direct way to do it, but I had also achieved what you needed through a script file. The trick is to get the absorption over the entire area that you have highlighted and delete the absorption values in the area that you do not need to take into consideration.

My approach was for nanostructures and so I executed additional processes to eliminate the numerical artifacts that occur due to the meshing technique in Lumerical FDTD solver.

The code below can be used for any monitor output as long as the input variable (in this case, Pabs) is changed accordingly.

My code:

## Setting absorption to areas only inside of the active layer

#Initialization

material="PPDT2FBTPC71BM"; # material

min_f=c/lam_nm(end,1)*e-9; # source min frequency

max_f=c/lam_nm(1,1)*e-9; # source max frequency

# Since the structure is built in the y-axis, let's get the y-axis data points

# (depends on meshing and so it is automated)

activeLayer_y = getdata("index_act_L","index_y",1);

# get FDTD fit of experimental data

n_fdtd=real(getfdtdindex(material,max_f,min_f,max_f));

## rounding to lower integer after multiplying by 100

# This will give an accuracy up to 2 decimal places

n_fdtd = floor(100*n_fdtd);

temp_v1 = floor(100*real(activeLayer_y(:,:,1,1)));

# Finding the matrix indices which should be made 0;

temp_v2 = abs(temp_v1 - n_fdtd) < 2;

## Remove numerical artifacts in FDTD for nanostructures

# This step requires matlab integration

matlabput(temp_v2);

matlab("temp_v3 = double(bwareaopen(temp_v2, 2));");

matlabget(temp_v3);

# Pabs is the output of the absorption monitor

removeNonActiveLayerPhotons = Pabs;

removeNonActiveLayerPhotons(:,:,1) = temp_v3;

Pabs = Pabs * removeNonActiveLayerPhotons;

temp_v1 = 0;

temp_v2 = 0;

temp_v3 = 0;

removeNonActiveLayerPhotons = 0;

# end absorption calculation