Exploring the Solvatochromism of Betaine 30 with Ab Initio Tools: From Accurate Gas-Phase Calculations to Implicit and Explicit Solvation Models

Benchmarking Charge-Transfer Excited States in TADF Emitters: ΔDFT Outperforms TD-DFT for Emission Energies

Charge-transfer (CT) excited states are crucial to organic light-emitting diodes (OLEDs), particularly to those based on thermally activated delayed fluorescence (TADF). However, accurately modeling CT states remains challenging, even with modern implementations of (time-dependent) density functional theory [(TD-)DFT], especially in a dielectric environment. To identify shortcomings and improve the methodology, we previously established the STGABS27 benchmark set with highly accurate experimental references for the adiabatic energy gap between the lowest singlet and triplet excited states (ΔEST). Here, we diversify this set to the STGABS27-EMS benchmark by including experimental emission energies (Eem) and use this new set to (re)-evaluate various DFT-based approaches. Surprisingly, these tests demonstrate that a state-specific (un)restricted open-shell Kohn-Sham (U/ROKS) DFT coupled with a polarizable continuum model for perturbative state-specific nonequilibrium solvation (ptSS-PCM) provides exceptional accuracy for predicting Eem over a wide range of density functionals. In contrast, the main workhorse of the field, Tamm-Dancoff-approximated TD-DFT (TDA-DFT) paired with the same ptSS-PCM, is distinctly less accurate and strongly functional-dependent. More importantly, while TDA-DFT requires the choice of two very different density functionals for good performance on either ΔEST or Eem, the time-independent U/ROKS/PCM approaches deliver excellent accuracy for both quantities with a wide variety of functionals.

The Simulation of Solvent Polarizabilities and Dipolarities with Polarizable Continuum Model

The ability of polarizable continuum models (PCM) to simulate nonspecific solvent effects (dipolarity and polarizability) was evaluated by calculating the transition energies of 1,1,10,10-tetrabutyldecanonaene (ttbp9) and 2-N,N-dimethylamino-7-nitrofluorene (DMANF), basis of Catalán's polarizability (SP) and dipolarity (SdP) solvent scales, respectively. Time-dependent density-functional theory (TD-DFT) calculations were performed at different levels of theory, employing four basis sets in 10 different solvents, covering the full range of the normalized SP and SdP scales. Transition energies were calculated using linear response (LR) and corrected linear response (cLR2) schemes. Although these methods yielded variable mean absolute errors, the LR-PCM calculations reproduced medium polarizability and dipolarity trends. While calculated ttbp9 transition energies correlated with SP and Laurence's dispersion-induced (DI) scales, the DMANF transition energies correlated poorly with SdP or Laurence's ES dipolarity scales. This result agrees with the fact that DMANF solvatochromism is "contaminated" by solvent polarizability and HB acidity. The incorporation of SP or DI contributions led to much better (r2 > 0.95) correlations with the DMANF-calculated transitions. The results offer a clearer picture of the limitations of continuum models in simulating the behavior of solvatochromic dyes in solution by pointing out their poor performance when specific solvent effects, such as hydrogen-bond interactions, play a significant role in their solvatochromism.

A Collection of Dispersion Induction DI, Electrostatic ES, and Hydrogen Bonding α1 and β1 Parameters for 380 Solvents and What They Say on Solvent Effects on Rates, Equilibria, and Spectra

The solvent parameters DI, ES, α1, and β1 are intended for the description of solute-solvent intermolecular forces, i.e., dispersion-induction, electrostatic, hydrogen-bond donation, and hydrogen-bond acceptance, respectively. An up-to-date collection of these parameters is presented for 380 solvents including green solvents, ionic liquids, and deep eutectic solvents. Their determination for additional solvents requires three commercial indicators, betaine dye B(30), 4-F-phenol, and 4-F-anisole (as well as the refractive index). The chemical significance of these parameters is outlined. Their use in the linear solvation energy relationship P = P° + di DI + e ES + a α1 + b β1 for 62 physicochemical properties P (reaction rates, equilibrium constants, and IR, UV, and NMR spectra) yields determination coefficients generally greater than 0.90 and regression coefficients whose sign and relative magnitude provide consistent information on the intermolecular forces acting on these properties.

Structural insight into piezo-solvatochromism of Reichardt's dye

To date, accurate modelling of the solvation process is challenging, often over-simplifying the solvent-solute interactions. The interplay between the molecular arrangement associated with the solvation process and crystal nucleation has been investigated by analysis of the piezo-solvatochromic behaviour of Reichardt's dye, ET(1), in methanol, ethanol and acetone under high pressure. High-pressure single-crystal X-ray diffraction and UV-Vis spectroscopy reveal the impact of solute-solvent interactions on the optical properties of ET(1). The study underscores the intricate relationship between solvent properties, molecular conformation and crystal packing. The connection between liquid and solid phases emphasizes the capabilities of high-pressure methods for expanding the field of crystal engineering. The high-pressure environment allowed the determination of the crystal structures reported here that are built from organic molecules fourfold solvated with ethanol or methanol: ET(1)·4CH3OH and ET(1)·4C2H5OH·H2O. The observed piezo-solvatochromic effects highlight the potential of ET(1) in nonlinear optoelectronics and expand the application of solvatochromic chemical indicators to pressure sensors.

