Huigang Wang*,
Hang Xu,
Qiuna Liu and
Xuming Zheng
Department of Chemistry, Key Laboratory of Advanced Textiles Materials and Manufacture Technology of the Ministry of Education, Engineering Research Center for Eco-Dyeing and Finishing of Textiles of the Ministry of Education, Zhejiang Sci-Tech University, Hangzhou 310018, P. R. China. E-mail: zdwhg@163.com; hugwang@ucdavis.edu; Tel: +86-571-8684-3627
First published on 21st August 2020
This article aims to correlate the noncoincidence effect phenomenon with the aggregation state of acetylacetone CO stretching in a binary mixture. CO stretching noncoincidence effect (NCE) was observed not only between IR and Raman spectra but also between the isotropic and anisotropic Raman spectra of acetonylacetone. The difference in CO stretching wavenumbers of the isotropic and anisotropic Raman spectra (NCE value) in a binary mixture at different concentrations has been calculated. We found that both isotropic and anisotropic Raman wavenumbers of CO stretching increase with the dilution of acetonylacetone by CCl4 while the NCE value decreases. These noncoincidence and concentration effect phenomena seem to go against the quantum theory. Herein, we proposed an aggregation-induced split (AIS) model to explain the NCE phenomenon and concentration effect. The experimental data were consistent with the DFT calculations performed at the B3LYP-D3/6-311++G (d,p) levels based on the proposed model. The dynamics of transformation from monomers to an aggregated structure can be easily controlled by tuning the concentration. Solvent dependent experiments show that the value of NCE decreased with the increase of the solvent dielectric constant at the same concentration, which is in accordance with Logan's theory.
Our group concentrated on these phenomena for several years. First, we discovered NCEs in CS vibrational modes.21,22 The difference between the isotropic and anisotropic peak frequencies of CS stretching for ethylene trithiocarbonate was determined to be 4.60 cm−1.21 This difference decreased upon dilution. These NCEs and concentration effects made us believe that CS stretching is not a single vibrational mode, but a complex vibrational mode that is beyond spectroscopic resolution.22
Matrix isolation is a powerful technique for the enhancement of spectroscopic resolution because the translational and rotational motions can freeze below 6 K.23,24 We applied this technique to study the CO stretching NCE behavior of acetone, the most investigated model molecule for NCE phenomena.25,26 Acetone was isolated in an argon matrix and the Raman spectra were collected at different annealing temperatures. Single, double, and triple peaks were detected separately at different temperatures for the CO vibration, and the isotropic and anisotropic spectra for each wavenumber overlapped fairly well with no NCE.26 Thus, an aggregation-induced split (AIS) model has been proposed to explain the acetone CO vibration NCE phenomenon and its concentration effect.25,26 The polar bond vibration coupling tends to align the molecules to reduce the potential energy and increase the attraction. The alignment and reorganization stabilized the aggregation structure to form dimers, trimers or clusters. Polar bonds such as CO, CS, SO, C–N, and C–O especially play an important role in these interactions. The vibration coupling between adjacent polar bonds split the degenerate vibrations to two vibrational modes: in phase vibration and out of phase vibration. The pairs with prominent vibrational wavenumber difference and depolarization ratio difference produced “NCE effects” with the limited spectroscopic resolution; however, in high-resolution spectroscopy, they were well separated.25,26 Acetylacetone is the simplest β-diketone that bears two tautomeric forms (keto and enol). Due to the existence of intramolecular hydrogen bonds,27,28 the enol form is more stable in the gas phase or in weak polar solvents, while the keto form is more stable in polar solvents, especially protic solvents. The keto structure will form intermolecular hydrogen bonds with protic solvents. Acetylacetone is a good molecule for studying solvent effects. This article aims to investigate the solvent polarity influence on its aggregation state.
The AIS model can explain most of the NCE phenomena.25,26,29–31 In this work, 2,5-hexanedione was chosen to investigate the CO coupling behavior and the NCE phenomenon. 2,5-Hexanedione contains two carbonyl groups, and there are two possible CO interactions,.intramolecular and intermolecular. Intramolecular CO interactions have no relation to their neighboring molecules, and do not show concentration effects, while intermolecular CO interactions show both NCE phenomena and concentration effects. In this work, the isotropic and anisotropic Raman bands of 2,5-hexanedione has been recorded and their concentration dependent behavior has been investigated to know the exact molecular interaction structure.
The Fourier transform (FT)-Raman and FT-IR spectra were obtained with 2 cm−1 resolution using a FT-Raman spectrometer at 1064 nm excitation (Thermo Nicolet 960, Thermo Fisher Nicolet, USA) and a FT-IR spectrometer (Thermo Nicolet avatar 370, Thermo Fisher Nicolet, USA).
We converted Raman intensity from scattering activities using the Multiwfn software (http://multiwfn.codeplex.com/releases). Multiwfn is an extremely powerful wavefunction analysis program and supports almost all of the most important wavefunction analysis methods.33 In our case, the excitation wavenumber is 20491.80 cm−1 (corresponding to 488 nm) and the conversion relationship is shown as34
With the parallel and perpendicular polarized Raman spectra we get the isotropic and anisotropic components of Raman spectra using the equation:27
Ianiso(vVH) = IVH(vVH) |
With the DFT calculated Gaussian Output File of acetonylacetone, we also get the isotropic and anisotropic spectrum theoretically from the dimer model DFT calculations. Gaussview or Multiwfn can transform the dimer frequencies into Raman spectra. The original output file can be regarded as IVV(vVV). Since we were interested in the CO vibration, we focused our attention on the CO frequency and their depolarization ratio. With the depolarization ratio band corresponding to the Raman activity of CO vibration, we can get the depolarized Raman activity of CO vibration. With Gaussview, we can transform the frequencies to the depolarized Raman spectra, IVH of the CO vibration. Finally, with we can get the isotropic Raman activity, Iisospectra (see ESI† for more demonstrations).
Modes | Computed/cm−1 | Experiment/cm−1 | Descriptions | ||||
---|---|---|---|---|---|---|---|
Monomer freq. | Dimer | Raman | IR | ||||
cis | anti | Freq. | D radio | ||||
ν1 | 3144 | 3138 | 3141/3138 | 0.71/0.72 | C–H stretch | ||
ν2 | 3139 | 3138 | 3137/3134 | 0.68/0.71 | C–H stretch | ||
ν3 | 3107 | 3085 | 3106/3104 | 0.32/0.72 | C–H stretch | ||
ν4 | 3094 | 3085 | 3097/3095 | 0.72/0.55 | C–H stretch | ||
ν5 | 3087 | 3064 | 3090/3088 | 0.75/0.65 | C–H stretch | ||
ν6 | 3053 | 3064 | 3083/3075 | 0.46/0.74 | C–H stretch | ||
ν7 | 3033 | 3034 | 3044/3038 | 0.03/0.04 | C–H stretch | ||
ν8 | 3030 | 3030 | 3035/3032 | 0.25/0.07 | C–H stretch | ||
ν9 | 3027 | 3024 | 3031/3027 | 0.08/0.07 | C–H stretch | ||
ν10 | 3017 | 3023 | 3027/3017 | 0.03/0.50 | C–H stretch | ||
ν11 | 1770 | 1781 | 1768/1760 | 0.41/0.32 | 1714S | OC stretch | |
ν12 | 1759 | 1775 | 1756/1750 | 0.38/0.18 | 1712vs | OC stretch | |
ν13 | 1476 | 1478 | 1478/1474 | 0.45/0.75 | H–C–H bend | ||
ν14 | 1472 | 1478 | 1472/1470 | 0.74/0.65 | H–C–H bend | ||
ν15 | 1467 | 1468 | 1466/1465 | 0.70/0.45 | H–C–H bend | ||
ν16 | 1458 | 1466 | 1464/1463 | 0.74/0.68 | H–C–H bend | ||
ν17 | 1457 | 1458 | 1460/1456 | 0.71/0.74 | 1455vs | 1400w | H–C–H bend + H–C–C–C tors |
ν18 | 1452 | 1454 | 1445/1430 | 0.75/0.74 | H–C–H bend | ||
ν19 | 1390 | 1401 | 1401/1394 | 0.31/0.63 | H–C–H bend | ||
ν20 | 1388 | 1388 | 1392/1390 | 0.69/0.70 | 1353vw | 1365s | H–C–H bend |
ν21 | 1374 | 1384 | 1383/1375 | 0.63/0.72 | H–C–H bend | ||
ν22 | 1355 | 1333 | 1373/1370 | 0.71/0.64 | H–C–C bend | ||
ν23 | 1284 | 1301 | 1262/1244 | 0.60/0.36 | H–C–C bend + H–C–C–C tors | ||
ν24 | 1215 | 1198 | 1224/1221 | 0.70/0.73 | H–C–C bend | ||
ν25 | 1208 | 1188 | 1203/1197 | 0.68/0.63 | H–C–C bend | ||
ν26 | 1175 | 1172 | 1183/1182 | 0.70/0.74 | C–C stretch | ||
ν27 | 1094 | 1098 | 1096/1090 | 0.10/0.05 | O=C–C–C out-of-plane bend | ||
ν28 | 1055 | 1096 | 1059/1050 | 0.29/0.72 | C–C stretch | ||
ν29 | 1021 | 993 | 1033/1026 | 0.75/0.59 | H–C–C–C tors | ||
ν30 | 995 | 968 | 984/973 | 0.33/0.71 | C–C stretch + H–C–C–C tors | ||
ν31 | 966 | 956 | 966/958 | 0.19/0.17 | C–C stretch | ||
ν32 | 891 | 872 | 908/894 | 0.42/0.28 | H–C–C bend + H–C–C–C tors + C–C–C–C tors | ||
ν33 | 835 | 869 | 862/862 | 0.68/0.44 | C–C stretch | ||
ν34 | 768 | 799 | 761/752 | 0.74/0.08 | 829vs | C–C stretch | |
ν35 | 734 | 729 | 749/749 | 0.58/0.02 | H–C–C–C tors | ||
ν36 | 621 | 602 | 639/637 | 0.60/0.73 | C–C stretch + OC–C bend | ||
ν37 | 556 | 588 | 576/574 | 0.43/0.47 | O=C–C bend + C–C–C bend | ||
ν38 | 530 | 479 | 505/485 | 0.75/0.75 | OC–C bend + H–C–C–C tors + OC–C–C out-of-plane bend | ||
ν39 | 467 | 477 | 473/455 | 0.68/0.71 | H–C–C–C tors + C–C–C–C tors + OC–C–C out-of-plane bend | ||
ν40 | 439 | 448 | 452/429 | 0.74/0.71 | C–C–C bend + OC–C–C out-of-plane bend | ||
ν41 | 348 | 316 | 360/343 | 0.24/0.34 | C–C–C bend | ||
ν42 | 262 | 260 | 271/270 | 0.27/0.74 | C–C–C bend | ||
ν43 | 234 | 130 | 246/232 | 0.57/0.70 | C–C–C bend | ||
ν44 | 188 | 109 | 163/150 | 0.55/0.66 | C–C–C bend | ||
ν45 | 120 | 107 | 139/126 | 0.70/0.54 | H–C–C–C tors | ||
ν46 | 102 | 89 | 114/113 | 0.71/0.68 | H–C–C–C tors + C–C–C–C tors | ||
ν47 | 56 | 61 | 108/97 | 0.65/0.73 | C–C–C–C tors | ||
ν48 | 45 | 28 | 87/81 | 0.75/0.74 | C–C–C–C tors | ||
ν49 | 72/62 | 0.69/0.70 | Relative rotation | ||||
ν50 | 51/36 | 0.63/0.75 | Relative rotation | ||||
ν51 | 33/26 | 0.56/0.71 | Relative rotation | ||||
Sum of electronic and thermal free energies (Kcal mol−1) | Momomer | G = −385.130890 | |||||
Dimer | G = −770.264944 | ΔG = G(dimer) − 2G(momomer) = −0.003164 | |||||
Remarks | In the dimer, there are in-phase and out-of-phase vibrational modes, in-phase vibrational frequency is lower than the out-of-phase vibrational frequency |
Fig. 2 Comparison of FT-IR and LCM-Raman spectra with calculated Raman and IR spectra of acetonylacetone. |
The 36 atoms of the acetonylacetone dimer give rise to 102 normal modes of vibration. The overall 102 normal modes of vibration for this dimer are considered to comprise 96 normal modes arising from the in-phase and out-of-phase coupling of these two acetonylacetone molecules and 6 modes associated with the relative translation and rotation of two acetonylacetone molecules. A detailed description is listed in Table 1. The converted calculated Raman spectra using Multiwfn software are shown in Fig. 2. The corresponding in-phase and out-of-phase vibrational modes may differ in wavenumbers and depolarization ratios and the magnitude of these splitting will depend on the strength of interactions between different parts of the neighboring molecules. The overall agreement between the DFT calculated vibrational wavenumbers and the experimental values is good for acetonylacetone. Fig. 2 displays a comparison of the calculated monomer and dimer Raman spectra with the FT-Raman spectrum and FT-IR spectrum of acetonylacetone. The dashed lines in Fig. 2 indicate the correlation between the vibrational modes of acetonylacetone in the calculated Raman spectra and those corresponding to the fundamental modes of acetonylacetone in the FT-Raman and IR spectrum, respectively. The 1712 cm−1 band is assigned to the CO stretch. The wavenumber difference between the FT-Raman spectrum and FT-IR spectrum is 2.3 cm−1. This wavenumber difference is a key characteristic of the noncoincidence effect; other experimental proofs including the isotropic and anisotropic Raman spectra at different concentrations will be presented later. The largest difference between the calculated monomer and dimer lies in the CO stretching wavenumber. The wavenumber of the monomer is higher than that of the dimer. Other modes are similar with the calculated monomer and dimer structures. The calculation results show that when acetonylacetone transforms from the dimer to the monomer, the frequency of the CO stretching increases while other modes remain consistent. From the spectral comparison between the FT-Raman and the calculated monomer and dimer, the spectral pattern of the dimer describes more accurately the liquid acetonylacetone. It needs to be stressed that the listed calculated wavenumbers describe the vibrational frequencies of the molecule in its gaseous phase. Hence, the experimentally observed spectrum of the liquid phase may differ to some extent from the calculated spectrum. In the DFT calculation, the B3LYP function tends to overestimate the wavenumbers of the fundamental modes compared to the experimentally observed values due to the combination of electron correlation effects and basis set deficiencies. By comprehensively considering the calculated and experimental FT-Raman spectra, acetonylacetone is prone to present a short dimer order induced by CO vibrational td–td interactions. While in dilute solution, due to solvation effects and cage effects, where solvent molecules are surrounded by solute molecules that stabilize the solvate, acetonylacetone may present the monomer formation. Thus, the frequency of CO stretching is influenced by concentration effects.
To characterize the concentration effects, we collected the isotropic and anisotropic spectra of acetonylacetone in a variety of volume fractions, as shown in Fig. 3. They are assigned to the ν11(CO) stretching mode. The isotropic peak frequencies at 1710.7 cm−1 and the anisotropic peak at 1718.4 cm−1 were assigned to the calculated wavenumbers at 1760 cm−1 and 1768 cm−1, respectively. The corresponding depolarization ratios are 0.32 and 0.41. The dimer model calculations are in good agreement with the experimental non-coincident isotropic and anisotropic Raman data. The isotropic and anisotropic spectra, at a variety of volume fractions of acetonylacetone in CCl4, are shown in Fig. 3. Fig. 3 demonstrates that both isotropic and anisotropic Raman wavenumbers of the CO stretch increase with the dilution of acetonylacetone by CCl4, while the separation between isotropic and anisotropic Raman wavenumbers decrease from 7.66 cm−1 in neat acetonylacetone to 1 cm−1 at χm = 0.05. The FWHM (full width at half maxima) of the CO stretching modes also gets smaller and the peak gets sharper with decreasing acetonylacetone concentrations. The peak frequencies abstracted from Fig. 3 for isotropic (Iiso) and anisotropic (Ianiso) CO stretching Raman spectra are shown in Fig. 4. The Raman peak frequencies of both components show an increase in wavenumber with decreasing solute concentrations. The difference between isotropic (Iiso) and anisotropic (Ianiso) CO stretching wavenumbers ΔυNCE decreases upon dilution with CCl4 and reduces to 1.00 cm−1 at χm = 0.05, as shown in Fig. 5. This experimental data can be explained with our aggregation-induced split theory. Normally, acetonylacetone presents a dimer pattern. The dilution of the solute alters the short-range order of the CO stretching normal coordinate and the relative alignment of the dimer structure. The solvent molecules diffuse towards the reference molecule and break its structure, thereby weakening the dipole–dipole interactions of solute molecules. Thus, this leads to a decrease in the non-coincidence effect. The breaking of the dimer structure also makes the CO vibrational wavenumber shift to a higher wavenumber because the calculation shows that the CO stretching of the monomer lies at a higher wavenumber than that of dimers. This is in accordance with the concentration effect observed in experiments shown in Fig. 3 and 4.
Fig. 4 Variation of isotropic and anisotropic Raman peak frequencies of the CO stretching mode of acetonylacetone as a function of the solute concentration. |
Fig. 5 Variation in NCE of the CO stretching mode of acetonylacetone as a function of solute volume fraction. |
In our study, the NCE of the CO stretching mode of acetonylacetone is positive, and may be due to the antiparallel CO side-by-side organization. The antiparallel CO coupling split the CO vibrational wavenumber to two. It has been widely known that the NCE represents a spectroscopic manifestation of resonant intermolecular interactions between nearby IR-active oscillators through the transition dipole–transition dipole interaction mechanism.8–10 The cis and trans forms of acetonylacetone have CO Raman activity vibrational frequencies at 1770 cm−1 and 1781 cm−1, respectively. We also carried out IRC/path scan for the cis to trans transformation, as shown in Fig. S1.† This is a barrierless process. In the solution, normally, we believe that acetonylacetone predominantly bears the anti-form but this is not the case for the Raman spectra. Fig. 3 shows a split in the isotropic and anisotropic parts in most of the concentrated solutions. The concentration effect demonstrates that the aggregated structure is formed by intermolecular interactions, instead of intramolecular.
To study the influence of dipole moment of the solvent upon NCE, we collected the isotropic and anisotropic Raman spectra of acetonylacetone in a series of solvents with different static dielectric constant, as shown in Fig. 6. The corresponding NCEs were calculated and are illustrated in Fig. 7. Generally, the value of the NCE declined with an increase in the solvent dielectric constant with the same concentration. This rule is consistent with Logan's theory.11,35 Especially in water, as shown in Fig. 6, the value of the NCE is nearly equal to zero. From ref. 27 and 28, we speculate that it may be due to the formation of intermolecular hydrogen bonds between protons of water and the carbonyl groups of acetonylacetone, which hinder the dimerization of acetonylacetone molecules, thus making the NCE disappear.
Fig. 6 The isotropic and anisotropic parts of ν12(CO) vibrational Raman spectra of acetonylacetone in the binary mixture with different solvents (φA = 0.500). |
Fig. 7 Variation in NCE of the CO stretching mode of acetonylacetone as a function of the solvent dielectric constant. |
To investigate NCE solvent polarity dependent properties and demonstrate the rationality of the dimer model, the polarizable continuum model (PCM) was applied to calculate the dimer structure at the hybrid B3LYP-D3 levels of theory with the 6-311++G(d,p) basis set using the Gaussian 09 program. The solvent polarity influence, the optimized geometry, and the corresponding vibrational frequencies were obtained. Table S1† shows the DFT/PCM calculated CO vibrational frequencies, depolarization ratios, dielectric constants (ε), dipole moments (μ)/D, and ΔνNCE in a variety of solvents. With a decrease in the solvent dielectric constant, the two monomers of the acetonylacetone dimer became closer and the value of the NCE increased. These results are consistent with the experimental results shown in Fig. 7. Simultaneously, this verifies that a strong polar solvent will weaken the dimer structure of acetonylacetone in the mixture while a nonpolar solvent can reinforce the dimer structure. All results that are obtained from experiments, the theoretical dimer model, and DFT calculations demonstrate a consistent picture of the relationship between the NCE behavior, a spectroscopic feature of vibrational Raman bands, and the effect of dipolar interactions in liquid mixtures at molecular level.
Footnote |
† Electronic supplementary information (ESI) available: Additional FT-Raman spectra for acetonylacetone, DFT/PCM calculation of CO vibrational frequencies, depolarization ratios, dielectric constant (ε), dipole moments (μ)/D, and ΔνNCE for acetonylacetone. See DOI: 10.1039/d0ra02932g |
This journal is © The Royal Society of Chemistry 2020 |