A methodology for the computational prediction of the crystal structures and resolution efficiency for diastereomeric salt pairs is developed by considering the polymorphic system of the diastereomeric salt pair (R)-1-phenylethylammonium (R/S)-2-phenylpropanoate. To alleviate the mathematical complexity of the search for minima in the lattice energy due to the presence of two flexible entities in the asymmetric unit, the range of rigid-body lattice energy global optimizations was guided by a statistical analysis of the Cambridge Structural Database for common ion-pair geometries and ion conformations. A distributed multipole model for the dominant electrostatic interactions and high-level ab initio calculations for the intramolecular energy penalty for conformational distortions are used to quantify the relative stabilities of the p- and n-salt forms. While the ab initio prediction of the known structure of the p-salt as the most stable structure was insensitive to minor changes in the rigid-ion conformations considered, the relative stabilities of the known polymorphs and hypothetical structures of the n-salt were very sensitive. Although this paper provides a significant advance over traditional search algorithms and empirical force fields in determining the structures and relative stabilities of diastereomeric salt pairs, the sensitivity of the computed lattice energies to the fine details of the ion conformations overtaxes current computational models and renders the design of diastereomeric resolution processes by computational chemistry a challenging problem.