Numeric and symbolic evaluation of the pfaffian of general skew-symmetric matrices

Abstract

Evaluation of pfaffians arises in a number of physics applications, and for some of them a direct method is preferable to using the determinantal formula. We discuss two methods for the numerical evaluation of pfaffians. The first is tridiagonalization based on Householder transformations. The main advantage of this method is its numerical stability that makes unnecessary the implementation of a pivoting strategy. The second method considered is based on Aitken's block diagonalization formula. It yields to a kind of LU (similar to Cholesky's factorization) decomposition (under congruence) of arbitrary skew-symmetric matrices that is well suited both for the numeric and symbolic evaluations of the pfaffian. Fortran subroutines (FORTRAN 77 and 90) implementing both methods are given. We also provide simple implementations in Python and Mathematica for purpose of testing, or for exploratory studies of methods that make use of pfaffians.

13 pages, Download links: http://gamma.ft.uam.es/robledo/Downloads.html and http://www.phys.washington.edu/users/bertsch/computer.html

References in corpus (2)

• The sign of the overlap of HFB wave functions
• Exact Algorithm for Sampling the 2D Ising Spin Glass

Cited by in corpus (18)

• Mean field and beyond description of nuclear structure with the Gogny force: A review
• Efficient numerical computation of the Pfaffian for dense and banded skew-symmetric matrices
• Symmetry restoration in mean-field approaches
• A microscopic benchmark-study of triaxiality in low-lying states of 76Kr
• Solution of the Skyrme-Hartree-Fock-Bogolyubov equations in the Cartesian deformed harmonic-oscillator basis. (VIII) hfodd (v2.73y): a new version of the program
• Evaluation of overlaps between arbitrary Fermionic quasiparticle vacua
• Effect of Strong Disorder in a 3-Dimensional Topological Insulator: Phase Diagram and Maps of the Z2 Invariant
• Toward extremes of angular momentum: Application of the Pfaffian algorithm in realistic calculations
• Many-body computations by stochastic sampling in Hartree-Fock-Bogoliubov space
• Efficient method to perform quantum number projection and configuration mixing for most general mean-field states
• The Faddeev-LeVerrier algorithm and the Pfaffian
• A Pseudo-BCS Wavefunction from Density Matrix Decomposition:Application in Auxiliary-Field Quantum Monte Carlo
• A Pfaffian formulation for matrix elements of three-body operators in multiple quasi-particle configurations
• A factorization algorithm to compute Pfaffians
• Lateral critical Casimir force in 2D Ising strip with inhomogeneous walls
• Counting statistics in finite Fermi systems: illustrations with the atomic nucleus
• Universality in the tripartite information after global quenches: spin flip and semilocal charges
• Fusion mechanism for quasiparticles and topological quantum order in the lowest Landau level