Intermolecular dissociation energies of 1-naphthol complexes with large dispersion-energy donors: Decalins and adamantane

The ground-state intermolecular dissociation energies D₀(S₀) of supersonic-jet cooled intermolecular complexes of 1-naphthol (1NpOH) with the bi- and tricycloalkanes trans-decalin, cis-decalin, and adamantane were measured using the stimulated-emission-pumping/resonant two-photon ionization (SEP-R2PI) method. Using UV/UV holeburning, we identified two isomers (A and B) of the adamantane and trans-decalin complexes and four isomers (A-D) of the cis-decalin complex. For 1NpOH·adamantane A and B, the D₀(S₀) values are 21.6 ± 0.15 kJ/mol and 21.2 ± 0.32 kJ/mol, those of 1NpOH·trans-decalin A and B are 28.7 ± 0.3 kJ/mol and 28.1 ± 0.9 kJ/mol, and those of 1NpOH·cis-decalin A and B are 28.9 ± 0.15 kJ/mol and 28.7 ± 0.3 kJ/mol. Upon S₀ → S₁ electronic excitation of the 1NpOH moiety, the dissociation energies of adamantane, trans-decalin, and the cis-decalin isomer C change by <1% and those of cis-decalin isomers A, B, and D increase only slightly (1%-3%). This implies that the hydrocarbons are dispersively adsorbed to a naphthalene "face." Calculations using the dispersion-corrected density functional theory methods B97-D3 and B3LYP-D3 indeed predict that the stable structures have face geometries. The B97-D3 calculated D₀(S₀) values are within 1 kJ/mol of the experiment, while B3LYP-D3 predicts D₀ values that are 1.4-3.3 kJ/mol larger. Although adamantane has been recommended as a "dispersion-energy donor," the binding energies of the trans- and cis-decalin adducts to 1NpOH are 30% larger than that of adamantane. In fact, the D₀ value of 1NpOH·adamantane is close to that of 1NpOH·cyclohexane, reflecting the nearly identical contact layer between the two molecules.

Excitonic splittings in molecular dimers: why static ab initio calculations cannot match them

After decades of research on molecular excitons, only few molecular dimers are available on which exciton and vibronic coupling theories can be rigorously tested. In centrosymmetric H-bonded dimers consisting of identical (hetero)aromatic chromophores, the monomer electronic transition dipole moment vectors subtract or add, yielding S₀ → S₁ and S₀ → S₂ transitions that are symmetry-forbidden or -allowed, respectively. Symmetry breaking by ¹²C/¹³C or H/D isotopic substitution renders the forbidden transition weakly allowed. The excitonic coupling (Davydov splitting) can then be measured between the S₀ → S₁ and S₀ → S₂ vibrationless bands. We discuss the mass-specific excitonic spectra of five H-bonded dimers that are supersonically cooled to a few K and investigated using two-color resonant two-photon ionization spectroscopy. The excitonic splittings Δ_calc predicted by ab initio methods are 5-25 times larger than the experimental excitonic splittings Δ_exp. The purely electronic ab initio splittings need to be reduced ("quenched"), reflecting the coupling of the electronic transition to the optically active vibrations of the monomers. The so-called quenching factors Γ < 1 can be determined from experiment (Γ_exp) and/or calculation (Γ_calc). The vibronically quenched splittings Γ·Δ_calc are found to nicely reproduce the experimental exciton splittings.

Structure and intermolecular vibrations of perylene·trans-1,2-dichloroethene, a weak charge-transfer complex

The vibronic spectra of strong charge-transfer complexes are often congested or diffuse and therefore difficult to analyze. We present the spectra of the π-stacked complex perylene trans-1,2-dichloroethene, which is in the limit of weak charge transfer, the electronic excitation remaining largely confined to the perylene moiety. The complex is formed in a supersonic jet, and its S0 ↔ S1 spectra are investigated by two-color resonant two-photon ionization (2C-R2PI) and fluorescence spectroscopies. Under optimized conditions, vibrationally cold (T(vib) ≈ 9 K) and well resolved spectra are obtained. These are dominated by vibrational progressions in the “hindered-rotation” Rc intermolecular vibration with very low frequencies of 11 (S0) and 13 cm(–1) (S1). The intermolecular Tz stretch and the Ra and Rb bend vibrations are also observed. The normally symmetry-forbidden intramolecular 1a(u) “twisting” vibration of perylene also appears, showing that the π- stacking interaction deforms the perylene moiety, lowering its local symmetry from D2h to D2. We calculate the structure and vibrations of this complex using six different density functional theory (DFT) methods (CAM-B3LYP, BH&HLYP, B97-D3, ωB97X-D, M06, and M06-2X) and compare the results to those calculated by correlated wave function methods (SCS-MP2 and SCS-CC2). The structures and vibrational frequencies predicted with the CAM-B3LYP and BH&HLYP methods disagree with the other calculations and with experiment. The other four DFT and the ab initio methods all predict a π-stacked “centered” structure with nearly coplanar perylene and dichloroethene moieties and intermolecular binding energies of D(e) = −20.8 to −26.1 kJ/mol. The 000 band of the S0 → S1 transition is red-shifted by δν = −301 cm(–1) relative to that of perylene, implying that the D(e) increases by 3.6 kJ/mol or 15% upon electronic excitation. The intermolecular vibrational frequencies are assigned to the calculated Rc, Tz, Ra, and Rb vibrations by comparing to the observed/calculated frequencies and S0 ↔ S1 Franck–Condon factors. Of the three TD-DFT methods tested, the hybrid-meta-GGA functional M06-2X shows the best agreement with the experimental electronic transition energies, spectral shifts, and vibronic spectra, closely followed by the ωB97X-D functional, while the M06 functional gives inferior results.

