Arxiv:2103.06749V1 [Physics.Chem-Ph] 11 Mar 2021
Total Page:16
File Type:pdf, Size:1020Kb
Collective vibrational strong coupling effects on molecular vibrational relaxation and energy transfer: Numerical insights via cavity molecular dynamics simulations Tao E. Li,1, ∗ Abraham Nitzan,1, 2, y and Joseph E. Subotnik1, z 1Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA 2School of Chemistry, Tel Aviv University, Tel Aviv 69978, Israel For a small fraction of hot CO2 molecules immersed in a liquid-phase CO2 thermal bath, classical cavity molecular dynamics simulations show that forming collective vibrational strong coupling (VSC) between the C O asymmetric stretch of CO2 molecules and a cavity mode accelerates hot- molecule relaxation. The physical mechanism underlying this acceleration is the fact that polaritons, especially the lower polariton, can be transiently excited during the nonequilibrium process, which facilitates intermolecular vibrational energy transfer. The VSC effects on these rates (i) resonantly depend on the cavity mode detuning, (ii) cooperatively depend on molecular concentration or Rabi splitting, and (iii) collectively scale with the number of hot molecules, which is similar to Dicke’s superradiance. For larger cavity volumes, due to a balance between this superradiant-like behavior and a smaller light-matter coupling, the total VSC effect on relaxation rates can scale slower than 1=N, and the average VSC effect per molecule can remain meaningful for up to N ∼ 104 molecules forming VSC. Moreover, we find that the transiently excited lower polariton prefers to relax by transferring its energy to the tail of the molecular energy distribution rather than equally distributing it to all thermal molecules. Finally, we highlight the similarities of parameter dependence between the current finding with VSC catalysis observed in Fabry–Pérot microcavities. I. INTRODUCTION In order to narrow the gap between theory and ex- periment, here we numerically investigate VSC effects on Collective vibrational strong coupling (VSC) can oc- two nonequilibrium processes — molecular vibrational cur if a macroscopic number of liquid-phase molecules energy relaxation and intermolecular vibrational energy are confined to a Fabry–Pérot microcavity and a molecu- transfer. These vibrational processes has been exten- lar vibrational mode is near resonant with a cavity mode sively studied both experimentally and theoretically out- [1, 2]. Under collective VSC, experimental reports indi- side a cavity, and the rates of which have been known cate not only a peak splitting, i.e., a Rabi splitting within to play an important role in many physical and chemi- molecular infrared (IR) spectroscopy, but also the mod- cal processes, including chemical reactions [15]. Inside a ification of chemical reaction rates [3–6] and crystalliza- cavity, a recent experiment [16] has studied the effect of tion processes [7] under thermal conditions. As pioneered VSC on intermolecular vibrational energy transfer rates first by Ebbesen [3] and co-workers, these observations by quantifying the response of hybrid light-matter states suggest that collective VSC might meaningfully modify (polaritons) after pumping the upper polariton (UP) for 12 13 individual molecular properties without external pump- a liquid mixture of W( CO)6 and W( CO)6. ing — although these intriguing experimental findings For the sake of simplicity, here our numerical study fo- cannot yet be well explained by current theory [8–13]. cuses on a pure liquid CO2 system when the C O asym- A simple example illustrating how conventional the- metric stretch forms VSC with a single optical cavity ory fails to explain the Ebbesen experiments is to con- mode (where two polarization directions are included). sider the case of N molecules forming VSC with a Rabi In such a system, instead of exciting polaritons, we will −1 splitting ΩN = 2g0pN 100 cm , where g0 denotes consider the case when a small fraction of uncorrelated the light-matter coupling∼ for individual molecules. Be- hot CO2 molecules dissipates and transfers vibrational cause g0 (= ΩN =2pN) is negligible when N becomes energy to the remaining thermal CO2 molecules at room macroscopic, one would guess that individual molecular temperature. Unlike many experiments and theoreti- properties (such as chemical reaction rates) cannot be cal studies concentrating on the polaritonic response, we meaningfully modified by a Fabry–Pérot microcavity, a will mainly focus on how individual molecules (which are theoretical prediction at odds with several experiments. arXiv:2103.06749v1 [physics.chem-ph] 11 Mar 2021 mostly composed of vibrational dark modes) relax and Recent efforts [10, 14] also suggest that, within a clas- transfer energy under VSC. In detail, we will extensively sical description of cavity photons and molecules, static study how vibrational energy relaxation and transfer de- properties of individual molecules during thermal equi- pend on cavity mode detuning, molecular concentration librium are entirely unchanged under usual VSC setups, (or Rabi splitting), and the number of hot molecules. Be- indicating a nonequilibrium (or perhaps quantum) origin cause there is no external polariton pumping, our investi- of the Ebbesen experiments. gation of how a cavity affects vibrational relaxation and energy transfer will hopefully yield insight into the VSC modifications of individual molecular properties (such ∗ [email protected] as chemical reaction rates) that are observed in experi- y [email protected] ments. In particular, by quantifying the asymptotic scal- z [email protected] ing of VSC effects with molecular system size or effective 2 ps). Below, we will report how VSC affects vibrational energy relaxation and transfer using CavMD simulations. II. RESULTS A. VSC effects on vibrational energy relaxation and transfer Fig. 2a plots the IR spectrum outside the cavity (black line; " = 0) or inside the cavity (red line; " = 2 10−4 −×1 a.u.) when the cavity mode (at !c = 2320 cm ; the FIG. 1. Sketch of the simulation setup where a large num- dashede vertical blue line) forms VSC withe the C O ber of liquid-phase carbon dioxide molecules forms VSC with asymmetric mode (peaked at ! = 2327 cm−1) of liquid a single cavity mode. The left cartoon shows that the vi- 0 CO2. Inside the cavity, a pair of lower (LP; peaked at brational energy relaxation and transfer from hot (upper) to −1 −1 thermal (bottom) molecules inside a cavity can be accelerated 2241 cm ) and upper (UP; peaked at 2428 cm ) polari- relative to that outside a cavity due to polariton-accelerated tons form and these polaritons are separated by a Rabi −1 intermolecular vibrational energy transfer. splitting of 187 cm . The IR spectrum is calculated by evaluating the Fourier transform of the dipole autocor- relation function from equilibrium trajectories; see Ap- pendix S-III for details. cavity volumes, our study will also partly address if VSC We now consider a nonequilibrium process where effects can persist and affect the properties of individ- Nhot = 10 uncorrelated hot molecules are immersed in a ual molecules in the limit that a very large number of thermal CO bath at room temperature (where in total molecules are present in a cavity. 2 there are Nsub = 216 molecules in the simulation cell); see The theoretical approach we will take is classical cav- Appendix S-II for details. Fig. 2b plots the average time- ity molecular dynamics (CavMD) simulations [14, 17], a resolved C O bond potential energy per hot molecule newly developed numerical tool implemented by the au- outside (black line) or inside (red line) the cavity, where thors to classically propagate the coupled dynamics be- a thermal energy kBT = 300 K has been subtracted from tween realistic molecules (assumed to stay in their elec- the C O bond potential energy; note that here we use tronic ground-state) and cavity photons in the dipole kBT instead of kBT=2 since each CO2 molecule contains gauge. Since the self-dipole term is included in the light- two C O bonds. As shown in Fig. 2b, the initial po- matter Hamiltonian of CavMD simulations, this numeri- tential energy in the two C O bonds per hot molecule cal approach preserves gauge invariance and maintains 3 is roughly 2~!c ( 6 10 K), i.e., the initial tempera- numerical stability [18]. Compared with VSC experi- ture of the hot molecules≈ × is 3 103 K. At later times, ments, this approach has reliably captured many VSC- the vibrational energy relaxation∼ × of the hot molecules in- induced phenomena, including an asymmetric Rabi split- side the cavity is accelerated compared with that outside ting [14, 19], polariton relaxation to vibrational dark the cavity. Meanwhile, as shown in Fig. 2c, the average modes on a time scale of ps and sub-ps [17, 20], and C O bond potential energy per thermal molecule inside a delay of population gain in the singly excited mani- (red line) the cavity increases faster than that outside fold of vibrational dark modes after pumping the lower (black line) the cavity. Here, "thermal molecules" refer polariton (LP), a process which stems from polariton en- to molecules that were prepared at thermal equilibrium. hanced molecular nonlinear absorption [17, 21]. Hence, During this nonequilibrium process, the total system en- CavMD simulations appears to be a promising tool to ergy is conserved: the simulation is performed under a study VSC-related phenomena. NVE (constant number, volume, and energy) ensemble; A brief introduction of CavMD is given in Appendix see simulation details in Appendix S-I. S-I and S-II; see Ref. [17] for more details regarding During the energy relaxation and transfer process, in- CavMD simulations of a liquid CO2 system and how the side the cavity, Fig. 2d plots the total (kinetic + po- −1 CO2 force field is defined. In short, as shown in Fig. 1, tential) energy of the cavity photons (!c = 2320 cm CavMD simulates a system with Nsub CO2 molecules in and with two polarization directions) subtracted by the a periodic cell coupled to a single cavity mode (with two thermal background 2kBT .