-
|
Hi everybody, I'm looking to reproduce the results of Chapter 4, Figure 10 from Molding the Flow of Light (2nd ed) for a (1D) multilayer structure. The image below shows (left) the text results and (right) my attempt at reproducing these with MPB. The bands along k = (0, ky, 0) (blue lines) seem to agree well with the text. The issue I'm seeing is with bands along k = (0, ky, pi/a) (green lines). Perhaps my code defines k incorrectly in this case, but I'm not sure that's the issue. It's definitely possible I am defining something incorrectly as I am relatively new to MPB. I would appreciate any insight the community is able to offer. Thank you!
|
Beta Was this translation helpful? Give feedback.
Replies: 3 comments
-
|
pi/a should be 0.5 in MPB units, assuming your period is a=1 |
Beta Was this translation helpful? Give feedback.
-
|
Also, for that figure in the book, the structure is periodic in z (see figure 9 in chapter 4), but in your code it is periodic in x, so your k directions don't match relative to the structure. |
Beta Was this translation helpful? Give feedback.
-
|
Awesome, thank you! I had to modify the orientation so that periodicity is along x since I saw that MPB has default orientation definitions for the Run functions (and possibly others), but it seems to be working now. Thanks again!
|
Beta Was this translation helpful? Give feedback.
pi/a should be 0.5 in MPB units, assuming your period is a=1