Simulations have been performed at 480 K for pure melts of atactic, isotactic and syndiotactic polypropylene chains with a degree of polymerization of 50 and composition C150H302, and also for the three 50:50 mixtures of two of these species. The simulations are performed on a high coordination lattice, with incorporation of short range intramolecular interactions from a rotational isomeric state model of polypropylene, and incorporation of long range interactions defined by a Lennard-Jones potential energy function for the interaction of pairs of molecules of propane. Both the rotational isomeric state model and Lennard-Jones parameters were taken directly from the literature (Suter and Prausnitz, respectively). The efficiency of the simulation on the sparsely occupied high coordination lattice facilitates the equilibration of the one- and two-component melts within accessible computer time. Onset of a tendency for demixing of isotactic and syndiotactic polypropylene is apparent in the intermolecular pair correlation functions. No such demixing occurs with isotactic and atactic polypropylene. Both of these predictions from the simulation are consistent with experimental results in the literature (Maier and Lohse, respectively). The simulation produces an ambigious prediction for the melt of atactic and syndiotactic polypropylene. This melt has been reported (Maier ) to exhibit phase separation, but less strongly than the isotactic-syndiotactic system, for which the simulation makes an unambiguous (and correct) prediction. The physical origin of the tendency for demixing in the simulations is identified as the differences in the preferred local conformations of polypropylene chains with various stereochemical sequences. This driving force is an example of "conformational asymmetry" induced solely by differences in stereochemical sequence.