Entropy-Enthalpy Compensation in Electron-Transfer Processes

Solvent reorganization energies, free energies, and entropies are obtained for photoexcitation of three molecules that exhibit strong solvatochromism [Nile red, 5-(dimethylamino)-5'-nitro-2,2-bisthiophene, and Reichardt's dye B30] by measuring their optical absorption spectra at temperatures between 150 and 300 K in solvents with a range of polarities. Energies, free energies, and entropies of solvent reorganization are also obtained from computer simulations of three intramolecular electron-transfer reactions (charge separation and recombination in Zn-porphyrin-quinone cyclophane and charge transfer in a bis-biphenylandrostane radical anion). Entropy-enthalpy compensation in the solvent reorganization free energy for photoexcitation or electron transfer is found to be essentially complete (having nearly equal and opposite contributions from entropic and enthalpic effects) for all of the processes with solvent reorganization energies less than about 0.1 eV. Compensation becomes less complete as the reorganization energy becomes larger. A semiclassical treatment of the solvent reorganization entropy can rationalize these results.

Temperature-dependent solvent reorganization entropies, free energies, and transition dipole strengths for the photoexcitation of Reichardt's dye B30

Absorption spectra of the solvatochromic dye 2,6-diphenyl-4-2,4,6-triphenyl-1-pyridinophenolate (B30) were measured in seven solvents of varying polarity over temperature ranging from each solvent's freezing point to 300 K. The excitation energies and their variances allowed calculations of the solvent reorganization energies, reorganization free energies and reorganization entropies as functions of temperature. The entropies of solvent packing around the chromophore are found to make major contributions to the reorganization free energies. The variances of the excitation energies depend only weakly on temperature, in disagreement with an expression that is often used for solvent reorganization free energies. Polar solvents reduce the transition dipole strength of B30's long-wavelength absorption band, probably because interactions with the solvent enhance the charge-transfer character of the transition. The dipole strength drops further at low temperatures.

Transient Terahertz Stark Effect: A Dynamic Probe of Electric Interactions in Polar Liquids

Electric forces acting on molecules in liquids at ambient temperature fluctuate at terahertz (THz) frequencies with a direct impact on their electronic and optical properties. We introduce the transient THz Stark effect to modify the electronic absorption spectra of dye molecules and, thus, elucidate and determine the underlying molecular interactions and dynamics. Picosecond electric fields of megavolts/cm induce a nonequilibrium response of the prototypical Betaine-30 in polar solution that is probed via transient absorption changes. The field-induced broadening of the absorption band follows the THz intensity in time, with a minor impact of solvent dynamics. The ground and excited state dipole energies in the THz field govern this response, allowing for a quantification of electric forces in a structurally frozen molecular environment.

The Global Polarity of Alcoholic Solvents and Water - Importance of the Collectively Acting Factors Density, Refractive Index and Hydrogen Bonding Forces

The DHBD quantity represents the hydroxyl group density of alcoholic solvents or water. DHBD is purely physically defined by the product of molar concentration of the solvent (N) and the factor Σn=n×f which reflects the number n and position (f-factor) of the alcoholic OH groups per molecule. Whether the hydroxyl group is either primary, secondary or tertiary is taken into account by f. Σn is clearly linearly correlated with the physical density or the refractive index of the alcohol derivative. Relationships of solvent-dependent UV/Vis absorption energies as ET (30) values, 129 Xe NMR shifts and kinetic data of 2-chloro-2-methylpropane solvolysis with DHBD are demonstrated. It can be shown that the ET (30) solvent parameter reflects the global polarity of the hydrogen bond network rather than specific H-bond acidity. Significant correlations of the log k1 rate constants of the solvolysis reaction of 2-chloro-2-methylpropane with DHBD show the physical reasoning of the approach.

