The singlet and pair distribution functions of a Fermi gas are evaluated for low temperatures. It is assumed that the particles interact with a hard-sphere potential and have an arbitrary spin, and the evaluation is made to first order in the hard-sphere diameter. The binary-kernel method developed by Siegert, Lee, and Yang is used and the cluster-expansion theory given recently by Isihara and Yee for the distribution functions is applied. The pair distribution function is determined by a dimensionless variable $\frac{r}{\ensuremath{\lambda}}$, with $\ensuremath{\lambda}$ being the thermal de Broglie wavelength. Thus, changing the distance corresponds to varying the temperature. Moreover, the pair distribution function depends on the density. Therefore, its asymptotic expression for $T\ensuremath{\rightarrow}0$ depends on the ways to approach the limit. If we keep the density finite and reduce the temperature, the gas will be degenerate when $\ensuremath{\lambda}$ is of order of ${n}^{\ensuremath{-}\frac{1}{3}}$, where n is the number density. The pair distribution function for this case is an oscillating function of the Fermi momentum. On the other hand, if we reduce the density first and then decrease the temperature, the pair distribution function decays as ${r}^{\ensuremath{-}4}$, and the free-particle diagrams give a contribution of the same decay.
The pair distribution function for a hard sphere Bose system has been calculated by using a method in which the hard sphere potential is replaced by the so-called pseudo-potential. The problem is carried through using the quasi-particle formalism of Bogoliubov. In this calculation one considers a system of bosons interacting via the pseudo-potential which can be regarded as a weak interaction when the system is very dilute. The Hamiltonian of the system is diagonalized by means of a canonical transformation which has the effect of separating the energy into two parts, one of which is the ground state energy and the other corresponds to an ideal Bose gas composed of ``elementary excitations'' or ``quasi-particles.'' This formalism is applied to the calculation of the pair distribution function. This quantity is calculated by averaging over a grand canonical ensemble constructed with the total number of elementary excitations. A result which is seen to be valid both for the condensed and gaseous phases of the system and also for any distance r between the particles is obtained. A discussion of the difference between these results and the ones obtained in a previous calculation and a comparison between both results and the experimental ones is also given.