I am a beginner user of the Pymatgen package. In order to process the data from VASP DFT calculation software, I use Pymatgen to visualize the output band structure. When I do spin polarised band calculation, I cannot separate the band structure of spin up and spin down into 2 different subplots or graphs. Can anyone tell me how can I deal with this? Because it looks messy when the package automatically plot it into the same graph.
Here I provide an example with my own python scripts to realize your purpose rather than using Pymatgen (You can save data firstly with Pymatgen and plot with python). I assume that you can perform correctly the band and dos calculations with the VASP code. The example I pick up is the monolayer FM NiBr$_2$ and the final result is the following:
- Example GitHub link: https://github.com/Kohn-X/VASP_two_spins