Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aeif_v2_0.tar.gz(3233 Kbytes)
Manuscript Title: Stochastic Hyperfine Interactions Modeling Library - Version 2
Authors: Matthew O. Zacate, William E. Evenson
Program title: SHIML
Catalogue identifier: AEIF_v2_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 199(2016)180
Programming language: C.
Computer: Any.
Operating system: LINUX, OS X.
RAM: Variable
Keywords: Perturbed angular correlation spectroscopy, Mössbauer effect, TDPAC, Stochastic models, Hyperfine methods.
PACS: 76.80.+y, 76.60.-k, 76.75.+i.
Classification: 7.4.

External routines: TAPP [1], BLAS [2], a C-interface to BLAS [3], and LAPACK [4]. Additionally, GSL [3] is needed to compile the example code that simulates Mössbauer spectra.

Does the new version supersede the previous version?: No

Nature of problem:
In condensed matter systems, hyperfine methods such as nuclear magnetic resonance (NMR), Mössbauer effect (ME), muon spin rotation (μSR), and perturbed angular correlation spectroscopy (PAC) measure electromagnetic fields due to electronic and magnetic structure within Angstroms of nuclear probes through the hyperfine interaction. When interactions fluctuate at rates comparable to the time scale of a hyperfine method, there is a loss in signal coherence, and spectra in the time domain are damped while spectra in the frequency domain are broadened. The degree of damping or broadening can be used to determine fluctuation rates, provided that theoretical expressions for spectra can be derived for relevant physical models of the fluctuations. SHIML provides routines to help researchers quickly develop code to incorporate stochastic models of fluctuating hyperfine interactions in calculations of hyperfine spectra.

Solution method:
Calculations are based on the method for modeling stochastic hyperfine interactions for PAC by Winkler and Gerdau [5]. The method is extended to include other hyperfine methods following the work of Dattagupta [6]. The code provides routines for reading model information from text files, allowing researchers to develop new models quickly without the need to modify computer code for each new model to be considered.

Reasons for new version:
The original version of the library provided support only for those methods that measure hyperfine interactions on one spin state of the nuclear probe. As such, it excluded important hyperfine methods that measure the interactions on two spin states such as Mössbauer spectroscopy and nuclear resonant scattering of synchrotron radiation. The present version of SHIML provides the necessary support for such double spin-state methods while maintaining backward compatibility for code already developed using the original version.

Summary of revisions:
Routines now check that values representing nuclear spins are positive integers or positive half-integers. Additional utility functions are provided to make it easier for code developers to calculate Hamiltonians of electric quadrupole interactions. A correction was made to the portion of code responsible for calculating the Blume matrix of single spin-state methods; however, this change will not alter results obtained from single spin-state simulations using version 1 of the library. The remaining revisions support calculations for double spin-state methods. (1) Model-file syntax is expanded in order to allow users to specify different hyperfine interactions for ground and excited spin states and to input isomer shifts. (2) New routines for initialization and for Blume-matrix calculations are included for the double spin-state case. (3) New example code is provided to illustrate how SHIML can be used to simulate Mössbauer spectra of polycrystalline samples for pure dipole or pure quadrupole transitions; background information about the Mössbauer examples can be found in ref. 7. Finally, updated software documentation is included in a User's Guide as a PDF file in the code distribution.

Running time:

[1] M. O. Zacate, The Adjustable Parameter Package, Technical Report 2, CINSAM Grant 2006-R7 (unpublished); available for download at http://tapp.nku.edu/.
[2] L. S. Blackford et al., ACM Trans. Math. Soft. 28 (2002) 135; J. Dongarra, International Journal of High Performance Applications and Supercomputing 16 (2002) 1; http://www.netlib.org/blas/.
[3] M. Galassi et al., GNU Scientific Library Reference Manual, third edition (2009); available for download at http://www.gnu.org/software/gsl/.
[4] E. Anderson et al., LAPACK Users' Guide, third ed. (Society for Industrial and Applied Mathematics, Philadelphia, PA, 1999); http://www.netlib.org/lapack/.
[5] H. Winkler, E. Gerdau, Z. Phys. 262 (1973) 363.
[6] S. Dattagupta, Hyperfine Interact. 11 (1981) 77.
[7] M. O. Zacate, W. E. Evenson, Hyperfine Interact. 231 (2015) 143.