The negative solvatochromism of Reichardt`s dye B30 - a complementary study

The UV/Vis spectra of a hypothetical negative solvatochromic dye in a solvent are theoretically calculated assuming the classical damped harmonic oscillator model and the Lorentz‐Lorenz relation. For the simulations, the oscillator strength of the solvent was varied, while for the solute all oscillator parameters were kept constant. As a result, a simple change of the oscillator strength of the solute can explain the redshift and intensity increase of the UV/Vis band of the solute. Simulated results are compared with measured UV/Vis spectroscopic data of 2,6‐diphenyl‐4‐(2,4,6‐triphenylpyridinium‐1‐yl) phenolate B30 (Reichardt‘s dye) Significant correlations of the absorption energy (1/λ_max) with the molar absorption coefficient ϵ as function of solvent polarity are demonstrated for several derivatives of B30. The approach presented is only applicable to negative solvatochromism. Therefore, while the approach is vital to fully understand solvatochromism, it needs to be complemented by other approaches, e. g., to describe the changes of the chemical interactions based on the nature of the solvent, to explain all its various aspects.

Theoretical, Semiempirical, and Experimental Solvatochromic Comparison Methods for the Construction of the α1 Scale of Hydrogen-Bond Donation of Solvents

Today, the hydrogen bonding donation (HBD) ability parameter of new solvents, α, is generally determined either by the Kamlet-Taft solvatochromic comparison of two probes, Reichardt betaine dye B(30) and 4-nitroanisole, or by the measurement of a single probe (e.g., solvatochromism of an iron coordination complex). This work highlights the shortcomings of these probes and recommends three replacement methods: (a) the theoretical comparison of the experimental and PCM-TD-DFT calculated transition energies E_T(30) of B(30), (b) the semiempirical comparison of the experimental and McRae calculated E_T(30), and, (c) for ionic liquids, the experimental comparison of E_T(30) and E_T(33) lying on the lower basicity of the betaine dye B(33) compared to B(30). These methods yield a new HBD parameter, α₁, for 101 molecular solvents and 30 ionic liquids. The novelty is emblematic for water, with α₁ = 1.54 instead of α (Kamlet-Taft) = 1.17. The solvent parameter α₁ is not equivalent to the solute hydrogen-bond acidity parameter α₂^H, partly because of the self-association of HBD solvents.

Solvation Effects in Organic Chemistry: A Short Historical Overview

This short overview describes the historical development of the physics and chemistry of organic solvents and solutions from the alchemist era until the present time based on some carefully selected examples that can be considered landmarks in the history of solution chemistry.

Micro- vs Macrosolvation in Reichardt's Dyes

Solvation is a complex phenomenon involving electrostatic and van der Waals forces as well as chemically more specific effects such as hydrogen bonding. To disentangle global solvent effects (macrosolvation) from local solvent effects (microsolvation), we studied the UV-vis and IR spectra of a solvatochromic pyridinium-N-phenolate dye (a derivative of Reichardt's dye) in rare gas matrices, in mixtures of argon and water, and in water ice. The π-π* transition of the betaine dye in the visible region and its C-O stretching vibration in the IR region are highly sensitive to solvent effects. By annealing argon matrices of the betaine dye doped with low concentrations of water, we were able to synthesize 1:1 water-dye complexes. Formation of hydrogen-bonded complexes leads to small shifts of the π-π* transition only, as long as the global polarity of the matrix environment does not change. In contrast, changes of the global polarity result in large spectral band shifts. Hydrogen-bonded complexes of the betaine dye are more sensitive to global polarity changes than the dye itself, explaining why E_T values determined with Reichardt's dyes are very different for protic and nonprotic solvents, even if the relative permittivities of these solvents are similar.

On the physical-chemical nature of solvent polarizability and dipolarity

The positive solvatochromism of three dyes, with a spectral behavior strongly dependents on the medium dipolarity/polarizability, was studied theoretically. Both a polarizable continuum-solvent model (CSM) and explicit solvent molecules were employed to model solvent effects. The CSM approach, coupled with ten different TDDFT methods, yielded unsatisfactory results in eleven solvents. The explicit-solvation calculations, thought of much higher computational cost, yielded excellent results. As CSM schemes are known correctly model non-specific electrostatic effects, our results indicate that the traditionally considered non-specific nature of solvent dipolarity needs to be reconsidered, requiring the explicit consideration of the solute-solvent interactions for their accurate theoretical description.

A computationally-derived model for the solvatochromism of p-phenolates with high predictive power

We report that the positive, reverse or negative solvatochromism of p-phenolate-based dyes is highly correlated with the multireferential (MR) character of their ground-state wave function, with negative compounds presenting the highest degeneracy. CASSCF/NEVPT2 calculations show that the high MR character of the wave-function in negative dyes allows those systems to increase the dipole moment of the ground state by breaking the degeneracy as a response to the field of a polar solvent. The resulting stabilization of the ground-state with increasing solvent polarity leads to the observed negative solvatochromic behavior. A computational indicator based on our results has been successfully used for determining the direction of the solvatochromic shift of 24 dyes. Thus, our work sheds light on the physical-chemical basis for solvatochromism while providing experimental chemists with a practical tool for the design of novel negative, positive or reverse solvatochromic dyes.

Solvation by Glycerol at Temperatures from 353 to 77 K: Its Solvatochromic Characterization and Use to Block the Molecular Structure of Conformationally Flexible Structures

Glycerol is UV/vis spectroscopically characterized by using suitable solvatochromic polarity probes spanning a wide temperature range of 353 to 77 K. For the first time we find experimental evidence that, when the solvent preserves its internal structure in a broad range of temperatures, all solvatochromically derived solvent parameters (i.e., SP, SdP, SA, and SB) also maintain their values unchanged. On the basis of these solvatochromic measurements, it is shown that, below 180 K, glycerol efficiently blocks the molecular structure of conformationally flexible solutes.