Problems of 'PAbs' in scuff-neq

[ZhihaoJia16]

Hi, I compute the fluxes of power transfer between two SiO2 spheres. http://homerreid.github.io/scuff-em-documentation/examples/SiO2Spheres/SiO2Spheres/

• Why the PAbs of 12 and 21 is different?
• PAbs of 12 is positive and PAbs of 21 is negtive in some angular frequencies, why is that?

In this example, I run:

#!/bin/bash
for N in 501 
do
  ARGS=""
  ARGS="${ARGS} --geometry SiO2Spheres_${N}.scuffgeo"
  ARGS="${ARGS} --OmegaFile OmegaFile"
  ARGS="${ARGS} --Power --ZForce"
  ARGS="${ARGS} --DSIRadius 5"
  ARGS="${ARGS} --DSIPoints 302"
  scuff-neq ${ARGS}
done

get

#1 transform tag
#2 omega 
#3 (sourceObject,destObject) 
#4 PAbs flux spectral density
#5 PRad flux spectral density
DEFAULT 1.00E-01 11 4.19E-06   -0.00000000e+00
DEFAULT 1.00E-01 12 -8.51E-08  -0.00000000e+00
DEFAULT 1.00E-01 21 4.26E-08   -0.00000000e+00
DEFAULT 1.00E-01 22 4.19E-06   -0.00000000e+00
...
DEFAULT 2.00E-01 11 1.44E-04   -0.00000000e+00
DEFAULT 2.00E-01 12 1.12E-07   -0.00000000e+00
DEFAULT 2.00E-01 21 -6.32E-07  -0.00000000e+00
DEFAULT 2.00E-01 22 1.44E-04   -0.00000000e+00
...

Thank you in advance.