A coarse-grained model of atactic polystyrene, in which meso and racemo diads are represented as single '' superatoms,'' parametrized using Iterative Boltzmann Inversion, has been subjected to connectivity-altering Monte Carlo simulations in order to simulate monodisperse atactic polystyrene melts of molar mass up to 210000 g mol(-1) at 500 or 413 K and 1 bar. Analysis of the Monte Carlo results reveals excellent equilibration of chain conformations at all length scales. Chain dimensions, as determined from the mean square end-to-end distance, the mean square radius of gyration, and simulated Kratky plots of the single-chain scattering function, are in excellent agreement with experiment. The equilibrated long-chain configurations are reduced to entanglement networks via topological analysis with the CReTA algorithm. The resulting Kuhn length of primitive paths provides an excellent estimate of the molar mass between entanglements and of the entanglement tube diameter extracted from plateau modulus measurements. The distribution of strand lengths between entanglements, when appropriately reduced, follows the same master curve as previously determined distributions of polyethylene, cis-1,4 polybutadiene, and poly(ethylene terephthalate). A new strategy is introduced for reverse mapping the long-chain coarse-grained configurations into detailed united-atom configurations in a manner that preserves the sequence of diad types along the chains. This strategy employs local '' flip '' Monte Carlo moves to relax the reverse-mapped configurations. Relaxation starts using bonded interactions only, and proceeds by gradually introducing nonbonded interactions. Final relaxation is achieved via short-time canonical molecular dynamics simulation. Predicted wide-angle X-ray diffraction patterns from reverse-mapped configurations are indistinguishable from those of short-chain melts equilibrated directly in the united atom representation using molecular dynamics, and in favorable agreement with experiment. Distributions of torsion angles and pairs of successive torsion angles in the reverse-mapped configurations exhibit some deviations from the corresponding distributions of directly equilibrated short-chain united atom melts and from experimental NMR measurements.