In this paper, we
present a new molecular model that can accurately
predict thermodynamic liquid state and phase-change properties for
organosilicon molecules including several functional groups (alkylsilane,
alkoxysilane, siloxane, and silanol). These molecules are of great
importance in geological processes, biological systems, and material
science, yet no force field currently exists that is widely applicable
to organosilicates. The model is parametrized according to the recent
Polarization-Consistent Approach (PolCA), which allows for polarization
effects to be incorporated into a nonpolarizable model through
post facto
correction terms and is therefore consistent
with previous parametrizations of the PolCA force field. Alkyl groups
are described by the United-Atom approach, bond and angle parameters
were taken from previous literature studies, dihedral parameters were
fitted to new quantum chemical energy profiles, point charges were
calculated from quantum chemical optimizations in a continuum solvent,
and Lennard-Jones dispersion/repulsion parameters were fitted to match
the density and enthalpy of vaporization of a small number of selected
compounds. Extensive validation efforts were carried out, after careful
collection and curation of experimental data for organosilicates.
Overall, the model performed quite well for the density, enthalpy
of vaporization, dielectric constant, and self-diffusion coefficient,
but it slightly overestimated the magnitude of self-solvation free
energies. The modular and transferable nature of the PolCA force field
allows for further extensions to other types of silicon-containing
compounds.