Vibrational quenching of excitonic splittings in H-bonded molecular dimers: the electronic Davydov splittings cannot match experiment

The S(1)/S(2) state exciton splittings of symmetric doubly hydrogen-bonded gas-phase dimers provide spectroscopic benchmarks for the excited-state electronic couplings between UV chromophores. These have important implications for electronic energy transfer in multichromophoric systems ranging from photosynthetic light-harvesting antennae to photosynthetic reaction centers, conjugated polymers, molecular crystals, and nucleic acids. We provide laser spectroscopic data on the S(1)/S(2) excitonic splitting Δ(exp) of the doubly H-bonded o-cyanophenol (oCP) dimer and compare to the splittings of the dimers of (2-aminopyridine)(2), [(2AP)(2)], (2-pyridone)(2), [(2PY)(2)], (benzoic acid)(2), [(BZA)(2)], and (benzonitrile)(2), [(BN)(2)]. The experimental S(1)/S(2) excitonic splittings are Δ(exp) = 16.4 cm(-1) for (oCP)(2), 11.5 cm(-1) for (2AP)(2), 43.5 cm(-1) for (2PY)(2), and <1 cm(-1) for (BZA)(2). In contrast, the vertical S(1)/S(2) energy gaps Δ(calc) calculated by the approximate second-order coupled cluster (CC2) method for the same dimers are 10-40 times larger than the Δ(exp) values. The qualitative failure of this and other ab initio methods to reproduce the exciton splitting Δ(exp) arises from the Born-Oppenheimer (BO) approximation, which implicitly assumes the strong-coupling case and cannot be employed to evaluate excitonic splittings of systems that are in the weak-coupling limit. Given typical H-bond distances and oscillator strengths, the majority of H-bonded dimers lie in the weak-coupling limit. In this case, the monomer electronic-vibrational coupling upon electronic excitation must be accounted for; the excitonic splittings arise between the vibronic (and not the electronic) transitions. The discrepancy between the BO-based splittings Δ(calc) and the much smaller experimental Δ(exp) values is resolved by taking into account the quenching of the BO splitting by the intramolecular vibronic coupling in the monomer S(1) ← S(0) excitation. The vibrational quenching factors Γ for the five dimers (oCP)(2), (2AP)(2), (2AP)(2), (BN)(2), and (BZA)(2) lie in the range Γ = 0.03-0.2. The quenched excitonic splittings Γ[middle dot]Δ(calc) are found to be in very good agreement with the observed splittings Δ(exp). The vibrational quenching approach predicts reliable Δ(exp) values for the investigated dimers, confirms the importance of vibrational quenching of the electronic Davydov splittings, and provides a sound basis for predicting realistic exciton splittings in multichromophoric systems.

The S1/S2 exciton interaction in 2-pyridone·6-methyl-2-pyridone: Davydov splitting, vibronic coupling, and vibronic quenching

The excitonic splitting between the S(1) and S(2) electronic states of the doubly hydrogen-bonded dimer 2-pyridone[middle dot]6-methyl-2-pyridone (2PY·6M2PY) is studied in a supersonic jet, applying two-color resonant two-photon ionization (2C-R2PI), UV-UV depletion, and dispersed fluorescence spectroscopies. In contrast to the C(2h) symmetric (2-pyridone)(2) homodimer, in which the S(1) ← S(0) transition is symmetry-forbidden but the S(2) ← S(0) transition is allowed, the symmetry-breaking by the additional methyl group in 2PY·6M2PY leads to the appearance of both the S(1) and S(2) origins, which are separated by Δ(exp) = 154 cm(-1). When combined with the separation of the S(1) ← S(0) excitations of 6M2PY and 2PY, which is δ = 102 cm(-1), one obtains an S(1)/S(2) exciton coupling matrix element of V(AB, el) = 57 cm(-1) in a Frenkel-Davydov exciton model. The vibronic couplings in the S(1)/S(2) ← S(0) spectrum of 2PY·6M2PY are treated by the Fulton-Gouterman single-mode model. We consider independent couplings to the intramolecular 6a(') vibration and to the intermolecular σ(') stretch, and obtain a semi-quantitative fit to the observed spectrum. The dimensionless excitonic couplings are C(6a(')) = 0.15 and C(σ(')) = 0.05, which places this dimer in the weak-coupling limit. However, the S(1)/S(2) state exciton splittings Δ(calc) calculated by the configuration interaction singles method (CIS), time-dependent Hartree-Fock (TD-HF), and approximate second-order coupled-cluster method (CC2) are between 1100 and 1450 cm(-1), or seven to nine times larger than observed. These huge errors result from the neglect of the coupling to the optically active intra- and intermolecular vibrations of the dimer, which lead to vibronic quenching of the purely electronic excitonic splitting. For 2PY·6M2PY the electronic splitting is quenched by a factor of ~30 (i.e., the vibronic quenching factor is Γ(exp) = 0.035), which brings the calculated splittings into close agreement with the experimentally observed value. The 2C-R2PI and fluorescence spectra of the tautomeric species 2-hydroxypyridine·6-methyl-2-pyridone (2HP·6M2PY) are also observed and assigned.