X-Ray Diffraction (XRD)

This example shows how to use pymatgen to quickly estimate a powder diffractogram.

Note

You can install pymatgen into the AMS Python environment by installing the m3gnet package in AMS2023 or later.

To follow along, either

Initial imports

from scm.plams import *
try:
    from ase import Atoms
    from pymatgen.core.structure import Structure
    from pymatgen.analysis.diffraction.xrd import XRDCalculator
except ImportError as e:
    print("You need ASE and pymatgen installed in the AMS python environment to run this example. Install the package for m3gnet to do this.")
    print(e)
    exit(1)

Create ASE atoms object for BaTiO3

at = Atoms(
    symbols=['Ba', 'Ti', 'O', 'O', 'O',],
    scaled_positions=[
        [0., 0., 0.,],
        [0.5, 0.5, 0.5],
        [0.0, 0.0, 0.5],
        [0.0, 0.5, 0.0],
        [0.5, 0.0, 0.0]
    ],
    cell=[4.01, 4.01, 4.01],
    pbc=(True, True, True),
)
plot_molecule(at, rotation='-5x,5y,0z') # show in Jupyter notebook
../../_images/xrd_3_0.png

Save ASE Atoms to .cif format

fname = "batio3.cif"
at.write(fname)

Load .cif in pymatgen and calculate XRD

structure = Structure.from_file(fname)
xrd_calc = XRDCalculator()
xrd_calc.show_plot(structure)
../../_images/xrd_7_0.png 
pattern = xrd_calc.get_pattern(structure)
print("2*Theta Intensity hkl d_hkl(angstrom)")
for two_theta, intensity, hkls, d_hkl in zip(pattern.x, pattern.y, pattern.hkls, pattern.d_hkls):
    hkl_tuples = [hkl["hkl"] for hkl in hkls]
    for hkl in hkl_tuples:
        label = ", ".join(map(str, hkl))
        print(f'{two_theta:.2f} {intensity:.2f} {hkl} {d_hkl:.3f}')

2*Theta Intensity hkl d_hkl(angstrom)
22.17 46.84 (1, 0, 0) 4.010
31.55 100.00 (1, 1, 0) 2.835
38.90 1.83 (1, 1, 1) 2.315
45.23 34.58 (2, 0, 0) 2.005
50.92 19.69 (2, 1, 0) 1.793
56.19 38.27 (2, 1, 1) 1.637
65.88 20.48 (2, 2, 0) 1.418
70.44 9.47 (2, 2, 1) 1.337
70.44 9.47 (3, 0, 0) 1.337
74.88 16.60 (3, 1, 0) 1.268
79.23 1.68 (3, 1, 1) 1.209
83.51 6.82 (2, 2, 2) 1.158
87.76 4.44 (3, 2, 0) 1.112

Complete Python code

#!/usr/bin/env amspython
# coding: utf-8

# ## Initial imports

from scm.plams import *
try:
    from ase import Atoms
    from pymatgen.core.structure import Structure
    from pymatgen.analysis.diffraction.xrd import XRDCalculator
except ImportError as e:
    print("You need ASE and pymatgen installed in the AMS python environment to run this example. Install the package for m3gnet to do this.")
    print(e)
    exit(1)


# ## Create ASE atoms object for BaTiO3

at = Atoms(
    symbols=['Ba', 'Ti', 'O', 'O', 'O',], 
    scaled_positions=[
        [0., 0., 0.,], 
        [0.5, 0.5, 0.5], 
        [0.0, 0.0, 0.5], 
        [0.0, 0.5, 0.0], 
        [0.5, 0.0, 0.0]
    ], 
    cell=[4.01, 4.01, 4.01],
    pbc=(True, True, True),
)
plot_molecule(at, rotation='-5x,5y,0z') # show in Jupyter notebook


# ## Save ASE Atoms to .cif format

fname = "batio3.cif"
at.write(fname)


# ## Load .cif in pymatgen and calculate XRD

structure = Structure.from_file(fname)
xrd_calc = XRDCalculator()
xrd_calc.show_plot(structure)


pattern = xrd_calc.get_pattern(structure)
print("2*Theta Intensity hkl d_hkl(angstrom)")
for two_theta, intensity, hkls, d_hkl in zip(pattern.x, pattern.y, pattern.hkls, pattern.d_hkls):
    hkl_tuples = [hkl["hkl"] for hkl in hkls]
    for hkl in hkl_tuples:
        label = ", ".join(map(str, hkl))  
        print(f'{two_theta:.2f} {intensity:.2f} {hkl} {d_hkl:.3f}')