Cavity molecular dynamics simulations of vibrational polariton-enhanced molecular nonlinear absorption

• Published in: The Journal of chemical physics 2021-03, Vol.154 (9), p.094124-094124
• Main Authors: Li, Tao E., Nitzan, Abraham, Subotnik, Joseph E.
• Format: Article
• Language: English
• Subjects: Carbon dioxide; Coupling (molecular); Excitation; Holes; Manifolds; Molecular dynamics; Multiphoton absorption; Nonlinear response; Numerical prediction; Optical properties; Polaritons; Subsystems
• ISSN: 0021-9606; 1089-7690
• DOI: 10.1063/5.0037623

Recent experiments have observed that the chemical and photophysical properties of molecules can be modified inside an optical Fabry–Pérot microcavity under collective vibrational strong coupling (VSC) conditions, and such modification is currently not well understood by theory. In an effort to understand the origin of such cavity-induced phenomena, some recent studies have focused on the effect of the cavity environment on the nonlinear optical response of the molecular subsystem. Here, we use a recently proposed protocol for classical cavity molecular dynamics simulations to numerically investigate the linear and the nonlinear response of liquid carbon dioxide under such VSC conditions following an optical pulse excitation. We find that applying a strong pulse of excitation to the lower hybrid light–matter state, i.e., the lower polariton (LP), can lead to an overall molecular nonlinear absorption that is enhanced by up to two orders of magnitude relative to the excitation outside the cavity. This polariton-enhanced multiphoton absorption also causes an ultrashort LP lifetime (0.2 ps) under strong illumination. Unlike usual polariton relaxation processes—whereby polaritonic energy transfers directly to the manifold of singly excited vibrational dark states—under the present mechanism, the LP transfers energy directly to the manifold of higher vibrationally excited dark states; these highly excited dark states subsequently relax to the manifold of singly excited states with a lifetime of tens of ps. Because the present mechanism is generic in nature, we expect these numerical predictions to be experimentally observed in different molecular systems and in cavities with different volumes.