D.
Bernhard
a,
M.
Fatima
b,
A.
Poblotzki
c,
A. L.
Steber
b,
C.
Pérez
b,
M. A.
Suhm
*c,
M.
Schnell
*b and
M.
Gerhards
*a
aTU Kaiserslautern, Fachbereich Chemie & Research Center Optimas, Erwin-Schroedinger-Str. 52, D-67663 Kaiserslautern, Germany. E-mail: gerhards@chemie.uni-kl.de
bDeutsches Elektronen-Synchrotron (DESY), Notkestr. 85, D-22607 Hamburg, Germany & Institute of Physical Chemistry, Christian-Albrechts-Universität zu Kiel, Max-Eyth-Str. 1, D-24118 Kiel, Germany. E-mail: melanie.schnell@desy.de
cInstitut für Physikalische Chemie, Universität Göttingen, Tammannstr. 6, D-37077 Göttingen, Germany. E-mail: msuhm@gwdg.de
First published on 11th July 2019
The structural preferences within a series of dibenzofuran–solvent complexes have been investigated by electronic, vibrational, and rotational spectroscopic methods probing supersonic jet expansions. The experimental study is accompanied by a detailed theoretical analysis including dispersion-corrected density functional theory, symmetry adapted perturbation theory, as well as coupled cluster approaches. The complementary, multi-spectroscopic results reveal a preferred OH⋯O structure for the smallest complex of dibenzofuran–water, whereas for the methanol complex an OH⋯π isomer is simultaneously observed. For the largest complex, dibenzofuran–tert-butyl alcohol, only a π-bound structure is found. These comprehensive investigations show that a completely inverse trend regarding the docking preference is observed by comparing the present results with the ones for analogous diphenyl ether complexes. This can be rationalized on the basis of the planarity/non-planarity and rigidity/flexibility of the different systems, providing valuable insight into the interplay between different non-covalent interactions. This analysis is a further step towards a quantitative description of very delicate energetic balances with the overall goal of yielding reliable structural predictions for non-covalently bound systems.
The structures of neutral, non-covalently bound molecular complexes involving aromatic moieties and water or alcohol molecules have been extensively studied in the gas phase (see e.g.ref. 7–9 and references therein). Among these studies, several complexes involving heteroaromatic moieties are found including works on indole–water,10 7-azaindole–water,11 pyrrole–12 and carbazole–solvent complexes.13,14 Moreover, several studies on furan derivatives were carried out, including a comparative FTIR jet and theoretical study on 2,5-dimethylfuran– as well as the 2,3-benzofuran–methanol complexes.15 The latter have additionally been studied by laser induced fluorescence and IR fluorescence dip spectroscopy, including the respective water complexes.16 For the 2,5-dimethylfuran–methanol complex, the OH⋯O binding motif was identified as the preferred structure, which exhibits additional CH⋯π stabilization by the interaction of the methyl group with the π-cloud. The OH⋯π-bound isomer was also observed as a slightly less stable structure, although it turned out to cause a stronger OH stretching red-shift than the OH⋯O-bound structure, i.e. the π-cloud causes a stronger distortion of the OH bond than the lone pairs of the ether oxygen. In the case of 2,3-benzofuran–water,16 a balanced situation was found with the coexisting, nearly isoenergetic OH⋯O and OH⋯π isomers. Within that study, no preferred site could be identified, which agrees with the theoretical predictions of less than 0.5 kJ mol−1 energy difference for the applied methods.16 However, the methanol complex of 2,3-benzofuran was found to prefer the π-docking site over the oxygen, which was clearly identified by linear FTIR jet spectroscopy.15 In this case, the π-bound structure exhibits a larger OH stretching red-shift than the OH⋯O structure, whereas theoretical predictions were ambiguous. The competing binding motifs in 2-methylfuran– and furan–methanol complexes are currently being examined in an experimental benchmarking study, with initial FTIR results already published,17 and a refined analysis with microwave support in preparation.
Several spectroscopic studies on dibenzofuran (DBF) have been performed in the condensed phase18–23 and in the gas phase,23–30 including works on the DBF dimer27,28,31 and mixed dimers of DBF with fluorene and benzofuran.28,32,33 Auty et al.24 recorded fluorescence excitation spectra of DBF and the DBF–water complex. They assumed that the complex is hydrogen-bonded to the oxygen atom of DBF based on the spectral shift of the fluorescence excitation bands. However, there is a lack of ab initio calculations supporting this assumption as well as further, structurally more sensitive spectroscopic experiments.
In previous studies, we established a multi-spectroscopic approach in order to elucidate the preferred binding sites in different aromatic ether–alcohol and –water complexes.34–38 Within the series of diphenyl ether (DPE) complexes, we have shown that water and small alcohol molecules prefer the π-docking site, whereas larger alcohols preferably bind to the ether oxygen atom. This observation contradicts the intuitive expectation of a preferred π docking, when the size of the alcohol increases. In that study, the respective contribution of London dispersion to the interaction energy for the different complexes was analyzed in order to explain the observed trend. In addition to that, the distortion of the DPE structure, caused by a twist of the phenyl rings upon aggregation of alcohol or water molecules, was identified as another major aspect influencing the trend.
In order to gain a deeper understanding of the influence of structural deformation upon complex formation, a systematic change in the structure of the ether is valuable. One possible change is the direct connection of the two phenyl rings of DPE, which leads to dibenzofuran (DBF). By doing so, the initial flexibility of DPE is entirely disabled since DBF is planar and rigid. Considering the fact that the π system is delocalized over both phenyl rings via the furan ring, the molecule is expected to remain planar upon solvent aggregation in order to maximize aromaticity. Therefore, there is no deformation of the ether geometry within the solvent complexes that might influence the binding preference, contrary to the case of DPE complexes in which deformation plays a substantial role.37 Moreover, the twisted structure of DPE was shown to enable CH⋯O contacts between ortho CH moieties and the oxygen atom of the alcohol or water molecule. This had a significant influence on the structural preference as well. Since for DBF only in-plane CH groups are available for CH⋯O contacts, structures mainly interacting via the π-cloud should not be affected by CH⋯O contacts. Whatever the relative importance of such qualitative concepts may be, they add up to a computable and experimentally verifiable energy difference between competing solvent docking sites.
Experimental verification of predicted structural preferences of such molecular complexes requires studies on a molecular level, where the isolated molecular aggregates can be investigated without the influence of any environment. These conditions can be fulfilled by molecular beam investigations, allowing the formation of molecular complexes and clusters in a supersonic expansion. A variety of spectroscopic methods can be combined with molecular beam experiments, including FTIR spectroscopy,15,34,35,38 mass- and isomer-selective IR/UV laser spectroscopy (IR/R2PI8,39–42) and chirped-pulse Fourier transform microwave (CP-FTMW43–45) spectroscopy. The combination of these different spectroscopic techniques yields complementary results, providing valuable experimental data ideally suited for benchmarking theoretical approaches.
In the present paper, we investigate a series of DBF complexes with water, methanol and tert-butyl alcohol (ROH with R = H, Me, t-Bu) by a multi-spectroscopic strategy, including FTIR, IR/UV and CP-FTMW spectroscopy. The experimental study is accompanied by a detailed theoretical analysis including dispersion-corrected density functional theory as well as wave function-based methods.
For the R2PI and IR/R2PI experiments, two tunable nanosecond laser systems were necessary, including one independent UV laser system and one IR laser system. The UV laser radiation is obtained via second harmonic generation in a BBO crystal using the output of a dye laser (Cobra-Stretch, Sirah). The latter is pumped by the second harmonic (532 nm) of a Nd:YAG laser (SpitLight 600, Innolas). The IR laser radiation in the range of 3520–3750 cm−1 is generated by difference frequency mixing (DFM) in a LiNbO3 crystal using the fundamental (1064 nm) of a Nd:YAG laser (Quanta-Ray Pro-230, Spectra-Physics) and the output of a second dye laser (PrecisionScan, Sirah), which is pumped by the second harmonic (532 nm) of the same Nd:YAG laser. The resulting IR radiation is amplified in an optical parametric amplification (OPA) process in another LiNbO3 crystal using the output of the DFM process and the fundamental (1064 nm) of the Nd:YAG laser. For the IR/R2PI spectra, the IR laser was irradiated 50 ns prior to the UV laser.
For each gas pulse, the ensemble of molecules was polarized with a series of eight microwave chirps with a 4 μs duration and spanning 2–8 GHz, following the fast-frame approach.50 The chirps were generated with an arbitrary waveform generator, amplified by a 300 W travelling wave tube amplifier, and transmitted into the vacuum chamber with a horn antenna. Following each excitation chirp, 40 μs of the free induction decay (FID) of the macroscopic ensemble of polarized molecules was recorded, yielding a frequency resolution of 25 kHz. For the experiments, a total of 5 million averages (for DBF–H2O) and 2 million averages (for DBF–MeOH and DBF–t-BuOH, respectively) were co-added and Fourier transformed with a Kaiser–Bessel window function to give the broadband rotational spectrum in the frequency domain.
All spectra were first fit to an asymmetric rotor Hamiltonian using the JB95 program.51 The transition frequencies were then refined using the AABS program suite, and the final asymmetric rotor Hamiltonian fits were completed with SPFIT.52 Line lists for all three dimers are provided in the ESI.† An analysis of the observed tunneling splitting arising from internal rotation of the methanol methyl group in the dibenzofuran–methanol complex was performed using the XIAM program.53 XIAM is a least squares fitting program specifically designed for analyzing spectra of molecules exhibiting internal rotors by employing the combined axis method of Woods to account for internal rotation through a potential barrier.
For the DFT calculations, the B3LYP functional58–60 including the D3 dispersion correction61 with Becke–Johnson (BJ) damping62 was used with the basis sets def2-TZVP,63 def2-QZVP63 and aug-cc-pVTZ64 (identical auxiliary basis sets for the RI approximation), while using the Berny optimization algorithm from Gaussian 0965 combined with energies and gradients calculated with Turbomole 7.3.66,67 Similarly, calculations were performed with the M06-2X functional68 including the D3 correction and the def2-TZVP basis set.
The SCS-CC2 calculations were carried out with the aug-cc-pVDZ64 and def2-TZVP basis sets using Turbomole 7.3, while correspondingly aug-cc-pVDZ-cbas69 and def2-TZVP-cbas69 were chosen as the auxiliary Coulomb fitting basis sets (cbas) required by the ricc2 module for the RI approximation. Harmonic vibrational frequencies at the SCS-CC2 level were calculated with the NumForce script of Turbomole 7.3.
All geometries were confirmed to be minima with only real harmonic vibrational frequencies. All DFT and SCS-CC2 energies were corrected for the basis set superposition error (BSSE) by applying the counterpoise correction method.70
DLPNO-CCSD(T) single-point calculations for the B3LYP-D3(BJ)/aug-cc-pVTZ geometries were carried out with ORCA 4.0.171 using the cc-pVTZ and cc-pVQZ basis sets64 with corresponding cc-pVTZ/C and cc-pVQZ auxiliary basis sets69 for the RI approximation. Additionally, the “TightPNO”72 and “TightSCF” options were applied. For comparison of zero-point-vibrational-energy (ZPE)-corrected energies, harmonic ZPE corrections obtained at the B3LYP-D3(BJ)/aug-cc-pVTZ level were added to the DLPNO-CCSD(T)/cc-pVTZ and DLPNO-CCSD(T)/cc-pVQZ energies. Furthermore, a local energy decomposition (LED) scheme73,74 was applied within the DLPNO-CCSD(T)/cc-pVQZ calculations. This was mainly used for extracting physically meaningful dispersion contributions to the total interaction energies. For comparison, second order SAPT(0) calculations75 were carried out with the truncated jun-cc-pVDZ basis set,64,76,77 using the PSI4 program.78
Transition state calculations were performed with the QST379 module of Gaussian 09 (B3LYP-D3(BJ)/def2-TZVP) with initial transition state and barrier guesses from the woelfling80 module of Turbomole 7.3 as well as the reaction path finder of the GFN-xTB program (option -path).
Fig. 1 Optimized minimum structures at the B3LYP-D3(BJ)/def2-TZVP level for DBF–ROH complexes (R = H, Me, t-Bu). Values in parentheses represent relative, ZPE-corrected energies in kJ mol−1 obtained at the DLPNO-CCSD(T)/cc-pVQZ level with ZPE from B3LYP-D3(BJ)/aug-cc-pVTZ calculations (see Table 1). Dashed lines qualitatively indicate different intermolecular contacts. |
Isomer | B3LYP-D3(BJ)/def2-TZVP | B3LYP-D3(BJ)/aug-cc-pVTZ | SCS-CC2/def2-TZVP | DLPNO-CCSD(T)/cc-pVQZa | |||||
---|---|---|---|---|---|---|---|---|---|
ΔE | ΔE0 | ΔE | ΔE0 | ΔE | ΔE0 | ΔE | ΔE0b | ||
a Geometries were used from B3LYP-D3(BJ)/aug-cc-pVTZ calculations. b Harmonic ZPE correction obtained from B3LYP-D3(BJ)/aug-cc-pVTZ calculations. | |||||||||
H2O | OH⋯Op | 0 | 0 | 0 | 0 | 0 | 1.58 | 0 | 0 |
OH⋯π5 | 2.42 | 1.18 | 2.92 | 1.88 | 0.58 | 0 | 3.38 | 2.34 | |
OH⋯π6 | 3.79 | 1.47 | 3.87 | 2.02 | — | — | 3.80 | 1.95 | |
MeOH | OH⋯Op | 0 | 0.57 | 0.39 | 1.06 | 1.62 | 3.63 | 0 | 0 |
OH⋯Ot | 1.33 | 0.96 | 0 | 0.05 | 0.88 | 1.88 | 0.85 | 0.21 | |
OH⋯π6 | 1.22 | 0 | 0.03 | 0 | 0 | 0 | 1.21 | 0.51 | |
t-BuOH | OH⋯Op | 1.66 | 1.57 | 0 | 0.21 | 0.60 | 1.38 | 0 | 0 |
OH⋯Ot | 0 | 0 | 0.16 | 0.84 | 0.64 | 0.85 | 0.98 | 1.46 | |
OH⋯π6 | 2.35 | 1.29 | 0.78 | 0 | 0 | 0 | 1.48 | 0.49 |
In order to analyze the respective contributions to the interaction energy of the investigated complexes, SAPT(0)/jun-cc-pVDZ calculations were performed (see Table S1, ESI†) as well as more sophisticated DLPNO-CCSD(T)/cc-pVQZ calculations for which a local energy decomposition (LED) scheme was applied (see Table S2, ESI†). As expected, both approaches yield a larger dispersion contribution in the OH⋯π motifs and significantly more electrostatic contribution for the OH⋯Op structure. This supports the finding that the OH⋯O hydrogen bond, dominated by electrostatics (see Table S1, ESI†), combined with the CH⋯O contact leads to a stronger stabilization than the two OH⋯π contacts within the other isomers.
For the DBF–MeOH complex, the calculated minimum structures at the B3LYP-D3(BJ)/def2-TZVP level are depicted in the second row of Fig. 1. In this case, two oxygen-bound isomers are found: within the Cs-symmetric OH⋯Ot isomer (on top), the methyl group of the methanol molecule is positioned above the furan ring, enabling CH⋯π interactions. The second oxygen-bound structure is denoted as OH⋯Op and exhibits an OH⋯O hydrogen bond in the DBF plane with the methyl group pointing away from the DBF plane. Therefore, the OH⋯Op isomer lacks CH⋯π interactions in contrast to the OH⋯Ot arrangement. However, the in-plane hydrogen bond allows for a stabilizing CH⋯O contact between the MeOH oxygen atom and a CH group of DBF (see Fig. 1), similar to the OH⋯Op isomer of DBF–H2O. The OH⋯π6 isomer is bound via an OH⋯π contact, and it is stabilized by CH⋯π interactions of the methyl group with the π-cloud. In contrast to the related systems 2,5-dimethylfuran–MeOH and 2,3-benzofuran–MeOH,15,16 no minimum structure is found with an OH⋯π interaction involving the five-membered furan ring. This might allow for the conclusion that within the extended π system of DBF, the six-membered benzene rings are better hydrogen bond acceptors than the furan moiety. As discussed in previous works on 2-methylfuran, 2,5-dimethylfuran, and 2,3-benzofuran,15–17 the furan oxygen acceptor site loses attractiveness upon the introduction of methyl groups or one phenyl moiety.
Analyzing the different energy contributions shows that the largest dispersion contribution is in the OH⋯π6 isomer, followed by the symmetric OH⋯Ot structure (see Tables S1 and S2, ESI†). Both arrangements contain CH⋯π stabilization. Accordingly, the dispersion contribution is significantly lower in OH⋯Op, and the structure is clearly dominated by the electrostatic contribution. The large difference in Eelst between the two OH⋯O structures might be explained by the almost coplanar arrangement of the hydrogen bond and the DBF plane in OH⋯Op in contrast to the clearly bent arrangement in OH⋯Ot.
In the latter case, this arrangement obviously reduces the electrostatic stabilization, while it simultaneously increases dispersion stabilization from CH⋯π interactions. Note that the OH⋯O hydrogen bonding angle itself is nearly identical in both structures (155°, Table S3, ESI†), indicating a clearly weakened hydrogen bond as it largely deviates from an ideal linear hydrogen bond. This significantly affects the OH stretching frequencies, which will be discussed in the Experimental results section.
The optimized minimum structures for the t-BuOH complex are shown in the last row of Fig. 1. Similar to the methanol complex, two OH⋯O structures and one OH⋯π arrangement are found as minimum geometries. The OH⋯Ot isomer is Cs-symmetric, identical to the OH⋯Ot isomer of the corresponding MeOH complex. The t-Bu moiety is positioned on top of the furan ring leading to CH⋯π interactions with the π-cloud. In contrast to the analogous MeOH complex, the hydrogen bond is less bent (177° for t-BuOH vs. 155° for MeOH, B3LYP-D3(BJ)/def2-TZVP level, see Table S3, ESI†).
In the OH⋯Op isomer, the t-BuOH moiety is tilted to one side, which indicates a slight CH⋯O interaction between the alcohol oxygen atom and a neighboring CH group, resembling a somehow distorted version of the OH⋯Op isomer of DBF–MeOH with the alcohol being located rather above the π plane due to stronger CH⋯π interactions. Comparing all non-Cs-symmetric OH⋯Op structures, the solvent molecule increasingly approaches the π-cloud above the DBF plane going from water to t-BuOH. In the respective OH⋯π6-bound isomer, the t-Bu moiety is in closer proximity to the DBF π-cloud, resulting in a larger interaction surface for CH⋯π interactions compared to the OH⋯Op isomer. However, CH⋯π interactions should be of similar magnitude in the OH⋯π and the Cs-symmetric OH⋯Ot isomer. This is in line with dispersion contributions obtained at the SAPT(0)/jun-cc-pVDZ and the DLPNO-CCSD(T)/cc-pVQZ levels, which are similar for the two isomers, but clearly smaller for the OH⋯Op isomer. As a hydrogen bond interaction is indicated by the OH stretching red-shift, calculated OH stretching frequencies can be compared for the competing structures. Calculations at the B3LYP-D3 level suggest a stronger OH⋯O hydrogen bond compared to OH⋯π. However, the contrary is predicted at the M06-2X/def2-TZVP and SCS-CC2 levels, indicating a stronger OH⋯π acceptor compared to OH⋯Ot. This aspect will be discussed later in the Experimental results section.
The calculated relative energies for all DBF–ROH complexes at different levels of theory are found in Table 1. The values for the B3LYP-D3(BJ) and SCS-CC2 levels result from geometry optimizations and harmonic frequency calculations, whereas single point calculations were performed at the DLPNO-CCSD(T)/cc-pVQZ level using the geometries obtained at the B3LYP-D3(BJ)/aug-cc-pVTZ level. No anharmonic treatments of the ZPE were used, as anharmonic corrections are expected to be small (assumed to be on the order of <0.5 kJ mol−1), and they have furthermore proven to perform non-systematically in relative energy predictions for similar systems.17 As shown in previous studies15 the structures with a rather localized OH⋯O hydrogen bond contain more ZPE than OH⋯π-bound structures. This is reflected in a consistent OH⋯Op destabilization on the order of 0.7–2.3 kJ mol−1 with respect to OH⋯π isomers when electronic (ΔE) and ZPE-corrected energies (ΔE0) are compared.
For the water complex, the prediction of the energetic order is almost uniform: the OH⋯Op structure is preferred by 0.6 up to 3.4 kJ mol−1, depending on the theoretical level and basis set. The highest applied level suggests the oxygen site to be preferred by about 2 kJ mol−1, which might raise questions about the population of a π-bound structure in molecular beam experiments unless major isomerization barriers prevent relaxation. The SCS-CC2 and M06-2X/def2-TZVP calculations (see Table S4, ESI†) prefer the OH⋯π5 structure. Regarding the DBF–MeOH complex, a rather undecided situation is found with an oscillation of the energetic order between an OH⋯π6, OH⋯Op and even OH⋯Ot preference. The ZPE destabilization of the oxygen-bound structures compared to the OH⋯π6 structure is even more relevant than in the case of water, since it switches the energetic order from a preferred OH⋯O isomer for pure electronic energies towards OH⋯π6 upon ZPE correction for B3LYP-D3(BJ)/def2-TZVP calculations. The most sophisticated theoretical approach at the DLPNO-CCSD(T)/cc-pVQZ level suggests an advantage for OH⋯Op of about 0.5 kJ mol−1 over OH⋯π6 and 0.2 kJ mol−1 over OH⋯Ot, whereas the B3LYP-D3(BJ)/aug-cc-pVTZ calculations even predict the OH⋯Ot and OH⋯π6 isomer to be isoenergetic (within 0.05 kJ mol−1) upon ZPE correction (see Table 1). Considering the energetic range of 0.5 kJ mol−1 for all three binding motifs obtained at the DLPNO-CCSD(T)/cc-pVQZ level – being certainly within the error bar of the method – would not exclude the presence of more than one isomer in molecular beam experiments. In the case of DBF–t-BuOH, the predicted binding preference is undecided as well among the different applied computational approaches: the symmetric OH⋯Ot isomer is preferred at the B3LYP-D3(BJ)/def2-TZVP level, whereas the larger aug-cc-pVTZ basis leads to an OH⋯π6 preference, together with the SCS-CC2 approach. Finally, the DLPNO-CCSD(T) approach favors the OH⋯Op structure by 0.5 kJ mol−1 over OH⋯π6. Overall, the relative ZPE-corrected energies of all three binding motifs are predicted to be within a range of 1.6 kJ mol−1. Thus, the simultaneous presence of more than one isomer cannot be excluded within supersonic jet experiments.
In order to elucidate the aspect of possibly co-existing isomers, being relevant for all investigated DBF–ROH complexes, the analysis of interconversion barriers can be helpful, aside from considering only the relative energies of the isomers. Therefore, transition state calculations were performed with the QST3 method as well as the woelfling module based on transition state guesses from the GFN-xTB method. The obtained interconversion barriers and transition state structures are shown in Fig. S1 (ESI†). In the case of the water complex, the calculated barrier of less than 1 kJ mol−1 between the two π-bound isomers OH⋯π5 and OH⋯π6 suggests that interconversion occurs under the supersonic expansion conditions. However, barriers of about 5 kJ mol−1 between the π-bound structures and the OH⋯Op isomer might allow a kinetic trapping of oxygen- and π-bound isomers, respectively, in the case where they are both initially populated.
For the MeOH complex, a low barrier of about 1 kJ mol−1 is predicted between the two oxygen-bound isomers OH⋯Ot and OH⋯Op, suggesting that interconversion occurs under the experimental conditions. Similar to the water complex, the isomerization barriers between OH⋯O and OH⋯π binding motifs are larger than the ones between the same binding motifs, yet they are slightly lower than that for DBF–H2O at about 3 kJ mol−1. Nevertheless, kinetic trapping of the respective lower energy isomer can be expected, in the case where more than one isomer is initially populated. Furthermore, the TS calculations suggest that the interconversion of the OH⋯Op structure into the OH⋯π6 isomer involves the Cs-symmetric OH⋯Ot structure as an intermediate state.
Regarding the tert-butyl alcohol complex, the barrier between the oxygen-bound isomers OH⋯Ot and OH⋯Op is calculated to be <1 kJ mol−1, suggesting interconversion. Similar to the case of DBF–MeOH, the TS calculations suggest that conversion of the OH⋯Op structure into the OH⋯π6 isomer occurs via the intermediate OH⋯Ot arrangement. The predicted isomerization barrier from the OH⋯Ot to the OH⋯π6 isomer is approximately 2 kJ mol−1. Hence, interconversion between the binding motifs should not be excluded as well. A discussion of these aspects with respect to the experimental findings will be continued in the Experimental results section.
For a comparison of theory and experiment, the structurally sensitive OH stretching vibration can serve as a spectroscopic probe to be compared to calculated harmonic OH stretching wavenumbers. In some cases, particularly if two competing structures with the same binding motif are present, the OH stretching vibrations might be indistinguishable. Therefore, the experimental rotational constants obtained from rotational spectroscopy combined with calculated dipole moment components can lead to an unambiguous structural assignment. All calculated values relevant for comparison to the experiments are found in Tables S5, S6 and S9–S11 (ESI†) and are discussed in the Experimental results section. In the end, comparison to the experiments will reveal the individual performance of each theoretical approach.
The R2PI spectrum of the DBF–H2O complex reveals a S1 ← S0 transition that is blue-shifted by +171 cm−1 compared to the DBF monomer transition (see Fig. S2, ESI†). This is in agreement with the findings of fluorescence excitation spectra.24 No additional, red-shifted transitions with respect to the monomer were detected. The experimentally observed shift of +171 cm−1 is in qualitative agreement with the predicted blue-shift of +174 cm−1 for the OH⋯Op isomer at the SCS-CC2/def2-TZVP level with respect to the DBF monomer transition (see Table S5, ESI†). In contrast, the OH⋯π6 isomer is expected to have a red-shifted S1 ← S0 transition compared to the DBF monomer (−138 cm−1; see Table S5, ESI†). Based on these findings, the assumption of an OH⋯O binding motif for the DBF–H2O complex based on the early non-mass-selective study by Auty et al.24 is confirmed by our mass-selective R2PI experiments combined with predicted S1 ← S0 excitation energies at the SCS-CC2/def2-TZVP level. These calculations have proven to yield robust predictions.57
In order to obtain additional structural information, the OH stretching vibration is used as a spectroscopic probe for identifying the docking motif of the complex. Therefore, an IR/R2PI spectrum was recorded via the electronic resonance at +171 cm−1, which is shown in Fig. 2 for the mass trace of the DBF–H2O complex compared to calculated harmonic OH stretching frequencies at the B3LYP-D3(BJ)/def2-TZVP level (additional computational results are found in Table S6, ESI†). The IR/R2PI spectrum obtained for the DBF–H2O complex (Fig. 2(a)) exhibits two OH stretching vibrations at 3623 and 3734 cm−1 (confirming preliminary observations in ref. 81). Calculations uniformly predict a red-shifted symmetric OH of the oxygen-bound structure compared to the π-bound structures. Mainly the symmetric OH is sensitive to the docking site, whereas the asymmetric OH is less affected. It is noticeable that the relative shifts of the symmetric OH between OH⋯O and OH⋯π structures are quite large for the B3LYP-D3(BJ) calculations in contrast to the ones for the M06-2X functional as well as the SCS-CC2 predictions. Upon comparison of symmetric OH frequencies to the experiment, the OH⋯O isomer matches better than the OH⋯π isomers for all applied theoretical approaches. In order to find a robust assignment, an additional aspect to be considered is the splitting between symmetric and asymmetric OH, which is significantly different for the OH⋯O and OH⋯π motifs (see Table S6, ESI†). The experimental splitting is found to be 111 cm−1, while the calculated splittings range from 105–131 cm−1 for the OH⋯Op isomer and only from 83 to 87 cm−1 for OH⋯π structures. Moreover, the changes of the splitting with respect to the splitting of free water clearly suggest the OH⋯Op isomer to be the observed structure. Note that the splittings obtained from the SCS-CC2 calculations are not considered here, as they are found to be unable to reproduce the frequency splitting of the free water molecule correctly.
Fig. 2 (a) IR/R2PI spectrum of DBF–H2O obtained via UV excitation at +171 cm−1 compared to calculated OH stretching frequencies (b) at the B3LYP-D3(BJ)/def2-TZVP level; scaling factor: 0.9677. |
Parametersa | Experiment | B3LYP-D3(BJ)/def2-TZVP | B3LYP-D3(BJ)/def2-TZVP |
---|---|---|---|
OH⋯Op | OH⋯Op | OH⋯π5 | |
a A, B, and C are the experimental rotational constants, ΔJ is the centrifugal distortion constant, and σ is the standard deviation of the fit. b Total number of fitted lines and as distributed over a-, b-, and c-type transitions. | |||
A (MHz) | 974.08 | 973 | 1339 |
B (MHz) | 575.23 | 581 | 479 |
C (MHz) | 362.00 | 365 | 447 |
Δ J (kHz) | 0.0334 | ||
σ (kHz) | 40 | ||
N lines (a/b/c) | 65(44/21/0) | ||
μ a/μb/μc (D) | 1.3/2.2/1.3 | 1.1/0.2/2.3 |
Fig. 3 FTIR spectra of DBF–MeOH (b–d) in comparison to 2,3-benzofuran–MeOH (a) taken from ref. 15. Spectrum (b) was measured using the double-slit nozzle (0.15% MeOH, T(DBF) = 120 °C), spectra (c) and (d) using the V-nozzle introducing MeOH via the saturator at temperatures of T(MeOH) = −25 °C (≈0.3%) and T(MeOH) = −10 °C (≈1%), respectively. The DBF concentration is slightly decreased in these spectra (T(DBF) = 110 °C). |
The two additional features marked with an asterisk (*) originate from ionization-induced fragmentation of a mixed DBF–MeOH–H2O cluster (see Fig. S3, ESI†). Given the overlap situation in the FTIR experiment and its relatively high nozzle temperature, as well as the different carrier gas, the wavenumber agreement between the two IR experiments is satisfactory. The observed isomer splitting of 5 cm−1 is probably more reliable than the 7 cm−1 deduced from the FTIR spectrum.
The relative order of calculated OH stretching wavenumbers for the different OH⋯O and OH⋯π structures turns out to be ambiguous: DFT calculations using the B3LYP-D3(BJ) functional predict a red-shifted OH stretch for both OH⋯O-bound structures compared to the OH⋯π6 isomer. Interestingly, calculations at the SCS-CC2 level suggest a switch of the order: the OH⋯π6 structure is predicted to have a lower OH stretching frequency than the OH⋯O equivalent, which suggests the π-cloud to be the stronger acceptor site. The same is observed for calculations with the functional M06-2X (see Table S6, ESI†). Note, however, that these two methods failed in predicting the correct complex with water.
Given the very small OH stretching frequency differences between the observed species, the prediction of the frequency order for a distinct theoretical method is ambiguous. Finally, based on the clear isomer assignment via the electronic resonances, the OH⋯Ot isomer is found to exhibit the less red-shifted OH stretching vibration compared to the OH⋯π6 isomer. This has been observed for similar systems as well (cf. discussion in the ESI† and ref. 15 and 16).
Drawing conclusions on relative populations from the electronic resonances in the R2PI spectrum is difficult in this specific case, as the R2PI spectrum of the DBF–MeOH mass trace is influenced by very strong resonances of the DBF monomer (for further explanation see Fig. S2, ESI†). Additional structural and abundance insight will be gained by rotational spectroscopy.
Fig. 5 Parts of the CP-FTMW broadband rotational spectrum, indicating the coexistence of the OH⋯π6 and the OH⋯Ot isomer. |
Parametersa | OH⋯π6 | OH⋯Ot | |
---|---|---|---|
Fit 1 A-states | Fit 2 XIAM | Fit 1 A-states | |
a A, B, and C are the rotational constants, ΔJ, ΔJK, ΔK, δJ and δK are the centrifugal distortion constants, and σ is the standard deviation of the fit. b V 3 is the barrier to internal rotation of the methanol methyl top, F0 is the corresponding rotational constant of the methyl top, and Dpi2J accounts for the internal motion-overall rotation centrifugal distortion, obtained from a global fit of both internal rotational components A and E with the program XIAM. c Total number of fitted lines and as distributed over a-, b-, and c-type transitions. | |||
A (MHz) | 987.3953(36) | 987.3736(27) | 808.91099(33) |
B (MHz) | 439.33213(28) | 439.32995(25) | 524.58247(15) |
C (MHz) | 417.58829(26) | 417.58775(21) | 375.30332(18) |
Δ J (kHz) | 0.0529(17) | 0.0554(13) | 0. 1434(15) |
Δ JK (kHz) | 0.738(10) | 0.7388(82) | 1.934(71) |
Δ K (kHz) | 56.91(80) | 56.85(75) | −0. 848(60) |
δ J (kHz) | — | — | 1.2416(45) |
F 0 (cm−1) | — | 159.0 | — |
V 3 (cm-1) | — | 338.99(89) | — |
V 3 (kJ mol−1) | — | 4.055(11) | — |
D pi2J (kHz) | — | 34.7(2.9) | — |
Δ (rad) | — | 3.883(13) | — |
ε (rad) | — | 0.579(23) | — |
# A state transitionsc | 65 (49/16/0) | 65 (49/16/0) | 122 (0/94/28) |
# E state transitions | — | 41 (29/12/0) | — |
σ (kHz) | 6.9 | 7.1 | 7.0 |
OH⋯Op | OH⋯Ot | OH⋯π6 | |
---|---|---|---|
A (MHz) | 737.4 | 831.9 | 1000.3 |
B (MHz) | 515.2 | 523.1 | 437.5 |
C (MHz) | 316.9 | 381.5 | 418.6 |
V 3 (kJ mol−1) | 3.8 | 3.1 | 5.1 |
μ a/μb/μc (D) | 1.9/1.6/1.3 | 0.01/3.1/0.5 | 1.1/1.2/0.2 |
Both complexes show internal rotation splitting due to the internal rotation of the methyl group of the methanol moiety, which results in characteristic doublets for each rotational line. The fact that this internal rotation leads to sizeable splittings and is not locked points to only a loose interaction of the methyl group with DBF. For OH⋯π6, two sets of rotational constants are presented. Fit 1 corresponds to a fit to an asymmetric rotor Hamiltonian including only the A states due to methyl group internal rotation, thus presenting effective rotational constants. Using the program XIAM, these line splittings can be analyzed, resulting in Fit 2. It includes the analysis of the methyl group internal rotation and thus also provides information about the torsional barrier V3. The V3 barrier determined from the experimental line splitting into A and E states is V3(exp) = 4.055(11) kJ mol−1 (Table 3), which is in decent agreement with the calculated barrier of about 5 kJ mol−1 (B3LYP-D3(BJ)/def2-TZVP). For OH⋯Ot, only an A-state fit is presented. An initial global fit including internal rotation via about 20 E-state lines resulted in standard deviations on the order of 100 kHz, potentially pointing to a second internal motion. The obtained rotational parameters of the A-state fit, however, allow a clear identification of the respective isomers.
The wealth of experimental information is thus well suited to identify and further characterize the respective molecules under study and can also be used to benchmark the theoretical models employed. Based on our experimental line intensities, the OH⋯π6 complex is found to be about 10 times stronger than the one for the OH⋯Ot complex. Considering the stronger dipole moment for the OH⋯Ot complex, this points to a clear energetic preference for the OH⋯π6 complex, which is also the global minimum at the B3LYP-D3(BJ)/def2-TZVP, B3LYP-D3(BJ)/aug-cc-pVTZ and the SCS-CC2/def2-TZVP levels. In the case of the B3LYP-D3(BJ)/aug-cc-pVTZ level, OH⋯π6 is nearly isoenergetic with the OH⋯Ot isomer. This finding agrees qualitatively with the FTIR evidence of a higher abundance of the more red-shifted species and the corresponding results for 2,3-benzofuran.15 It also agrees with the IR/UV experiment, considering the mentioned intensity uncertainty within the R2PI spectra. Interestingly, the third complex, OH⋯Op, was not observed in our microwave study despite intense analysis, even though it is predicted to be the global minimum by the DLPNO-CCSD(T) approach (cf.Table 1). Thus, a low interconversion barrier from the OH⋯Op complex to one or both of the other complexes, OH⋯Ot and OH⋯π6, can be assumed.
Fig. 6 FTIR spectrum of DBF–t-BuOH. The V-nozzle was used and the temperatures set to T(DBF) = 110 °C and T(t-BuOH) = 10 °C. |
As pointed out in the theoretical results section, the isomerization barrier between the OH⋯Op and the OH⋯π structure of the DBF-t-BuOH complex is calculated to be about 2 kJ mol−1, whereas the barrier between OH⋯Ot and OH⋯Op is smaller than 1 kJ mol−1. Since only one isomer is found in the experiment with neon as the carrier gas, it might be concluded that the isomerization barriers are too low for both binding motifs to be stabilized during the supersonic expansion. This would lead to the exclusive population of the global minimum structure in the molecular beam. Comparing the IR/UV and FTIR investigations (see Fig. 6 and 7), the shoulder at 3613 cm−1 exclusively observed in the FTIR spectrum seems to originate either from a less stable DBF–t-BuOH isomer, populated due to different expansion conditions, or from a larger cluster.
Fig. 8 Part of the CP-FTMW broadband rotational spectrum for DBF–t-BuOH, illustrating the presence of the OH⋯π6 isomer. |
Parametersa | Experiment | B3LYP-D3/def2-TZVP | B3LYP-D3/def2-TZVP | B3LYP-D3/def2-TZVP |
---|---|---|---|---|
OH⋯π6 | OH⋯π6 | OH⋯Ot | OH⋯Op | |
a A, B, and C are the rotational constants, ΔJ, ΔJK, ΔK and δJ are the centrifugal distortion constants, and σ is the standard deviation of the fit. b Total number of fitted lines and as distributed over a-, b-, and c-type transitions. | ||||
A (MHz) | 513.73023(20) | 515 | 453 | 473 |
B (MHz) | 351.76117(15) | 353 | 404 | 376 |
C (MHz) | 269.56523(14) | 269 | 263 | 246 |
Δ J (kHz) | 0.10324(73) | |||
Δ JK (kHz) | −0.3604(29) | |||
Δ K (kHz) | 0.2601(43) | |||
δ J (kHz) | 0.03650(54) | |||
σ (kHz) | 7.1 | |||
N lines (a/b/c) | 162 (83/36/43) | |||
μ a/μb/μc (D) | 1.3/0.1/0.4 | 2.2/0/2.0 | 2.1/0.6/1.7 |
For the DBF–H2O complex an oxygen bound structure was identified by electronic, vibrational and rotational spectroscopy, building on the early work of Auty et al.24 Despite the prediction of nearly isoenergetic π-bound structures, no second isomer is found. The DLPNO-CCSD(T)/cc-pVQZ method as well as B3LYP-D3(BJ)/aug-cc-pVTZ calculations yield reasonable relative energies that explain the experimental observations. Regarding the methanol complex, two isomers were identified in the molecular beam experiments. The species were identified as the OH⋯π6 isomer and the OH⋯Ot isomer. For the oxygen-bound structure, an interconversion of OH⋯Op to OH⋯Ot is expected due to a low isomerization barrier. Based on the broadband rotational spectroscopic results, the OH⋯π6 isomer is found to be more strongly populated than the OH⋯Ot structure, which is confirmed by the FTIR results and is also reasonable within the uncertainty of the R2PI signal intensities. Nevertheless, within the error of the methods, the predicted relative energies at the DLPNO-CCSD(T)/cc-pVQZ and B3LYP-D3(BJ)/aug-cc-pVTZ level are in agreement with the experimental findings since both theoretical methods indicate two nearly isoenergetic structures. Furthermore, the chosen theoretical approaches largely deviate in predicting OH stretching wavenumbers, which even leads to a switch in the order between the OH stretches of the two docking motifs. This has also been observed for the related 2,5-dimethylfuran–methanol15 and 2,3-benzofuran–methanol complexes.15,16 The only approaches that correctly predict a red-shifted OH stretching vibration for the OH⋯π6 isomer compared to OH⋯Ot are the M06-2X functional and calculations at the SCS-CC2 level. This probably indicates a deficiency of the established theoretical approaches including the harmonic approximation, which should be considered in future developments.
The tert-butyl alcohol complex, representing the largest ether–solvent complex in this study, was shown to form only one stable isomer in the molecular beam experiments. Based on a red-shifted S1 ← S0 transition, a π-bound structure was identified. Rotational spectroscopy clearly confirmed the observed structure to be the OH⋯π6 isomer.
By comparing all investigated DBF–solvent complexes, we observe a binding preference that switches from oxygen via a balanced situation to the π site when going from small solvent molecules to larger ones (see Fig. 9). This is inverse to the trend that has been observed for the previously studied series of diphenyl ether–solvent complexes. The stabilization due to London dispersion is found to be more pronounced in π-bound structures than in oxygen motifs, indicated by dispersion contributions extracted from both SAPT(0) and DLPNO-CCSD(T) calculations. An exception is found for the t-BuOH complex, where the Cs-symmetric OH⋯Ot and the OH⋯π6 isomers are found to have similar dispersion contributions. We moreover found the influence of CH⋯O contacts on OH⋯Op structures to decrease from the small solvent molecules to the larger ones, while simultaneously enabling stronger CH⋯π interactions in this series. However, their magnitude in the OH⋯O arrangements is always outweighed by the one in the respective OH⋯π structures. Thus, the additional CH⋯O contact leads to preferred oxygen-binding for the smaller solvent molecules, whereas for the larger t-BuOH, London dispersion finally outbalances the CH⋯O-stabilization, leading to preferred π-binding. Regarding relative energies, the overall performance of the B3LYP-D3(BJ)/aug-cc-pVTZ approach is satisfactory and it is, at a first glance surprisingly, even superior to the DLPNO-CCSD(T) approach that seems to slightly underestimate the stability of OH⋯π6 complexes compared to OH⋯Op. Given the fact that the relative energies are mostly below 1 kJ mol−1 and therefore in the order of ZPE corrections, this slight inconsistency might be attributed to the usually less relevant fact that the geometries are not optimized at the DLPNO-CCSD(T) level. Furthermore, neglected anharmonic contributions to the ZPE do not allow for a safe relative ranking of the two electronic structure methods at this subtle level.
Fig. 9 Schematic overview of the transition of binding motif preferences from oxygen-bound to π-bound structures with increasingly larger solvent molecules in DBF–ROH complexes. |
Upon comparison to the series of diphenyl ether–solvent complexes37 and also the phenyl vinyl ether–methanol complex,38 the absence of backbone deformation in the DBF complexes proves to be true, as no structures were found involving a non-planar DBF structure, which is not surprising since aromaticity is preserved. The influence of CH⋯O contacts – playing a decisive role in DPE and PVE complexes – partly remains in the DBF–solvent structures as well. However, it is constrained to the oxygen-bound motifs, as the CH moieties of DBF available for CH⋯O contacts are within the DBF plane and therefore not in vicinity to the alcohol oxygen atoms within OH⋯π arrangements. This finally leads to a trend regarding the binding preference that is contrary to the series of DPE–solvent complexes. A further interesting aspect is that the central, five-membered furan moiety in DBF is only a competitive π docking site for the water complex. In the MeOH and t-BuOH complexes, only structures involving the six-membered benzene moieties as π acceptors were observed.
One parameter, namely the rigidity of the aromatic chromophore (DBF vs. DPE), completely changes the preferred binding site within a series of solvent complexes. The big challenge is to quantify the low-temperature abundance by small energy differences between different structures. We succeeded by a comprehensive combination of different experimental and theoretical methods, which finally offers a clear structural assignment. Within the investigated series of DBF–solvent complexes, we found a variable interplay between non-covalent interactions among which London dispersion forces make an important difference in determining the final docking preference.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9cp02635e |
This journal is © the Owner Societies 2019 |