Multi-criteria optimization is used for developing molecular models for ethylene carbonate (EC) and propylene carbonate (PC), organic solvents commonly used in Li-ion batteries. The molecular geometry and partial charges of the solvents are obtained from quantum mechanical calculations. Using a novel optimization strategy that combines systematic variations of the Lennard-Jones parameters with a reduced units approach, the models are fitted to experimental data on the liquid density, vapor pressure, relative permittivity, and self-diffusion coefficient. Since no experimental data for the self-diffusion coefficient of pure EC were available in the literature, they are measured in this work using a gradient-based nuclear magnetic resonance technique. For all pure component properties, excellent agreement between experiment and simulation is obtained. Moreover, the predictive capabilities of the new solvent models are assessed by comparison to experimental data for the liquid density and relative permittivity of mixtures of EC and PC. In addition, molecular models for the anions PF6-, BF4-, and ClO4- in solutions of their lithium electrolytes in PC are developed using experimental data on the solution densities. Finally, the self-diffusion coefficients of LiPF6 in PC and in aqueous solution are predicted and compared, showing that diffusion is much slower in the organic solution due to the formation of larger solvent shells around the ions. Furthermore, an analysis of the radial distribution functions in these solutions suggests that the ions have much less impact on the structure of the solvent PC than on water.