We present an approach for parametrizing spectroscopic maps of carbonyl groups against experimental IR absorption spectra. The model correlates electric fields sampled from molecular dynamics simulations with vibrational frequencies and line shapes in different solvents. We perform an exhaustive search of parameter combinations and optimize the parameter values for the ester carbonyl stretching mode in ethyl acetate by comparing to experimental FTIR spectra of the small molecule in eight different solvents of varying polarities. Hydrogen-bonding solvents require that the peaks are fit independently for each hydrogen bond ensemble to compensate for improper sampling in molecular dynamics simulations. Spectra simulated using the optimized electrostatic map reproduce C═O IR absorption spectra of ethyl acetate with a line center RMSD error of 4.9 cm(-1) over 12 different solvents whose measured line centers span a 45 cm(-1) range. In combination with molecular dynamics simulations, this spectroscopic map will be useful in interpreting spectra of ester groups in heterogeneous environments such as lipid membranes.