×

Parallel implementation of the MAGPACK package for the analysis of high-nuclearity spin clusters. (English) Zbl 1219.81273

Summary: Molecular clusters are formed by a finite number of exchange-coupled paramagnetic centers and they are model systems between molecules and extended solids. In order to simulate their properties and extrapolate to solids, the size of the systems to be treated should be as large as possible. In this context, the use of efficient parallel codes is essential. We present the parallel programs ParAni and ParIso, for anisotropic and isotropic models, that enable the calculation of large energy matrices in parallel and the subsequent computation of the relevant spectral information. The evaluation of the matrix elements is based on the serial package MAGPACK that uses the irreducible tensor operators technique and takes into account all kinds of anisotropic and isotropic magnetic interactions. To obtain the eigenvalues, the energy matrix is partially diagonalized by means of the SLEPc library. The calculation of eigenvalues and eigenvectors of these spin clusters enables us to evaluate the bulk magnetic properties (magnetic susceptibility and magnetization) as well as the spectroscopic properties (inelastic neutron scattering spectra). The results are encouraging in terms of parallel efficiency and open the way to address very challenging problems.

MSC:

81V55 Molecular physics
81Q05 Closed and approximate solutions to the Schrödinger, Dirac, Klein-Gordon and other equations of quantum mechanics
82D40 Statistical mechanics of magnetic materials
81Q10 Selfadjoint operator theory in quantum theory, including spectral analysis
68T10 Pattern recognition, speech recognition
68Q10 Modes of computation (nondeterministic, parallel, interactive, probabilistic, etc.)
81T80 Simulation and numerical modelling (quantum field theory) (MSC2010)
81-08 Computational methods for problems pertaining to quantum theory
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Coronado, E.; Dunbar, K. R., Preface for the forum on molecular magnetism: The role of inorganic chemistry, Inorg. Chem., 48, 3293-3295 (2009)
[2] Borras-Almenar, J. J.; Clemente-Juan, J. M.; Coronado, E.; Tsukerblat, B. S., MAGPACK: A package to calculate the energy levels, bulk magnetic properties, and inelastic neutron scattering spectra of high nuclearity spin clusters, J. Comput. Chem., 22, 9, 985-991 (2001)
[3] Hernandez, V.; Roman, J. E.; Vidal, V., SLEPc: A scalable and flexible toolkit for the solution of eigenvalue problems, ACM Trans. Math. Software, 31, 3, 351-362 (2005) · Zbl 1136.65315
[4] Borras-Almenar, J. J.; Clemente-Juan, J. M.; Coronado, E.; Tsukerblat, B. S., High-nuclearity magnetic clusters: Generalized spin Hamiltonian and its use for the calculation of the energy levels, bulk magnetic properties, and inelastic neutron scattering spectra, Inorg. Chem., 38, 6081-6088 (1999)
[5] Parlett, B. N., The Symmetric Eigenvalue Problem (1980), Prentice-Hall: Prentice-Hall Englewood Cliffs, NJ, reissued with revisions by SIAM, Philadelphia, 1998 · Zbl 0431.65016
[6] (Bai, Z.; Demmel, J.; Dongarra, J.; Ruhe, A.; van der Vorst, H., Templates for the Solution of Algebraic Eigenvalue Problems: A Practical Guide (2000), Society for Industrial and Applied Mathematics: Society for Industrial and Applied Mathematics Philadelphia, PA) · Zbl 0965.65058
[7] Wu, K.; Simon, H., Thick-restart Lanczos method for large symmetric eigenvalue problems, SIAM J. Matrix Anal. Appl., 22, 2, 602-616 (2000) · Zbl 0969.65030
[8] Gilbert, D.; Roman, J. E.; Garland, W. J.; Poehlman, W. F.S., Simulating control rod and fuel assembly motion using moving meshes, Ann. Nucl. Energy, 35, 2, 291-303 (2008)
[9] Roman, J. E.; Kammerer, M.; Merz, F.; Jenko, F., Fast eigenvalue calculations in a massively parallel plasma turbulence code, Parallel Comput., 36, 5-6, 339-358 (2010) · Zbl 1204.68267
[10] Stewart, G. W., A Krylov-Schur algorithm for large eigenproblems, SIAM J. Matrix Anal. Appl., 23, 3, 601-614 (2001) · Zbl 1003.65045
[11] Lehoucq, R. B.; Sorensen, D. C.; Yang, C., ARPACK Users’ Guide, Solution of Large-Scale Eigenvalue Problems by Implicitly Restarted Arnoldi Methods (1998), Society for Industrial and Applied Mathematics: Society for Industrial and Applied Mathematics Philadelphia, PA · Zbl 0901.65021
[12] S. Balay, K. Buschelman, V. Eijkhout, W. Gropp, D. Kaushik, M. Knepley, L.C. McInnes, B. Smith, H. Zhang, PETSc users manual, Tech. Rep. ANL-95/11, Revision 3.0.0, Argonne National Laboratory, December 2008.; S. Balay, K. Buschelman, V. Eijkhout, W. Gropp, D. Kaushik, M. Knepley, L.C. McInnes, B. Smith, H. Zhang, PETSc users manual, Tech. Rep. ANL-95/11, Revision 3.0.0, Argonne National Laboratory, December 2008.
[13] Hernandez, V.; Roman, J. E.; Tomas, A., Parallel Arnoldi eigensolvers with enhanced scalability via global communications rearrangement, Parallel Comput., 33, 7-8, 521-540 (2007)
[14] Karypis, G.; Kumar, V., Parallel multilevel \(k\)-way partitioning scheme for irregular graphs, SIAM Rev., 41, 2, 278-300 (1999) · Zbl 0918.68073
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.