Calculate XRD pattern for 2D material

Is there to calculate xrd pattern of single or Double layer 2D material using pymatgen or any other software. I am working on Li intercalated Phosphorene. I want to calculate XRD for different Li concentration.

Thanks in Advance

Please post pymatgen questions on the pymatgen discourse:

Look at the XRDCalculator in pymatgen: