A vectorized implementation of the Debye Equation on CPU and GPU
Project description
ChemRxiv | [Paper]
DebyeCalculator
Welcome to DebyeCalculator
! This is a simple tool for calculating the scattering intensity $I(Q)$ through the Debye scattering equation, the Total Scattering Structure Function $S(Q)$, the Reduced Total Scattering Function $F(Q)$, and the Reduced Atomic Pair Distribution Function $G(r)$ from an atomic structure.
The Debye scattering equation can be used to compute the scattering pattern of any atomic structure and is commonly used to study both crystalline and non-crystalline materials with a range of scattering techniques like powder diffraction, total scattering with pair distribution function and small-angle scattering. Although the Debye scattering equation is extremely versatile, the computation of the double sum, which scales O(N2), has limited the practical use of the equation.
Here, we provide an optimised code for the calculation of the Debye scattering equation on Graphics processing units (GPUs) which accelerate the calculations with orders of magnitudes.
Installation
Install with pip
Run the following command to install the DebyeCalculator package. (Requires: Python >=3.9, <3.12)
pip install DebyeCalculator
Install locally
Run the following command to install the DebyeCalculator package.
pip install .
or
python setup.py install
Usage
Interactive mode at Google Colab
Follow the instructions in our Interactive Google Colab notebook and try to play around.
Example usage
from debyecalculator import DebyeCalculator
# Initialise calculator object
calc = DebyeCalculator(qmin=1.0, qmax=8.0, qstep=0.01)
# Define structure path
XYZ_path = "some_path/some_file.xyz"
# Print object to view all parameters
print(calc)
## [OUTPUT] DebyeCalculator{'qmin': 1.0, 'qmax': 8.0, 'qstep': 0.01, 'rmin': 0.0, 'rmax': 20.0, ...}
# Calculate Powder (X-ray) Diffraction
Q, I = calc.iq(structure_path=XYZ_path)
# Update parameters for Small Angle (Neutron) Scattering
calc.update_parameters(qmin=0.0, qmax=3.0, qstep=0.01, radiation_type="neutron")
# Calculate Small Angle (Neutron) Scattering
calc.iq(structure_path=XYZ_path)
# Update parameters for Total Scattering with Pair Distribution Function analysis
calc.update_parameters(qmin=1.0, qmax=30.0, qstep=0.1, radiation_type="xray")
# Calculate Pair Distribution Function
r, G = calc.gr(structure_path=XYZ_path)
.....
Authors
Frederik L. Johansen1, 2
Andy S. Anker1
1 Department of Chemistry & Nano-Science Center, University of Copenhagen, Denmark
2 Department of Computer Science, University of Copenhagen, Denmark
Should there be any questions, desired improvements or bugs please contact us on GitHub or through email: frjo@di.ku.dk and andy@chem.ku.dk.
Cite
If you use our code or our results, please consider citing our paper. Thanks in advance!
@article{Johansen_anker2023debye,
title={A GPU-Accelerated Open-Source Python Package for Calculating Powder Diffraction, Small-Angle-, and Total Scattering with the Debye Scattering Equation},
author={Frederik L. Johansen, Andy S. Anker, Ulrik Friis-Jensen, Erik B. Dam, Kirsten M. Ø. Jensen, Raghavendra Selvan},
journal={ChemRxiv}
year={2023}}
Contributing to the software
We welcome contributions to our software! To contribute, please follow these steps:
- Fork the repository.
- Make your changes in a new branch.
- Submit a pull request.
We'll review your changes and merge them if they meet our quality standards, including passing all unit tests. To ensure that your changes pass the unit tests, please run the tests locally before submitting your pull request. You can also view the test results on our GitHub repository using GitHub Actions.
Reporting issues
If you encounter any issues or problems with our software, please report them by opening an issue on our GitHub repository. Please include as much detail as possible, including steps to reproduce the issue and any error messages you received.
Seeking support
If you need help using our software, please reach out to us on our GitHub repository. We'll do our best to assist you and answer any questions you have.
Project details
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
Hashes for debyecalculator-1.0.1-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 9d576078665c67bb2a1f1398ac4d526b06c058dba8c2af712946c843d021ace5 |
|
MD5 | f7785c7f145937fe02590b49b367c2d2 |
|
BLAKE2b-256 | 338e541069c1aa08d757c8dd9cfcf7a095fe61916b60c59fb17bc603d12db47f |