In this theoretical work we evaluate how the chemical environment influences some features presented in the infrared spectrum, such as band intensities and band location of embedded species in icy matrices. The calculations were performed employing the Polarized Continuum Model (PCM) approach with the second-order Møller–Plesset perturbation theory (MP2) level using the Gaussian 09 package. Here, we simulate the effects of molecular vicinity around embedded species in terms of the effects of the dielectric constant (ε) of the icy and solid samples. Gas phase calculation was also performed for comparison purpose. The investigated embedded single molecules were CO, CO₂, CH₄, NH₃, SO₂ HCOOH, CH₃OH and also H₂O. The results suggest that for most vibrational modes, the strengths of IR bands show an increase with ε, which implies they also decrease with respect to porosity. The frequency shifts showed opposite behavior in relation to the band strengths, with few exceptions. A correlation between calculated band intensities with the band strengths A (taken from literature) was determined and described by a linear function I ∼ 6 × 10¹⁸ A [km mol⁻¹], with A in unity of cm per molecule. In addition, an associative exponential function was adjusted to the studied dataset to characterize the evolution of frequency-shift and intensity-shift and band strength ratio as function of the dielectric constant. Since astrophysical ice mantles over cold dust grains can vastly vary in composition in space (having different dielectric constants) they are a challenge to be well characterized. Therefore, this work can help the astrochemistry community to better understand astrophysical ices and its observations in the infrared.