Dear all, I'm trying to simulate the photonic bands of a 2D structure made of blocks (i.e. squares) arranged in a triangular lattice. I'm using the python interface and the latest release of MPB.
I'm a pretty new user so I've started by modifying the mpb_tri_rods.py file from the Data Anlysis tutorial (see https://mpb.readthedocs.io/en/latest/Python_Data_Analysis_Tutorial/). That tutorial shows that the dielectric function returned by get_epsilon() is distorted due to the fact that the cell of our lattice is no longer a square, and uses the convert() function from a MPBData instance (with rectify=True) to correct that. However, now that I'm working on squares, that doesn't happen: if I plot the raw dielectric costant data (e.g. by matplotlib's matshow()) I get a perfect square, while rectify=True distorts the squares into parallelograms. Is this the correct behaviour? Should I see (before any conversion is made) a square or a parallelogram? I'm worried that I'm running my simulations with a wrong geometry. I've already tried manually defining the same structure using a prism insted of a block, but the result is the same. Best regards, Andrea Barzaghi
_______________________________________________ mpb-discuss mailing list mpb-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
