Hello,

How to compute the energy of a field in a supercell? Take e.g. MPB's tutorial script for triangular rods: https://github.com/stevengj/mpb/blob/master/examples/tri-rods.ctl. I can calculate the energy in one rod by invoking:

((output-dpwr-in-objects output-dpwr 0.001 (make cylinder (center 0 0 0)
                             (radius 0.2) (height infinity)
                       (material (make dielectric (epsilon 12)))
                       )) 4)

But how to do it for all rods in the supercell? Making a list doesn't work, as (output-dpwr-in-objects) doesn't work with this type of argument:

ERROR: wrong type for argument type (list . geometric-object) in compute-energy-in-object-list

Regards
Adam

_______________________________________________
mpb-discuss mailing list
mpb-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss

Reply via email to