  We investigate quantum entanglement between two (spin-1/2) fermions inside a
cylindrical harmonic trap, making use of the von Neumann entropy for the
reduced single particle density matrix as the pure state entanglement measure.
We explore the dependence of pair entanglement on the geometry and strength of
the trap and on the strength of the pairing interaction over the complete range
of the effective BCS to BEC crossover. Our result elucidates an interesting
connection between our model system of two fermions and that of two interacting
bosons.
