We used simulations to compute the mass densities and dielectric constants of 17 polar and nonpolar organic liquids using the AM1-BCC-dPol electrostatic model and compared the results to those obtained from matched simulations using the fixed-charge AM1-BCC electrostatic model. Condensed-phase simulations of these liquids were performed with OpenMM . The MPID plugin was used to enable calculations with AM1-BCC-dPol and typed polarizabilities, using the keyword "direct" to enable direct polarization. For all liquid systems, the valence terms and van der Waals terms were assigned from the OpenFF Sage force field . Cubic liquid boxes composed of 256 molecules were built with Packmol package as simulation systems. Scripts to set up simulations are available online. Each system was equilibrated with 5 ns simulations in the NVT ensemble, followed by 15 ns NPT production using a time step of 2 fs. MD trajectories for post-processing were saved every 1 ps. For a given MD trajectory, the mass density was evaluated as: