  The purpose of this article is to compute the radial distribution function of
liquid CO2 at 250 K and 50 bar and its expansion coefficients by spherical
harmonics as a function of intermolecular distance up to 20 angstroms for all
possible combinations of angular indices l, l prime, and m less than or equal
to 8, in contrast to past works of limited combinations. The positions and
orientations of CO2 molecules were determined by pressure-constant Monte Carlo
simulation using a Kihara potential model. It was found that the radial
distribution function had an intense peak at 4.05 angstroms and a broad band at
around 7.85 angstroms. The expansion coefficients had prominent structures
below approximately 7.5 angstroms, showing positive and/or negative peaks,
which indicated that the orientational correlations of the central molecule
were intense with the molecules present in that region. It was also found that
the peaks of some expansion coefficients were not small even when l = 8.
Furthermore, molecular orientations had a tendency to orient along one
direction presumably because of the linear form of CO2 molecule.
