The X-R ay D iffraction SIMULATOR package defines polycrystals as a mesh of tetrahedral single crystals and simulates diffraction as collected by a 2D discretized detector array while the sample is rocked around an arbitrary rotation axis. The full journal paper associated to the release of this code can be found here:
xrd_simulator: 3D X-ray diffraction simulation software supporting 3D polycrystalline microstructure morphology descriptions Henningsson, A. & Hall, S. A. (2023). J. Appl. Cryst. 56, 282-292. https://doi.org/10.1107/S1600576722011001
xrd_simulator was originally developed with the hope to answer questions about measurement optimization in scanning x-ray diffraction experiments. However, xrd_simulator can simulate a wide range of experimental diffraction setups. The essential idea is that the sample and beam topology can be arbitrarily specified, and their interaction simulated as the sample is rocked. This means that standard "non-powder" experiments such as scanning-3dxrd and full-field 3dxrd (or HEDM if you like) can be simulated as well as more advanced measurement sequences such as helical scans for instance. It is also possible to simulate powder like scenarios using orientation density functions as input.
Before reading all the boring documentation (which is hosted here) let's dive into some end to end examples to get us started on a good flavour.
The xrd_simulator is built around four python objects which reflect a diffraction experiment:
- A beam of x-rays (using the
xrd_simulator.beammodule)- A 2D area detector (using the
xrd_simulator.detectormodule)- A 3D polycrystal sample (using the
xrd_simulator.polycrystalmodule)- A rigid body sample motion (using the
xrd_simulator.motionmodule)
Once these objects are defined it is possible to let the detector collect scattering of the polycrystal as the sample undergoes the prescribed rigid body motion while being illuminated by the xray beam.
Let's go ahead and build ourselves some x-rays:
import numpy as np from xrd_simulator.beam import Beam # The beam of xrays is represented as a convex polyhedron # We specify the vertices in a numpy array. beam_vertices = np.array([ [-1e6, -500., -500.], [-1e6, 500., -500.], [-1e6, 500., 500.], [-1e6, -500., 500.], [1e6, -500., -500.], [1e6, 500., -500.], [1e6, 500., 500.], [1e6, -500., 500.]]) beam = Beam( beam_vertices, xray_propagation_direction=np.array([1., 0., 0.]), wavelength=0.28523, polarization_vector=np.array([0., 1., 0.]))
We will also need to define a detector:
from xrd_simulator.detector import Detector # The detector plane is defined by it's corner coordinates det_corner_0,det_corner_1,det_corner_2 detector = Detector(pixel_size_z=75.0, pixel_size_y=55.0, det_corner_0=np.array([142938.3, -38400., -38400.]), det_corner_1=np.array([142938.3, 38400., -38400.]), det_corner_2=np.array([142938.3, -38400., 38400.]))
Next we go ahead and produce a sample, to do this we need to first define a mesh that describes the topology of the sample, in this example we make the sample shaped as a ball:
from xrd_simulator.mesh import TetraMesh # xrd_simulator supports several ways to generate a mesh, here we # generate meshed solid sphere using a level set. mesh = TetraMesh.generate_mesh_from_levelset( level_set=lambda x: np.linalg.norm(x) - 768.0, bounding_radius=769.0, max_cell_circumradius=450.)
Every element in the sample is composed of some material, or "phase", we define the present phases in a list of xrd_simulator.phase.Phase objects, in this example only a single phase is present:
from xrd_simulator.phase import Phase quartz = Phase(unit_cell=[4.926, 4.926, 5.4189, 90., 90., 120.], sgname='P3221', # (Quartz) path_to_cif_file=None # phases can be defined from crystalographic information files )
The polycrystal sample can now be created. In this example the crystal elements have random orientations and the strain is uniformly zero in the sample:
from scipy.spatial.transform import Rotation as R from xrd_simulator.polycrystal import Polycrystal orientation = R.random(mesh.number_of_elements).as_matrix() polycrystal = Polycrystal(mesh, orientation, strain=np.zeros((3, 3)), phases=quartz, element_phase_map=None)
We may save the polycrystal to disc by using the builtin save() command as
polycrystal.save('my_polycrystal', save_mesh_as_xdmf=True)
We can visualize the sample by loading the .xdmf file into your favorite 3D rendering program. In paraview the sampled colored by one of its Euler angles looks like this:
We can now define some motion of the sample over which to integrate the diffraction signal:
from xrd_simulator.motion import RigidBodyMotion motion = RigidBodyMotion(rotation_axis=np.array([0, 1/np.sqrt(2), -1/np.sqrt(2)]), rotation_angle=np.radians(1.0), translation=np.array([123, -153.3, 3.42]))
Now that we have an experimental setup we may collect diffraction by letting the beam and detector interact with the sample:
polycrystal.diffract(beam, detector, motion) diffraction_pattern = detector.render(frames_to_render=0, lorentz=False, polarization=False, structure_factor=False, method="project")
The resulting rendered detector frame will look something like the below. Note that the positions of the diffraction spots may vary as the crystal orientations were randomly generated!:
import matplotlib.pyplot as plt fig,ax = plt.subplots(1,1) ax.imshow(diffraction_pattern, cmap='gray') plt.show()
To compute several frames simply change the motion and collect the diffraction again. The sample may be moved before each computation using the same or another motion.
polycrystal.transform(motion, time=1.0) polycrystal.diffract(beam, detector, motion)
Many more options for experimental setups and intensity rendering exist, have fun experimenting! The above example code can be found as a single .py file here.
xrd_simulator is distributed on the conda-forge channel and the preferred way to install the xrd_simulator package is via Anaconda:
conda create -n xrd_simulator
conda activate xrd_simulator
conda install -c conda-forge xrd_simulatorThis is meant to work across OS-systems and requires an Anaconda installation.
(The conda-forge feedstock of xrd_simulator can be found here.)
To install with anaconda on windows you must make sure that external dependencies of pygalmesh are preinstalled on your system. Documentation on installing these package can be found elsewhere.
Pip installation is possible, however, external dependencies of pygalmesh must the be preinstalled on your system. Installation of these will be OS dependent and documentation can be found elsewhere.:
pip install xrd-simulatorNaturally one may also install from the sources:
git clone https://github.com/FABLE-3DXRD/xrd_simulator.git
cd xrd_simulator
python setup.py installThis will then again require the pygalmesh dependencies to be resolved beforehand.
xrd_simulator makes good use of xfab and pygalmesh. The source code of these repos can be found here:
If you feel that xrd_simulator was helpful in your research we would love for you to cite us.
xrd_simulator: 3D X-ray diffraction simulation software supporting 3D polycrystalline microstructure morphology descriptions Henningsson, A. & Hall, S. A. (2023). J. Appl. Cryst. 56, 282-292. https://doi.org/10.1107/S1600576722011001