Light–matter strong coupling is evolving as a physical tool for chemists to control chemical reactivity , , , , , , , , , . This was triggered by the modified spiropyran–merocyanine photoisomerisation dynamics under the electronic strong coupling (ESC)  and further fuelled by the development of vibrational strong coupling (VSC), where the vibrational mode of a molecule is strongly coupled , , , , , , , . Vibrational strong coupling was envisioned to alter the rates and yields of ground state chemical reactions as the latter involves the scission and formation of chemical bonds and therefore vibrations . The retarded kinetics of the desilylation of 1-phenyl-2-trimethylsilylacetylene (PTA) provided the first evidence for a VSC-controlled chemical reaction . Subsequent experiments under VSC showed that the reactions can be accelerated , , and the chemical landscape could be modified in such a way that the product branching ratio is altered . Remarkably, the reactivity can also be controlled by the strong coupling of solvent vibrational modes , . Various theoretical models have been developed recently, which support the modified chemical dynamics under strong coupling , , , , , , , . However, the mechanistic details of chemical reactions under VSC are not well understood, although a clear change in the activation barrier is observed in experimental , , ,  and in theoretical studies , , . The present study aims to provide further insights into the transition state (TS) changes under VSC by monitoring the thermodynamic parameters of the desilylation reaction of PTA as a function of the strong coupling strength.
Light–matter strong coupling, for example VSC, occurs when the resonance energy exchange between the molecular vibration and the eigen mode of an electromagnetic cavity (Fabry–Pérot cavity) proceeds faster than any other decay process . This results in hybrid light–matter states, called the polaritonic states (P+ and P−), which are visible as split vibrational modes separated by the Rabi splitting energy (ħΩR) in the absorption spectrum of the molecule. The formation of polaritonic states is schematically shown in Figure 1A. The coupling strength, ħΩR, is often decisive in the modification of molecular properties, and its magnitude for a molecule can be calculated with the Jaynes–Cummings two-state model, which reduces to (1) in the absence of any dissipation:
where d is the transition dipole moment of the molecule (which is determined by the derivative of the permanent dipole moment of the vibrational mode), ħω is the optical resonance, ε0 is the vacuum permittivity, ν is the volume of the electromagnetic mode, and nph gives the number photons involved in the strong coupling. Detailed description of strong light–matter coupling, and its impact on molecular and material properties can be found in , , , , , , , . For discussions specific to the physics of VSC, the reader is referred to  and . The last part of (1) points out that ħΩR has a residual value known as vacuum Rabi splitting energy even in the absence of real photons. In other words, the zero-point energy fluctuations of the cavity enable the formation of light–matter hybrid states in the dark. The experiments reported here are carried out under this regime, such that the weak IR light used to intermittently probe the system does not increase the ħΩR in this collective coupling condition. The magnitude of ħΩR is varied only when the number of molecules (N) involved in the coupling process is changed. This is because the volume of a single electromagnetic mode can accommodate a large number of molecules. The wave functions of the hybrid states are delocalised over all the molecules interacting with the cavity mode, and therefore the where C is the concentration of the molecule. This scenario is schematically shown in Figure 1A with VSC of PTA as an example.
2 Results and discussion
The desilylation of PTA (Supplementary Material, Figure S1) was chosen as the model reaction for the present study. The VSC of the Si–C stretching vibrational mode (νSi−C; 860 cm−1, blue curve in Figure 1B) of PTA has shown to retard the Si–C bond scission dynamics , prompting us to do a detailed thermodynamic study, varying the ħΩR. Vibrational strong coupling of the νSi−C of PTA is achieved using a Fabry–Pérot cavity assembled in a temperature-controllable and demountable IR cell (purchased from Specac) equipped with holes to inject the liquids. An air-filled Fabry–Pérot cavity (empty cavity) is built in the IR cell by placing a Mylar spacer (6±1 μm thick with a hollow central region) on top of an Au-coated ZnSe window (bottom mirror), which is then covered by a top mirror of Au sputtered on a ZnSe substrate having an inlet and outlet hole for sample injection. To protect the Au from the reactants, a 200-nm-thick SiO2 (glass) was deposited on top of the Au mirrors. The IR cell is then closed by tightening the screws on the top lid. Such an empty cavity has many equally spaced modes in the IR whose periodicity is determined by the path length of the cavity (Figure S2). Using a screwdriver, the empty cavity path length is then fine-tuned by adjusting the screws on the cell such that upon injection of the reaction mixture one of the modes of the cavity will be at 860 cm−1. In other words, during the reaction under VSC, the cavity mode will be on resonance with the νSi−C of PTA at normal incidence. For more details on cavity preparation, see Supplementary Materials and earlier reports , .
The VSC of the νSi−C as a function of the concentration of PTA (CPTA) was analysed first to determine a lower limit for VSC in the experiments. Different solutions of varying CPTA (5–0.5 m) were prepared by the addition of methanol to a neat PTA and were then injected to a tuned cavity described earlier. The VSC of the νSi−C of PTA in each experiment was evident from appearance of the split peaks (P+ and P−) in the IR transmission spectra shown in Figure 1B (red curves). The ħΩR of the strongly coupled νSi−C reduced from 128 cm−1 to 55 cm−1 with decreasing CPTA. The lowest CPTA used (0.5 m; 10% by volume in methanol) also satisfies the strong coupling conditions as its Rabi splitting energy (55 cm−1) is larger than the full width at half maximum of the cavity mode (28 cm−1) and the width of the νSi−C (39 cm−1). Further dilution of PTA with methanol resulted in reduced signal and was difficult to distinguish in the noisy background. The variation of ħΩR scaled linearly with square root of the CPTA, as shown in Figure 1C, confirming strong coupling for all the concentrations studied.
To analyse the impact of ħΩR on the thermodynamic parameters of the desilylation under VSC, the CPTA in the reaction mixture [a freshly prepared homogeneous solution of PTA and tetrabutyl ammonium fluoride (TBAF) in methanol] was varied from 0.5 to 2.5 m in steps of 0.5 m, while maintaining the fixed concentration of TABF (0.5 m). The presence of TBAF facilitates the desilylation reaction at room temperature . A higher CPTA resulted in immiscible solutions when mixed with TBAF and methanol and is therefore not suitable for experiments under VSC. The TBAF concentration was fixed based on the half-life of the reaction outside the cavity for the lowest CPTA and the homogeneity of the reaction mixtures for all the CPTA used.
In order to extract the thermodynamic parameters of the desilylation of PTA under VSC, the rate constants for the reaction with a fixed ħΩR (CPTA) in a small range of temperature (20°C–40°C) were determined as follows. The IR cell containing the tuned empty cavity is placed in a temperature-controlled jacket and is allowed to equilibrate to the set temperature. Subsequently, the reaction mixture is injected to the cavity without removing it from the temperature-controlled jacket. Precise tuning of the cavity mode to the νSi−C of PTA is ensured in every experiment and confirmed by the appearance of the P+ and P− in the IR transmission spectra. The rate constant of the silyl bond scission was determined following the relative shift of a higher-order cavity mode as explained in the earlier reports ,  and in the Supplementary Material (Figure S3). For comparison, reactions were also carried out outside the cavity (noncavity), as exactly described for experiments under VSC, using a cell made of SiO2-coated ZnSe. Plots showing the comparison of desilylation kinetics at 25°C, under VSC and noncavity, are given in the Supplementary Material (Figure S4).
The temperature-dependent rate constants of the desilylation for reaction mixtures with different ħΩR were then analysed based on the TS Eyring equation:
where k is the rate constant, kB is the Boltzmann constant, T is the temperature, h is Plank’s constant, and R is the universal gas constant. From the slope and the intercept of the Eyring plot, a plot of rate constant against inverse of temperature, enthalpy (ΔH‡), and entropy (ΔS‡) of activation of the TS is determined. The Eyring plots for various values of ħΩR and its comparison for a noncavity reaction for the same CPTA are shown in Figure 2. A clear change in slope of the Eyring plot is visible for reactions under VSC with respect to the variation in the ħΩR, whereas those for the noncavity experiments remained more or less the same for the different CPTA. The plots were then analysed quantitatively to determine the values of ΔH‡, ΔS‡, and the corresponding free energy of the TS (ΔG‡).
The red squares in Figure 3A and B show the variation of the ΔH‡ and TΔS‡ values under VSC as a function of ħΩR. As clearly seen in Figure 3A, ΔH‡ (35±4 kJ mol−1, lowest ħΩR) more than doubled to 84±8 kJ mol−1 for the reaction proceeding under the largest ħΩR, indicating a higher enthalpy requirement for the desilylation reaction with increasing coupling strength. The noncavity experiments did not exhibit noticeable variation in ΔH‡ as a function of CPTA (Figure 3A, blue circles), revealing that the change under VSC is not a concentration effect. The TΔS‡ (−49±6 kJ mol−1 for the lowest ħΩR, T=298 K) underwent a decrease in magnitude and nearly changed its sign to −1±7 kJ mol−1 for the largest ħΩR studied here, representing a gradual change from an organised to a less-ordered TS as the strong coupling strength is increased. Again, it can be seen that in the case of the noncavity conditions no such change in the TΔS‡ is observed (Figure 3B, blue circles). Figure 3C shows Δ(ΔG‡), the difference between ΔG‡ of the reaction under VSC and noncavity, as a function of ħΩR. The variation in Δ(ΔG‡) is relatively small compared to the large changes observed for ΔH‡ and ΔS‡, because the latter compensate each other. The upward trend of Δ(ΔG‡) with the increase of ħΩR is representative of the higher energy requirement for the desilylation under VSC compared to the noncavity and leads to the retarded reactivity up to a factor 4.5 as shown in Figure 3D and as reported earlier .
It is remarkable that although the changes in the thermodynamic parameters follow the Rabi splitting, their magnitude is much greater than the VSC induced ħΩR. This indicates that the effect of VSC is not only an energetic one on the potential energy surface leading to the products, it must also imply a change in reaction mechanism, which is confirmed by the evolution of ΔH‡ and ΔS‡ as discussed elsewhere . It is interesting to see that the VSC has a retarding effect on the reaction, which could be understood by looking into the formation and organisation of the TS of the rate-limiting step.
The thermodynamic parameter ΔH‡ mostly represents the reorganisation energy requirement for the formation of the TS, and ΔS‡ primarily indicates the degree of organisation of the TS . Under ambient conditions, the reaction studied here proceeds through a pentacoordinate intermediate and follows an associative mechanism with an ordered TS , , . The lower value of ΔH‡ (30±4 kJ mol−1) and the large negative value of TΔS‡ (−54±6 kJ mol−1, T=298 K) for the noncavity experiments agree well with the associative mechanism and with a TS that is more reactant/intermediate-like as schematically displayed in Figure 4A. Interestingly, the thermodynamic data show that the reaction proceeds through a more product-like TS as shown in Figure 4B under VSC. The dissociative TS requires higher energy as it is more bond-breaking, and it is less strained because of the tetrahedral-like structure. This could mean that the intermediate generated under VSC is easier to form but more stable than the noncavity case either by a charge stabilisation or by anisotropic interactions induced by light–matter coupling. For example, mode selective chemistry experiments have shown that TS control of the reaction can be modulated by anisotropic forces either to enhance or suppress a particular reaction path . In other words, the role of strong coupling could be more important in determining the position and structure of the TS as shown in recent experimental ,  and theoretical studies , .
In summary, we have shown that the thermodynamic parameters of a chemical reaction under VSC are indeed dependent on the strong coupling strength. This further illustrates that the polaritonic states are delocalised over all the molecules involved in the strong coupling and act collectively, even though the reactions are ultimately localised on single molecules. The observed changes in chemical dynamics of the present system, despite the relatively weak Rabi splitting energy compared to the TS energy of the reaction, point to structural changes in the TS under VSC. With more experiments and theoretical understanding about the TS changes under strong coupling, it should be possible to predict the catalytic activity and the product selectivity of reactions controlled by light–matter hybrid states.
Thomas A, George J, Shalabney A, et al. Ground-state chemical reactivity under vibrational coupling to the vacuum electromagnetic field. Angew Chem Int Ed 2016;55:11462–6. Web of ScienceCrossrefGoogle Scholar
Munkhbat B, Wersäll M, Baranov DG, Antosiewicz TJ, ShegaiT. Suppression of photo-oxidation of organic chromophores by strong coupling to plasmonic nanoantennas. Sci Adv 2018;4:eaas9552. Web of ScienceCrossrefPubMedGoogle Scholar
Lather J, Bhatt P, Thomas A, Ebbesen TW, George J. Cavity catalysis by cooperative vibrational strong coupling of reactant and solvent molecules. Angew Chem Int Ed 2019;58:10635–8. CrossrefWeb of ScienceGoogle Scholar
Stranius K, Hertzog M, Börjesson K. Selective manipulation of electronically excited states through strong light–matter interactions. Nat Commun 2018;9:2273. CrossrefWeb of SciencePubMedGoogle Scholar
Shalabney A, George J, Hutchison J, Pupillo G, Genet C, Ebbesen TW. Coherent coupling of molecular resonators with a microcavity mode. Nat Commun 2015;6:5981. Web of SciencePubMedCrossrefGoogle Scholar
Long JP, Simpkins BS. Coherent coupling between a molecular vibration and Fabry–Perot optical cavity to give hybridized states in the strong coupling limit. ACS Photonics 2015;2: 130–6. CrossrefWeb of ScienceGoogle Scholar
Muallem M, Palatnik A, Nessim GD, Tischler YR. Strong light–matter coupling and hybridization of molecular vibrations in a low-loss infrared microcavity. J Phys Chem Lett 2016;7: 2002–8. CrossrefWeb of ScienceGoogle Scholar
Campos-Gonzalez-Angulo J, Ribeiro RF, Yuen-Zhou J. Resonant catalysis of thermally-activated chemical reactions via vibrational polaritons. Nat Commun 2019;10:article no. 4685. Google Scholar
Flick J, Narang P. Excited-state nanophotonic and polaritonic chemistry with ab initio potential-energy surfaces. arXiv:1907.04646v1 [Physics, Physics:Quant-Ph] 2019. Google Scholar
Ribeiro RF, Martínez-Martínez LA, Du M, Campos-Gonzalez-Angulo J, Yuen-Zhou J. Polariton chemistry: controlling molecular dynamics with optical cavities. Chem Sci 2018;9:6325–39. PubMedCrossrefWeb of ScienceGoogle Scholar
Ruggenthaler M, Tancogne-Dejean N, Flick J, Appel H, Rubio A. From a quantum-electrodynamical light–matter description to novel spectroscopies. Nat Rev Chem 2018;2:0118. CrossrefWeb of ScienceGoogle Scholar
Dovzhenko DS, Ryabchuk SV, Rakovich YP, Nabiev IR. Light–matter interaction in the strong coupling regime: configurations, conditions, and applications. Nanoscale 2018;10:3589–605. CrossrefPubMedWeb of ScienceGoogle Scholar
Carey FA, Sundberg RJ. Advanced organic chemistry: Part A: structure and mechanisms. 5th ed. New York, NY, USA, Springer US, 2007. Google Scholar
Clayden J, Greeves N, Warren SG. Organic chemistry. New York, USA, Oxford University Press Inc.,2012. Google Scholar
Bassindale AR, Taylor PG. Reaction mechanisms of nucleophilic attack at silicon. In: Organic Silicon Compounds (1989). Wiley-Blackwell, 2004:839–92. Google Scholar
The online version of this article offers supplementary material (https://doi.org/10.1515/nanoph-2019-0340).
About the article
Published Online: 2019-11-29
Funding Source: International Center for Frontier Research
Award identifier / Grant number: ANR-10-EQPX-52-01
The authors acknowledge support of the International Center for Frontier Research in Chemistry (icFRC, Strasbourg), the ANR Equipex Union (ANR-10-EQPX-52-01, Funder Id: http://dx.doi.org/10.13039/501100001665), the Labex NIE projects (ANR-11-LABX-0058 NIE, Funder Id: http://dx.doi.org/10.13039/501100001665), and CSC (ANR-10-LABX- 0026 CSC, Funder Id: http://dx.doi.org/10.13039/501100001665) within the Investissement d’Avenir program ANR-10-IDEX-0002-02, Funder Id: http://dx.doi.org/10.13039/501100001665, the ERC (project no. 788482 MOLUSC and no. 639170 CARBONFIX), and QuantERA projet RouTe.