GenL: An extensible fitting program for Laue oscillations and whole pattern fitting

GenL is a flexible, extensible, and open-source MATLAB-based program that utilizes a genetic algorithm to simulate and fit X-ray reflectivity and diffraction data from epitaxial thin films, offering both source code and pre-compiled binary options for extracting structural parameters like strain profiles and crystal roughness.

The Gist

Imagine you have built a very tall, perfect tower out of microscopic Lego bricks. You want to know exactly how many bricks are in the tower, how smooth the walls are, and if any of the bricks are slightly squished or stretched. To find out, you shine a special flashlight (X-rays) at the tower. The light bounces off the layers, creating a complex pattern of ripples and waves, much like the ripples you see when you drop a stone in a pond.

This paper introduces a new computer tool called GenL (which stands for "Genetic Laue") that helps scientists decode these ripples.

Here is a simple breakdown of what the paper says:

1. The Problem: Reading the Ripples

When scientists look at high-quality thin films (like the Lego tower), they see something called Laue oscillations. These are the ripples in the light pattern.

• The good news: If the ripples are clear and strong, it means the film is very high quality, smooth, and uniform.
• The bad news: Reading these ripples by hand is incredibly hard. It's like trying to guess the exact recipe of a cake just by tasting a crumb. You need to know the thickness, the roughness of the layers, and if the material is under stress (strain).

2. The Solution: GenL (The "Smart Detective")

The authors created GenL, a program that acts like a super-smart detective. Instead of just guessing, it uses a method called a genetic algorithm.

• The Analogy: Imagine you are trying to tune a radio to find a specific station. You turn the dial a little bit, listen, and if it's not clear, you turn it a different way. GenL does this millions of times in seconds. It tries thousands of different "recipes" (different thicknesses, roughness levels, and strains) until it finds the one that perfectly matches the ripples seen in the experiment.

3. How It Works: Two Ways to Look at the Tower

GenL can analyze the data in two different ways, depending on how deep you want to look:

• The "Quick Sketch" (Kinematic Theory): This is like looking at the tower from a distance. It's fast and good for a general idea. It counts the layers and checks the spacing.
• The "Microscope" (Dynamic Theory): This is like looking at the tower through a powerful microscope. It calculates exactly how the light bounces off every single atom, including how the light gets absorbed or bent. This is slower but much more precise for complex situations.

4. What GenL Can Tell You

Once GenL matches the "recipe" to the data, it tells the scientists:

• How thick the film is: Exactly how many atomic layers are there?
• How rough the surface is: Are the layers smooth like glass, or bumpy like sandpaper?
• Is it under stress? Are the atoms being pulled apart (tension) or squished together (compression)?
• Interface details: What happens right where the film meets the material underneath it?

5. Real-World Examples in the Paper

The authors tested GenL on real materials to show it works:

• Tungsten on Sapphire: They looked at a tungsten film and found it had a roughness of about 4 atoms. GenL confirmed this matched what other tools found, proving it was accurate.
• Iron on a Special Crystal: They looked at an iron film that was changing its shape as it got thicker. GenL was able to spot that the iron started as one shape (bct) and then changed into another (bcc) as it grew. It even figured out that the "ripples" in the light were asymmetrical because of this change.
• Superlattices: They showed that GenL can simulate complex stacks of different materials (like a Fe/V superlattice), proving it can handle complicated Lego towers made of different colored bricks.

6. Who Can Use It?

• It's Free: The program is open-source (free for everyone to use and change).
• Easy to Use: It comes with a "Graphical User Interface" (GUI). Think of this as a user-friendly dashboard with buttons and sliders, so you don't have to be a computer programmer to use it.
• Flexible: If you are a programmer, you can also run it as a command-line tool to build your own custom models.

Summary

GenL is a new, free, and flexible computer program that helps scientists understand the hidden details of ultra-thin films. By using a "smart search" method to match theoretical models with real X-ray data, it reveals the thickness, roughness, and internal stress of these materials with high precision. It is designed to be an easy-to-use tool that can be expanded by users to study even more complex structures in the future.

Technical Summary

Technical Summary of GenL: An Extensible Fitting Program for Laue Oscillations and Whole Pattern Fitting

Problem Statement
High-quality epitaxial thin films are critical model systems, where the presence of Laue oscillations in X-ray diffraction (XRD) patterns serves as a key indicator of crystallinity, thickness uniformity, and interface smoothness. While the qualitative presence of these oscillations is well-established, quantitative analysis of their angular spacing, intensity decay, and asymmetry offers deeper structural insights regarding coherent thickness, strain profiles, and roughness.

