I am trying to simulate IR and Raman of a catalyst and trying to see the effect of Na+ on the frequency. For the catalyst I am getting a very good match for Raman (after converting to intensity), but for IR some strong experimentally observed modes (like S=O str.) are appearing as medium or weak band in theory. I am taking a dimer (obtained from crystal structure and then optimized) in gas phase to simulate the solid. I know IR intensity strongly depends upon intermolecular interaction which is not described properly by taking dimer. How to match the intensity better?
Then I tried to add Na+ cation to the catalyst to see the shift in frequency. Shift in frequency is qualitatively changing as observed in experiment. Now I want to add a counter ion (BPh4), to get the contribution from it as seen in experiment. Probably since I am doing the calculation in gas phase adding BPh4- close to Na+ causing the Na to go little distant from catalyst and shift in catalyst frequency is diminished. The crystal structure of the complex is not known, so I do not have any idea of where to add it. I do see interaction of BPh4 with the catalyst, 4 Phenyl ring breathing mode is split in the complex. Could you suggest how to solve the problem?
Thanks and regards.
Partha