SpinW / spinw

SpinW Matlab library for spin wave calculation
http://www.spinw.org
GNU General Public License v3.0
35 stars 15 forks source link

User request: preferred orientation in powder simulations #195

Open mducle opened 2 months ago

mducle commented 2 months ago

A user has requested the ability to perform a powder simulation for a sample with preferred orientation of crystallites, using the formulism used in powder crystallography, specifically the March formula (article in German - see also this article) as implemented in FullProf which posits a distribution of preferred Q vectors defined by an angle alpha between Q and a preferred orientation direction n (called PREF by FullProf) as:

W(\alpha) = f + (1 - f)*\left( r^2\cos^2(\alpha) + \frac{1}{r}\sin^2(\alpha) \right)^{-3/2}

Where the March parameter r (called G1 by FullProf) in the range [-1,1] determines the habit of the preferred orientation, and the parameter f (called G2 by FullProf) is the fraction of sample which has normal (random) orientations.

It would be nice if in powspec we can implement a preferred orientation mechanism where the user can input the preferred direction n, the r and f parameters and we can then sample on the sphere surface according to the March formula.