Nuclear-Electronic Quantum Dynamics in a Plasmonic Nanocavity

This paper demonstrates that real-time nuclear-electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) coupled to multimode, lossy plasmonic nanocavities can efficiently simulate and reveal how strong light-matter coupling modifies ultrafast excited-state proton transfer dynamics and generates distinct spectroscopic signatures, such as Rabi-like oscillations and resonance evolution.

The Gist

Imagine you have a tiny, super-fast chemical reaction happening inside a molecule—specifically, a proton (a hydrogen nucleus) jumping from one side of the molecule to the other. This happens in the blink of an eye, faster than a camera shutter could ever capture. Now, imagine you want to watch this jump happen, but you also want to see if you can influence the jump by putting the molecule inside a special "box" made of metal that traps light.

This paper is about building a virtual microscope to watch these ultra-fast jumps and seeing how the "light box" (called a plasmonic nanocavity) changes the game.

Here is the breakdown of their work using simple analogies:

1. The Stage: The "Gold Mirror" Box

Think of a plasmonic nanocavity like a tiny, high-tech echo chamber made of gold.

• Normal mirrors reflect light like a billiard ball hitting a wall.
• This gold box is different. It's so small that it traps light in a space smaller than the light wave itself. It creates a "storm" of electric fields.
• The Problem: This box is messy. It's not just one clear echo; it's a chaotic mix of many different frequencies (like a drum kit being hit all at once), and the light leaks out very quickly (it's "lossy"). Modeling this on a computer is usually a nightmare because there are too many variables.

2. The Actors: The "Dancing Proton"

The scientists are watching a specific molecule called oHBA (and later a similar one called AMIEP).

• Inside this molecule, a proton is like a dancer waiting to jump from one partner (the donor) to another (the acceptor).
• When you shine light on the molecule, the dancer gets excited and starts moving. In a normal world, this jump happens in about 10 to 100 femtoseconds (that's a quadrillionth of a second).
• The scientists use a special computer method called RT-NEO-TDDFT. Think of this as a super-slow-motion camera that treats the proton not as a tiny ball, but as a fuzzy cloud of probability. This allows them to see exactly how the "cloud" moves and jumps.

3. The Experiment: Watching the Light "Sing"

The researchers put their dancing proton inside their virtual gold box. They asked two main questions:

Scenario A: The "Passive Observer" (Weak Coupling)

Imagine the gold box is just a quiet audience. The molecule jumps, and the box listens.

• What happened: As the proton jumps, the molecule changes its shape and energy.
• The Result: The gold box "sings" back. Because the box has many different frequencies, different parts of the box start vibrating at different times.
• The Analogy: It's like a rainbow of sound. At first, the molecule is high-energy, so the box sings a high note. As the proton jumps and settles, the molecule's energy drops, and the box's song shifts to a lower note.
• Why it matters: By listening to which note the box is singing at what time, the scientists can reconstruct the movie of the proton's jump without ever looking at the molecule directly. It's like deducing a dancer's moves by listening to the changing pitch of the music they are dancing to.

Scenario B: The "Active Partner" (Strong Coupling)

Now, imagine they turn up the volume in the gold box so much that the light and the molecule become best friends. They stop being two separate things and become a hybrid creature called a polariton.

• What happened: The light is so strong that it grabs the proton and says, "Wait, don't jump yet!"
• The Result: The proton transfer gets slowed down or even stopped. The system starts doing a "Rabi oscillation"—think of it like a swing. The energy swings back and forth between the molecule and the light, over and over, before the proton finally moves.
• The Analogy: It's like a tug-of-war. The light pulls the proton one way, the molecule pulls it the other. They get stuck in a loop, unable to finish the jump they were supposed to make.

4. The Real-World Test: The "Gold Nanoparticle"

To make sure this wasn't just a fantasy, they tested it with a setup that actually exists in real labs: a Gold Nanoparticle on a Mirror (NPoM).

• This is a tiny gold ball sitting on a flat gold mirror.
• They used a different molecule (AMIEP) that jumps at a lower energy.
• The Surprise: Even though the gold box wasn't perfectly tuned to the molecule at the start, as the molecule relaxed and changed its energy, it naturally "slid" into resonance with the box. The box started singing loudly only after the molecule had moved.
• The Conclusion: They found that for a single molecule, the effect is subtle. But if you have a small group of molecules (like 4 or 9) in the box, the effect becomes huge. The light and matter lock together, creating clear "polariton" peaks.

The Big Picture

This paper is a breakthrough because it bridges the gap between theory and reality.

• Before: Scientists could only model simple, perfect boxes or tiny molecules.
• Now: They have a tool that can simulate messy, real-world gold boxes with many frequencies and watch how they control chemical reactions in real-time.

In short: They built a digital time machine that lets us watch a proton jump in slow motion, and they discovered that by tuning the "light box" around it, we can either listen to the jump to understand chemistry better, or stop the jump entirely to control chemical reactions. This could lead to new ways to design solar cells, sensors, or even control how drugs work inside the body.

Technical Summary

1. Problem Statement

Plasmonic nanocavities offer a unique platform for achieving strong light–matter coupling at the single-molecule level and enhancing spectroscopies. However, modeling these systems is theoretically challenging due to two primary factors:

• Multimodal Character: Unlike micron-scale Fabry-Perot cavities, plasmonic nanocavities possess strongly multimodal electromagnetic fields with complex spectral densities.
• High Loss Rates: These nanocavities have extremely short lifetimes (high loss), making them significantly more dissipative than traditional optical cavities.

Existing theoretical methods face a trade-off:

• Quantum Electrodynamical (QED) electronic structure theory captures electron-photon interactions accurately but is limited to small systems, simple lossless environments, and often neglects non-adiabatic nuclear-electronic effects.
• Non-adiabatic dynamical methods (e.g., surface hopping) often struggle to accurately describe complex, multimode nanocavity environments coupled with nuclear quantum effects.

There is a critical need for a framework that can simultaneously treat nuclear quantum effects (specifically proton transfer), electronic dynamics, and complex, lossy, multimode electromagnetic environments to understand how these cavities influence chemical reactivity and serve as probes.

2. Methodology

The authors employ Real-Time Nuclear–Electronic Orbital Time-Dependent Density Functional Theory (RT-NEO-TDDFT) coupled with classical cavity modes.

• RT-NEO-TDDFT Framework:

  • Treats electrons and specific nuclei (typically protons involved in transfer) quantum mechanically on an equal footing.
  • Propagates the quantum mechanical densities of electrons and protons in real time without invoking the Born-Oppenheimer approximation between them.
  • Uses the epc17-2 electron-proton correlation functional and PB4D protonic basis sets within the Q-Chem software package.

• Cavity Modeling:

  • Classical Modes: Cavity modes are treated as classical harmonic oscillators coupled to the molecular dipole. This approximation allows for the inclusion of many modes without a prohibitive increase in computational cost.
  • Loss and Spectral Density: Cavity loss is incorporated via a mode-specific inverse lifetime ($\gamma_\alpha$). The environment is defined by a spectral density $J(\omega)$, which describes the lineshape of modes due to finite cavity lifetimes (modeled as Lorentzians).
  • Two Cavity Models:
    1. Gaussian Multimode Environment: A tunable, discretized Gaussian distribution of modes used to study general principles of coupling strength and resonance.
    2. Nanoparticle-on-Mirror (NPoM): A physically realistic model based on transformation optics for a gold nanoparticle (30 nm radius) separated from a mirror by a 0.9 nm gap. This generates a specific, experimentally relevant multimodal spectral density.

• Observables:

  • Proton Dynamics: Tracked via the expectation value of proton position ($H-O_D$ and $H-O_A$ distances).
  • Cavity Emission: Calculated as the mode-specific emission rate $I_\alpha(t)$, proportional to the mode occupation number and loss rate. This serves as a time- and energy-resolved probe of the molecular dynamics.

3. Key Contributions

1. Development of a Hybrid Framework: The paper establishes a robust RT-NEO-TDDFT framework coupled to multiple classical cavity modes that includes cavity loss. This enables the simulation of chemical reactions in physically realistic, lossy, multimode electromagnetic environments.
2. Cavity as a Spectroscopic Probe: The authors demonstrate that the time- and energy-resolved cavity emission from a multimode cavity can act as a passive probe for ultrafast nuclear-electronic dynamics. The emission spectrum shifts as the molecular system evolves, reflecting changes in the instantaneous transition energy.
3. Mechanism of Strong Coupling in Lossy Systems: The study elucidates how strong coupling affects proton transfer in plasmonic cavities, showing that polariton formation can suppress chemical reactions (proton transfer) and induce Rabi-like oscillations in the cavity emission.
4. Resonance Evolution: The work reveals that a molecule undergoing excited-state proton transfer (ESIPT) can dynamically evolve into resonance with a cavity, even if it starts out of resonance with the dominant cavity peak.

4. Results

A. Intermediate Coupling Regime (Passive Probing)

• System: o-hydroxybenzaldehyde (oHBA) in a Gaussian multimode cavity.
• Observation: When the light-matter coupling is weak ($g_0/\mu_0 = 0.03$ eV·nm$^{-1}$e$^{-1}$), the cavity does not alter the proton transfer dynamics.
• Probing Capability: The cavity emission spectrum acts as a real-time monitor.
  • In the S0 geometry (no proton transfer), emission remains centered near the initial vertical excitation energy ($\Delta E_{S1}$).
  • In the constrained S1 geometry (ultrafast proton transfer), the emission maximum shifts rapidly to lower energies (red-shift) as the proton moves toward the acceptor. This shift corresponds to the lowering of the instantaneous S0-S1 transition energy.
  • The emission intensity increases dramatically after proton transfer due to the localization of the proton density, which narrows the distribution of excitation energies and enhances resonance with lower-energy cavity modes.

B. Strong Coupling Regime (Active Modification)

• System: oHBA with increased coupling strength ($g_0/\mu_0 = 0.3$ eV·nm$^{-1}$e$^{-1}$).
• Observation: The cavity actively modifies the dynamics.
  • Suppression: In the constrained S1 geometry, the cavity suppresses the proton transfer. The proton is inhibited from fully transferring and recrosses back to the donor side.
  • Rabi-like Oscillations: The cavity emission exhibits complex, multi-timescale behavior, including rapid modulations and longer-timescale Rabi-like oscillations in the frequency of maximum emission. This indicates the formation of polaritons (hybrid light-matter states).

C. Realistic NPoM Cavity Simulation

• System: Schiff base molecule (AMIEP) in a gold NPoM cavity.
• Initial Mismatch: The initial vertical excitation of AMIEP (2.61 eV) is not resonant with the dominant NPoM peak (2.41 eV).
• Dynamic Resonance: As the molecule undergoes ESIPT, its excitation energy lowers. By ~15 fs, the molecular subsystem evolves into resonance with the dominant cavity mode, leading to significant cavity emission.
• Polariton Formation:
  • For a single molecule, even at resonance, well-defined upper and lower polariton peaks were not clearly resolved due to the broad linewidth of the NPoM spectral density.
  • By scaling the coupling by $\sqrt{N_{mol}}$ (simulating 4 and 9 molecules), clear Rabi splitting (30 meV for $N=4$, 50 meV for $N=9$) was observed, demonstrating that strong coupling is achievable with small ensembles in realistic nanocavities.

5. Significance

• Bridging Theory and Experiment: The study provides a theoretical bridge between complex experimental setups (plasmonic nanocavities) and quantum chemical dynamics, offering a way to interpret experimental spectroscopic signatures.
• Control of Reactivity: It demonstrates that plasmonic nanocavities can not only probe but also actively control chemical reaction pathways (e.g., suppressing proton transfer) through strong coupling, even in highly lossy environments.
• New Spectroscopic Paradigm: The work proposes a new method for single-molecule sensing where the cavity emission itself provides time- and energy-resolved information about the evolving quantum state of a molecule, bypassing the need for traditional fluorescence detection which may be too slow for ultrafast processes.
• Scalability: The RT-NEO framework with classical cavity modes is shown to be computationally efficient enough to handle the vast number of degrees of freedom required for realistic multimode simulations, paving the way for studying more complex chemical systems in engineered electromagnetic environments.