Pragya Verma,
Arnulf Rosspeintner and
Tatu Kumpulainen*
Department of Physical Chemistry, University of Geneva, 30 Quai Ernest Ansermet, Geneva, Switzerland. E-mail: tatu.kumpulainen@unige.ch; Fax: +41 22 379 6518; Tel: +41 22 379 6530
First published on 22nd June 2020
Characterization of propyl acetate/butyronitrile (PA/BuCN) mixtures by various spectroscopic techniques is described. The neat solvents have identical viscosities and refractive indices but their dielectric constants differ significantly. Detailed solvatochromic and titration data show that the mixtures do not exhibit specific solute–solvent interactions or significant dielectric enrichment effects. Therefore, the mixtures are ideally suited for investigating the effect of dielectric stabilization on (photo)chemical reactions. Dynamic Stokes shift experiments performed on two push–pull probes demonstrate that the solvation dynamics are significantly decelerated in the mixtures as compared to the neat solvents. Therefore, the mixtures allow for varying both the extent and time scale of the dielectric stabilization in a predictable manner.
Solvent mixtures of varying dielectric constant enable systematic studies on the influence of solvent stabilization on photochemical reactions.18–21 Ideally, other solvent parameters such as viscosity and refractive index should remain constant. In addition, the mixtures should not exhibit any specific hydrogen-bonding (HB) interactions or dielectric enrichment effects that might influence the reactivity. Specific solvent interactions can disrupt hydrogen bonds between the reaction partners or alter the solvent stabilization in a non-systematic way.22–26 On the other hand, solvents with a large difference in the dielectric constants (εr) are particularly prone to dielectric enrichment effects.27–29
A binary mixture between propyl acetate (PA) and butyronitrile (BuCN) is a promising candidate for investigating the effect of dielectric stabilization. The dielectric constant increases from 6.0 to 24.8 upon going from PA to BuCN and varies linearly with the weight fractions of the solvents according to:
εr = wPA·εr(PA) + (1 − wPA)·εr(BuCN), | (1) |
Solvent | η (cP) | εr | n | α | β | π* | μ (D) |
---|---|---|---|---|---|---|---|
a α, β, and π* are the Kamlet–Taft parameters related to HB donating, HB accepting ability, and dipolarity/polarizability respectively and μ is the dipole moment of the solvent. All parameters were taken from ref. 45. | |||||||
PA | 0.551 | 6.0 | 1.382 | 0.00 | 0.40 | 0.53 | 1.78 |
BuCN | 0.549 | 24.8 | 1.382 | 0.00 | 0.45 | 0.63 | 4.07 |
The article is structured as follow: first, we describe the steady-state solvatochromic behavior of betaine-30 (ref. 23) (B30), coumarin 153 (C153), and coumarin 152A (C152A) (Chart 1). The extremely large ground-state dipole moment of ca. 15 D (ref. 31) and the presence of the negatively charged phenolate oxygen render B30 absorption band frequency very sensitive towards both HB donors and dielectric enrichment effects.32,33 On the other hand, B30 is nearly non-fluorescent.34–36 Therefore, the two coumarins were used and their solvatochromism was investigated by using both steady-state absorption and emission. Additionally, we present a detailed investigation on the solvent relaxation parameters utilizing the time-dependent Stokes shift method37 on broadband fluorescence spectra with a sub-100 fs time-resolution.38,39 Time constants for the solvent relaxation were determined for both C153 and C152A in the mixtures. In the third part, we present a detailed investigation on the microscopic viscosity by investigating both the rotational and translational diffusion of C153 in the mixtures. Both of these coefficients depend on the viscosity and hence should remain constant in isoviscous mixtures in the absence of specific solute–solvent interactions.40–42
Chart 1 Chemical structures of (A) solvatochromic probes, (B) photoacid and organic base, and (C) solvents used in the study. Same abbreviations are used throughout the text. |
Since our main goal is to use the mixtures for proton transfer reactions, we demonstrate in the last part that the mixture does not disrupt the HB interactions between the reaction partners and that the dielectric constant has a significant impact on the driving force of the proton-transfer reaction. For this purpose, we use a 1,8-naphthalimide-based photoacid (Chart 1) complexed with either a weak base (N-methylimidazole, NMI) or deprotonated by a strong base (1,8-diazabicyclo[5.4.0]undec-7-ene, DBU).43,44
(2) |
(3) |
ET(30) [kcal mol−1] = hc0NA = 2.8591·0 [cm−1], | (4) |
The transition energies of the coumarins were determined from the peak frequencies of the log-normal fits in transition dipole moment (TDM) representation since this corresponds to the true lineshape of the underlying transitions.51 This is supported by the bandwidths of the absorption and emission spectra, presented in Fig. S7 (ESI†). The absorption and emission bandwidths are nearly identical in the TDM representation whereas they differ by more than 500 cm−1 in the wavenumber representation.
All transition energies and Stokes shifts are represented as a function of the reaction field factor, Δf, defined as:
(5) |
The factor, derived from Onsager's reaction-field theory, accounts for the dielectric stabilization of a non-polarizable point dipole in a spherical cavity due orientational dipole–dipole interactions.55 Since the refractive index of the mixtures is constant, the transition energies are expected to show linear behavior with respect to Δf with slopes proportional to a combination of the ground- and excited-state permanent dipole moments.56,57 Any deviations from linearity would indicate the presence of specific solute–solvent interactions or dielectric enrichment effects. For consistency, all solvent parameters are plotted against Δf and the weight fractions of PA were chosen to be equally spaced in Δf. Relation between Δf and the weight fraction (Fig. S4, ESI†) is nonlinear.
Rotational diffusion times of C153 were determined from the fluorescence anisotropy decays. These decays were calculated from intensities measured at a parallel (I∥) and perpendicular (I⊥) pump polarization with respect to the detected signal according to:
(6) |
The rotational times were obtained by fitting the anisotropy decays with a bi-exponential function. The goodness of fit was judged based on the weighted residual and the χr2 value. Translational diffusion coefficients of C153 were determined from the 19F-DOSY NMR experiments (see ESI† for details).
Solvation dynamics were determined from the fs-broadband fluorescence spectra using the time-dependent Stokes shift method.37 The spectra were analyzed globally using the approach discussed in detail in ref. 58. In brief, the fluorescence band is fitted with a time-dependent log-normal function where the dynamic Stokes shift is expected to follow a multi-exponential function.
Fig. 1 (A) Absorption spectra of B30 in the PA/BuCN solvent mixtures. (B) Corresponding ET(30) values as a function of Δf. The black line represents ideal behavior connecting the two extreme points corresponding to the neat solvents and deviation of the experimental values is shown in the top panel in kcal mol−1. The dashed gray line represents a linear regression line to all points. Red crosses indicate the ET(30) values of the neat solvents taken from ref. 45. |
The ET(30) values show a good linear correlation as a function Δf (R2 = 0.997) represented by the dashed gray line in Fig. 1B. Minor deviations are, however, visible at large weight fractions of PA (Δf < 0.5). In the ideal case, the ET(30) values would vary linearly between the extreme points corresponding to the neat solvents represented by the black solid line in Fig. 1B. The solvation energy increases more rapidly in PA rich mixtures but returns to ideal behavior at around Δf = 0.52 corresponding to a nearly 50/50 mixture (wPA = 0.56). Deviation of the experimental values from the ideal behavior is presented in the top panel of Fig. 1B. It should be noted, that B30 is particularly susceptible to both dielectric enrichment as well HB donor interactions. We observed significant deviations from linearity before purification of the solvents. Therefore, the purity of the solvents was confirmed by comparing the measured ET(30) values of the neat solvent to the literature values (red crosses in Fig. 1B) prior any further measurements. After purification the agreement between the measured and literature values was around 0.1 kcal mol−1. Therefore, the observed deviation from linearity is unlikely due to any impurities but can be attributed to the dielectric enrichment. However, the maximum deviation of 0.25 kcal mol−1 at Δf = 0.48 represents only 5% of the total variation and overall linear correlation between the solvation energy and Δf is excellent.
The peak frequencies of the absorption and emission spectra of C153 and C152A were determined in the TDM representation and the Stokes shifts were calculated as a difference between the absorption and emission frequencies. Experimental spectra and spectral parameters are given in ESI.† The peak frequencies and Stokes shifts as a difference to the values determined in neat PA are plotted as a function of Δf in Fig. 2. The negative slopes of the peak frequencies indicate red shifts of the bands whereas the positive slopes of the Stokes shifts indicate an increase.
For both coumarins, the linear correlation between the transition energies and Δf is reasonably good. C153 shows slight scattering of the points, which are overemphasized in the Stokes shifts. Deviations from the linear regression lines are shown in the top panels. However, this largely originates from uncertainties in the fitted peak maxima due to the log-normal function and the deviations do not show any clear patterns. For C153 in particular, the log-normal function is a sub-optimal band shape function in solvents of medium polarity. For C152A, the log-normal function reproduces the spectra better and the spectral shifts are larger resulting in an excellent linear correlation. Therefore, the observed solvatochromism of both compounds can be attributed to non-specific dielectric stabilization.
In summary, dielectric stabilization in the mixtures follows nearly ideal behavior for all compounds investigated. Although slight dielectric enrichment was observed for B30, the linear correlation between the solvation energies and Δf is remarkably good. Therefore, the mixture is ideally suited for controlling the extent of the dielectric stabilization.
Δf | α1 | τ1/ps | α2 | τ2/ps | α3 | τ3/ps | 〈τsolv.〉/ps |
---|---|---|---|---|---|---|---|
a Average solvation time was calculated according to 〈τsolv.〉 = Σαiτi. | |||||||
0.39 | 0.54 | 0.25 | 0.37 | 2.1 | 0.09 | 15.8 | 2.3 ± 0.2 |
0.43 | 0.54 | 0.24 | 0.37 | 2.3 | 0.09 | 25.2 | 3.2 ± 0.2 |
0.48 | 0.53 | 0.25 | 0.39 | 2.1 | 0.08 | 36.1 | 3.9 ± 0.3 |
0.52 | 0.55 | 0.24 | 0.38 | 2.0 | 0.07 | 27.4 | 2.7 ± 0.1 |
0.56 | 0.61 | 0.26 | 0.39 | 1.6 | — | — | 0.8 ± 0.1 |
Similarly to the steady-state spectra, magnitudes of the time-dependent Stokes shifts increase upon increasing Δf, C152A showing slightly larger shifts compared to C153. The relaxed fluorescence band maxima resulting from the global fit decrease linearly with Δf and are within 100 cm−1 from the peak maxima of the steady-state spectra. However, the band maxima of the initial time-zero spectra (t=0) as well as the total shifts (Δtot) are significantly overestimated in the global fits for both compounds. This is largely due to the non-exponential behavior of the spectral shifts at t < 150 fs that is sometimes modeled as a Gaussian decay.26,59,60 In this region, the exponential fit significantly deviates from the experimental band maxima as can be seen in Fig. 3A. Therefore, the total shifts were additionally calculated using the time-zero spectra estimated from the steady-state spectra in non-polar solvent according to a previously published procedure (see ESI† for details and representative spectra).58,61
The peak frequencies of the time-zero and relaxed fluorescence estimated from the steady-state spectra (black markers) and those resulting from the global fit (red markers) are presented in Fig. 3B. The total shifts, reflecting the solvation free energy,62 increase linearly with Δf in the estimated spectra whereas the shifts resulting from the global fit are overestimated by ca. 200–300 cm−1. However, the experimental band positions (markers in Fig. 3A) at t = 0 agree within 100 cm−1 with the estimated time-zero spectra. This further supports that the dielectric stabilization follows ideal behavior and is not complicated by specific solute–solvent interactions or dielectric enrichment effects.62
Despite the different magnitude of the total shifts, the time constants and relative amplitudes of the two coumarins are very similar. This demonstrates that the solvation dynamics are independent of the coumarin, in agreement with previous results.58 Therefore, the characteristic solvent relaxation parameters, presented in Table 2, are given as an average of the individual parameters obtained for C153 and C152A. The mean solvation times correspond to amplitude-weighted average lifetimes and are presented in Fig. 3C as a function of Δf.
Surprisingly, the solvation dynamics are significantly decelerated in the mixtures. Due to this unexpected behavior, the experiments were performed in duplicate reproducing the same behavior. The solvent relaxation becomes gradually slower upon increasing the amount of BuCN reaching a maximum at Δf = 0.48 corresponding to wPA = 0.79. The first two lifetimes do not show significant differences, but the third time constant increases significantly in the mixtures. After this point, the solvation dynamics become faster upon increasing fraction of BuCN. In neat BuCN, the solvation dynamics are the fastest and slowest decay component (τ3) disappears. The physical origin of this effect remains unclear and is beyond the scope of the present study. However, similar deceleration has been observed in other binary mixtures where it has been attributed to deceleration of the intrinsic solvent reorientation dynamics or translational components.18,63,64 The rather long time scale of the deceleration observed in PA/BuCN mixtures points towards translational components that likely involve a significant reorientation of the first solvation shell.
The rotational diffusion coefficient is proportional to the inverse of the rotational correlation time (τrot.). For a spherical molecule, τrot. is directly proportional to the viscosity of solvent according to the Stokes–Einstein–Debye equation:41,42
(7) |
(8) |
The rotational correlation times of C153 were determined from the fluorescence anisotropy decays in the PA/BuCN mixtures.42 Representative fluorescence decays measured in parallel and perpendicular polarizations together with the constructed anisotropy decay are shown in Fig. S12 (ESI†). The anisotropy decays were analyzed with a bi-exponential function and the best-fit parameters are summarized in Table S9 (ESI†). The bi-exponential decay originates from the non-spherical shape of C153 (ref. 65) but both individual decay times are expected to show linear correlation with viscosity, as has been observed in other binary mixtures.16 The mean rotational correlation times, 〈τrot.〉, were calculated as an amplitude-weighted average lifetime and are presented as a function of Δf in Fig. 4A. The translational diffusion coefficients of C153 were determined from the 19F-DOSY experiments in the PA/BuCN mixtures. Experimental details and representative data are given in ESI.† The diffusion coefficients are presented as a function of Δf in Fig. 4B.
Fig. 4 (A) Rotational correlation times and (B) translational diffusion coefficients of C153 in the PA/BuCN mixtures as a function of Δf. The lines indicate the mean values in the mixtures. |
Both the rotational correlation times and translational diffusion coefficients are nearly constant in the mixtures. This unambiguously demonstrates that the viscosity of the mixtures remains constant. The invariant rotational correlation times are additionally a further indication for the absence of strong solute–solvent interactions. Both specific interactions and dielectric enrichment has been shown to influence the rotational correlation times.42,66 Our results show that the mixture behaves ideally and does not show any deviations from hydrodynamic theories.
Since C4-dHONI has two binding sites, a model based on 1:2 host–guest stoichiometry was used. The association constants, K1 and K2, were obtained from the steady-state absorption spectra using a previously described global analysis procedure (see ESI† for details).44 Representative absorption spectra of C4-dHONI upon addition of NMI and DBU in a PA/BuCN mixture with Δf = 0.48 are presented in Fig. 5A and C and the association or deprotonation constants as a function of Δf in all the mixtures in Fig. 5B and D.
The absorption spectra of C4-dHONI shows a broadening and a red shift upon addition of NMI. We have previously shown that NMI forms neutral hydrogen-bonded 1:1 and 1:2 complexes with C4-dHONI.44 The association constants, K1 and K2, determined from the global analysis remain constant in the mixtures. This clearly demonstrates that the solvent mixtures do not influence the HB equilibrium between the reaction partners. Therefore, the mixtures could be additionally used to investigate the effect of dielectric constant on different HB association reactions.
The absorption spectra upon addition of DBU show a significant decrease of the neutral band around 380 nm and an appearance of a long-wavelength band at 450 nm. This has been attributed to a deprotonation of a single hydroxyl group by the strong base followed by HB association of another base molecule on the remaining hydroxyl group of C4-dHONI. This is due to the fact that deprotonation of the first hydroxyl group significantly increases the pKa of the second. The deprotonation constants, K1, increase drastically upon increasing Δf, being 30 times higher in BuCN compared to PA. This shows that the reaction free energy for proton transfer reactions significantly decreases (more negative ΔG) upon increasing Δf. On the other hand, the association constant for the binding of a second DBU, K2, remains nearly constant in neat PA and the mixtures but increases by a factor of four neat BuCN. This further demonstrates that the mixtures do not significantly influence HB equilibrium. The increase in BuCN might originate from incomplete description of the experimental spectra by the association model because different ion pair species (contact and non-contact) can further complicate the spectral behavior in polar solvents.44 In any case, the impact of the dielectric stabilization on the ground-state deprotonation constant is very clear. The next step is to investigate the impact on photoinduced excited-state proton-transfer reactions and work in this direction is under way in our laboratory.
Footnote |
† Electronic supplementary information (ESI) available: Additional experimental details, characterization of the macroscopic solvent properties, supplementary experimental data, fits and parameters. See DOI: 10.1039/d0ra04525j |
This journal is © The Royal Society of Chemistry 2020 |