Mass flows and angular momentum density for p(x)+ip(y) paired fermions in a harmonic trap

We present a simple two-dimensional model of a p(x)+ip(y) superfluid in which the mass flow that gives rise to the intrinsic angular momentum is easily calculated by numerical diagonalization of the Bogoliubov-dc Gennes operator. We find that, at zero temperature and for constant director 1, the mass flow closely follows the Ishikawa-Mermin-Muzikar formula j(mass)=1/2 curl/(rho hl/2). (C) 2007 Elsevier Inc. All rights reserved.