A new formulation to calculate general HFB matrix elements through Pfaffian

A new formula is presented for the calculation of matrix elements between multi-quasiparticle Hartree-Fock-Bogoliubov (HFB) states. The formula is expressed in terms of the Pfaffian, and is derived by using the Fermion coherent states with Grassmann numbers. It turns out that the formula corresponds to an extension of generalized Wick's theorem and simplifies the combinatorial complexity resulting from practical applications of generalized Wick's theorem by unifying the transition density and the transition pairing tensor in the HFB theory. The resultant formula is simpler and more compact than the traditional description of matrix elements of general many-body operators. In addition, through the derivation of our new formula, we found that the Pfaffian version of the Lewis Carroll formula corresponds to the relation conjectrured by Balian and Brezin for the HFB theory in 1969.