Existing open-source software (e.g., SUPREX, GenX, CADEM, InteractiveXRDFit) and commercial packages often focus primarily on fitting superlattice peaks or X-ray reflectivity (sensitive to electron density profiles) rather than the detailed shape of Laue oscillations. Many existing tools lack the capability to simultaneously simulate and fit experimental diffraction patterns with specific attention to atomic ordering, strain profiles, and crystal roughness, or they require significant reconfiguration to handle high-angle scattering data. There is a identified need for a flexible, open-source program capable of both simulating and fitting these specific features.

Methodology
GenL is a MATLAB-based program designed to simulate and fit X-ray diffraction data from epitaxial thin films. It employs a modular approach utilizing two theoretical frameworks:

1. Kinematic Theory: The primary approach for fitting, calculating diffracted intensity based on the summation of scattering amplitudes from coherently scattering lattice planes. It incorporates:

   • Form Factors: Element-specific, Q-dependent form factors including anomalous dispersion and absorption corrections.
   • Physical Corrections: Polarization factors (accounting for monochromator position), absorption factors (angle-dependent attenuation), and Lorentz factors (instrumental resolution and mosaicity).
   • Structural Parameters: Interplanar spacings, number of coherently scattering planes, and Debye-Waller factors.
   • Defects and Interfaces: Models for layer roughness (Gaussian distribution of thickness), strain profiles (exponential or linear decay from interfaces), and substrate peaks (modeled as Lorentzian functions or via coherent amplitude addition).
   • Background: Linear or polynomial background modeling to account for thermal diffuse scattering, fluorescence, and detector noise.

2. Dynamic Theory: An exact calculation based on Holý and Fewster's formulation using Parratt's algorithm for matrix propagation of complex electron density. Unlike slab models used in reflectivity codes, GenL divides unit cells into numerous slabs (typically 100) to calculate density in real space, allowing for the conceptualization of strain and roughness. This approach accounts for refraction, Darwin broadening, and total reflection, effects absent in the kinematic approximation.

Fitting Algorithm
The program utilizes a differential evolution algorithm within a genetic algorithm framework to minimize the figure of merit between experimental data and simulation. The fitting process tunes parameters such as interplanar spacing, coherent thickness, strain parameters, roughness, and background coefficients. The algorithm iterates through a population of parameter vectors until a stopping criterion (maximum iterations) is met.

Key Contributions and Results

• Dual Mode Operation: GenL offers both a Graphical User Interface (GUI) for fitting single layers using the kinematic approximation and a command-line version capable of handling arbitrary structures (including multilayers and superlattices) using either kinematic or dynamic approaches.
• Strain and Roughness Modeling: The program explicitly models out-of-plane strain profiles (exponential or linear) and layer roughness, allowing for the analysis of asymmetric Laue oscillations often associated with strain gradients or interface defects.
• Validation and Comparison:
  • GaAs Bulk Crystal: Dynamic calculations for GaAs showed excellent agreement with GenX (Parratt formalism) in the reflectivity region and correctly reproduced Darwin widths, validating the dynamic implementation.
  • Vanadium on MgO: A fit of a 105 Å V film demonstrated GenL's ability to account for tensile out-of-plane strain and roughness. The resulting coherent thickness (86.3 Å) was comparable to previous CADEM simulations (84 Å), with GenL providing a more detailed strain profile.
  • Tungsten on Al₂O₃: Fits of W (110) and W (220) peaks demonstrated the program's ability to handle overlapping substrate peaks and estimate interface roughness (4 Å), consistent with GenX reflectivity fits.
  • Iron on MgAl₂O₄: Analysis of Fe films revealed that a simple strain profile could not fully capture the asymmetry of oscillations. A bilayer model (bct Fe/bcc Fe) successfully reproduced the asymmetry, highlighting the program's utility in investigating growth modes at the thin-film limit.
  • Superlattices: While fitting superlattice data is not yet implemented in the GUI, the command-line version successfully simulated diffraction patterns for Fe/V superlattices, demonstrating extensibility to complex multilayer stacks.

Significance
The authors present GenL as a versatile, open-source tool (released under the GNU General Public License) that fills a gap in the analysis of epitaxial thin films. By combining the ability to simulate atomic ordering with robust fitting capabilities for Laue oscillations, GenL complements existing reflectivity tools like GenX. The program enables a more detailed structural picture of highly crystalline films by extracting parameters such as atomic interplanar spacings, coherent scattering plane counts, strain profiles, and crystal roughness. The authors suggest that the ability to include and fit interface and structural parameters (such as atomic steps and superlattices) will serve as a useful asset for the research community, encouraging detailed and quantitative analysis of diffraction patterns.