Probing voltage-induced chemical reactions and anharmonicity with a confined vacuum light field

In this work, we present a proof-of-concept investigation of non-equilibrium chemical reaction dynamics at a molecule-electrode interface, driven out of equilibrium by an applied votage bias and mediated by a confined, enhanced vacuum electromagnetic field inside an optical cavity. The coupled electron-vibration-photon system, together with the electrodes and a dissipative environment, is described within an open quantum system framework and solved using a numerically exact quantum dynamical approach. The reaction coordinate is modeled with a Morse potential, enabling explicit treatment of molecular anharmonicity and bond-breaking behavior. By varying the cavity frequency across the infrared regime to cover typical nuclear vibrational energies, we observe multiple resonant rate suppression features that emerge whenever the cavity mode is brought into resonance with a dipole-allowed vibrational transition along the anharmonic ladder up to the dissociation threshold. These findings open the door to extending polaritonic chemistry into genuinely nonequilibrium scenarios relevant to molecule-electrode interfaces. Moreover, building on these results, we further propose a multi-mode vibrational strong coupling strategy in which several cavity modes are individually matched to distinct vibrational transitions. This engineered multi-resonant cavity induces a stepwise vibrational ladder descending process that efficiently drains vibrational excited energy. The resulting cavity-assisted cooling suggests a potential route toward mitigating voltage-induced bond rupture and the long-standing instability issues of molecular junctions operating under high bias.

💡 Research Summary

In this paper the authors present a theoretical proof‑of‑concept study of how a confined, enhanced vacuum electromagnetic field inside an optical cavity can influence voltage‑driven chemical reactions at a single‑molecule junction. The system consists of a molecule bridging two electrodes under a bias voltage, a single reactive vibrational mode that serves as the reaction coordinate, and one or several quantized cavity modes. The reaction coordinate is modeled by a Morse potential for both the neutral and the charged electronic states, thereby incorporating anharmonicity and a realistic dissociation barrier. Electron transport between the molecule and the electrodes is described by fermionic reservoirs with chemical potentials ±Φ/2, while vibrational relaxation and cavity loss are modeled by bosonic baths.

All degrees of freedom—electronic, vibrational, photonic, and the surrounding environments—are treated on an equal quantum footing within an open‑quantum‑system framework. The authors employ the hierarchical equations of motion (HEOM) method, combined with a compact tensor‑network state representation, to obtain numerically exact dynamics. The fermionic and bosonic bath correlation functions are decomposed into sums of exponentials using Lorentzian spectral densities and Padé pole expansions, which allows the construction of an extended wavefunction that includes auxiliary pseudomodes for both electrons and bosons.

Key findings emerge when the cavity frequency is swept across the infrared region. Whenever a cavity mode becomes resonant with a dipole‑allowed vibrational transition—whether a fundamental (Δv = 1) or an overtone (Δv = 2, 3)—the reaction rate for bond rupture is sharply suppressed. This “resonant suppression” originates from the cavity efficiently extracting vibrational energy and converting it into photons, thereby reducing the population of highly excited vibrational states that would otherwise cross the dissociation barrier. The effect is highly selective: off‑resonant cavity frequencies produce negligible changes in the rate.

Building on this single‑mode insight, the authors propose a multi‑mode strong‑coupling strategy. By designing a cavity that hosts several modes, each tuned to a distinct vibrational transition along the anharmonic ladder, a stepwise vibrational‑ladder‑descending pathway is opened. Excited vibrational quanta cascade down the ladder, emitting photons at each step, which results in an effective cavity‑assisted cooling of the molecular vibration. This cooling dramatically lowers the probability of voltage‑induced bond rupture and suggests a practical route to stabilizing molecular junctions operating under high bias.

The study highlights the crucial role of anharmonicity: the Morse potential yields a dense, non‑uniform set of vibrational energy gaps, enabling multiple resonances that would be absent in a purely harmonic model. Consequently, multi‑mode designs become feasible without requiring an impractically large number of cavity modes.

Finally, the authors discuss experimental prospects. They suggest that nanoparticle‑on‑mirror or metasurface‑based plasmonic cavities could provide the required field confinement and tunability, while graphene or thin‑film electrodes could serve as the leads. They also outline future directions, including real‑time spectroscopic monitoring of cavity‑mediated reactions, extension to multiple electronic orbitals and vibrational modes, and algorithmic improvements to scale the HEOM‑tensor‑network approach to larger systems.

Overall, the paper extends polaritonic chemistry from equilibrium, thermally activated reactions to genuinely nonequilibrium, voltage‑driven processes, offering a clear mechanistic picture of how vacuum light fields can be harnessed to control bond breaking at the nanoscale.