search

boyuanliuoptics / DOS-calculation

Programs for density-of-states calculation in band theory, especially in photonic crystals.
MIT License
13 stars 6 forks source link

DOS in one direction!! #13

Open tarafabio opened 6 years ago

tarafabio commented 6 years ago

Hello, Thanks to the group for the code. I have tried to derive DOS of the 3Dexample in just one direction (here G-N). I changed the ctl code and remove the other symmetry points except for G and N. But the result is the same as when I have DOS for whole the BZ ( G-N-H...). Noted, the figure of Bandstructure is correct and shows the bandstructure belong to the G-N direction. Using these code, Is it possible to load DOS in some specific region of BZ?

Thanks, Regards

boyuanliuoptics commented 6 years ago

The symmetry points settings are only related to the band structure plot in the figure. If you want to calculate the DOS in one direction, you need to modify the part of defining "kpts" into one dimension. While at present the .m codes can only calculate the total DOS of 3D and 2D, it is not difficult to extend the DOS calculation codes (especially GGR code) into 1D. But I am not sure if DOS in one direction of a 3D structure has any physical meaning. If the specific region is a parallelepiped, then the DOS in this region could be calculated from the codes.