Calculation of Vibrational Spectra of Linear Tetrapyrroles. 3. Hydrogen-Bonded Hexamethylpyrromethene Dimers

The structure and vibrational spectra of hexamethylpyrromethene (HMPM) have been investigated by X-ray crystallography, IR and Raman spectroscopies, and density functional theory calculations. HMPM crystallizes in the form of dimers, which are held together by bifurcated N−H(···N)2 hydrogen bonds, involving one intramolecular and one intermolecular N−H···N interaction. The monomers are essentially planar, and the mean planes of the monomers lie approximately perpendicular to one another, so that the four N atoms in the dimer form a distorted tetrahedron. The structure of the HMPM dimer is well-reproduced by B3LYP/6-31G* calculations. A comparison of the calculated geometry of the dimer with that of the monomer reveals only small changes in the N−H···N entity and the methine bridge angles upon dimerization. These are a result of weakening of the intramolecular N−H···N hydrogen bond and the formation of a more linear N−H···N intermolecular hydrogen bond. Using an empirical relation between the shift of the N−H stretching frequency of pyrrole and the enthalpy of adduct formation with bases [Nozari, M. S.; Drago, R. S. J. Am. Chem. Soc. 1970, 92, 7086−7090], estimates of the strength of the intra- and intermolecular hydrogen bonds are obtained. IR and Raman spectroscopies of HMPM and its isotopomers deuterated at the pyrrolic nitrogen atom and at the methine bridge reveal that the molecule is monomeric in nonpolar organic solvents but dimeric in a solid Ar matrix and in KBr pellets. The matrix IR spectra show a splitting of vibrational modes for the dimer, particularly those involving the N−H coordinates. Due to intrinsic deficiencies of the B3LYP/6-31G* approximation, a satisfactory reproduction of these modes of the monomeric and dimeric HMPM requires specific adjustments of the NH scaling factors for the calculated force constants and, in the case of the NH out-of-plane modes of HMPM dimers, also of intra- and intermolecular coupling constants. This parametrization does not significantly affect the other calculated modes, which in general reveal a very good agreement with the experimental data.