Watching Molecules Dance with Light

Author: Denis Avetisyan


New simulations reveal how light trapped in nanoscale cavities can influence and track the ultrafast movements of atoms within molecules.

The interplay between a gold nanoparticle and a mirror surface creates a constrained environment where emission frequencies <span class="katex-eq" data-katex-display="false">\omega_k</span> shift over time, evidenced by a normalized emission intensity <span class="katex-eq" data-katex-display="false">I_k</span> dependent on the dielectric functions of gold <span class="katex-eq" data-katex-display="false">\epsilon_{Au}(\omega)</span> and the surrounding solvent <span class="katex-eq" data-katex-display="false">\epsilon_{env}</span>, as well as the distance δ between the nanoparticle and mirror-a phenomenon further characterized by a spectral density <span class="katex-eq" data-katex-display="false">J(\omega)</span> comprised of summed Lorentzian modes and exhibiting negligible difference in proton position expectation values whether or not a cavity is present.
The interplay between a gold nanoparticle and a mirror surface creates a constrained environment where emission frequencies \omega_k shift over time, evidenced by a normalized emission intensity I_k dependent on the dielectric functions of gold \epsilon_{Au}(\omega) and the surrounding solvent \epsilon_{env}, as well as the distance δ between the nanoparticle and mirror-a phenomenon further characterized by a spectral density J(\omega) comprised of summed Lorentzian modes and exhibiting negligible difference in proton position expectation values whether or not a cavity is present.

Real-time modeling of nuclear-electronic dynamics in plasmonic nanocavities demonstrates control over excited-state proton transfer.

Confining light and matter at the nanoscale presents a significant challenge to accurately modeling quantum dynamics. This is addressed in ‘Nuclear-Electronic Quantum Dynamics in a Plasmonic Nanocavity’, which details a theoretical framework combining real-time nuclear-electronic orbital time-dependent density functional theory with multimode plasmonic nanocavities to simulate strong light-matter coupling. Our simulations reveal that cavity emission can not only probe, but also potentially control, ultrafast excited-state proton transfer reactions, even when the system is initially offresonant. Could this approach pave the way for designing nanoscale environments that actively steer chemical reactions with unprecedented precision?


The Inevitable Interplay: Bridging Quantum and Classical Realms

Simulating the dynamic behavior of molecules demands a simultaneous consideration of both the electrons and the nuclei that comprise them, presenting a significant computational challenge. Molecular motion isn’t simply the nuclei traversing a static electronic landscape; rather, the electrons and nuclei constantly influence each other. Accurately representing this interplay requires solving the Schrödinger equation for all particles, a task that scales exponentially with the number of constituents. Consequently, even relatively small molecules can overwhelm conventional computing resources, necessitating approximations or the utilization of high-performance computing infrastructure. The complexity arises from the differing timescales of electronic and nuclear motion – electrons move much faster – demanding computational methods capable of resolving these disparate rates to faithfully capture molecular phenomena like chemical reactions and energy transfer.

Molecular simulations conventionally employ the Born-Oppenheimer approximation, a simplification that treats the nuclei as stationary while the electrons adjust instantaneously. This approach dramatically reduces computational demands, yet it falters when electronic and nuclear motions become intimately linked-a scenario common in photochemistry, chemical reactions, and certain materials. When nuclei move rapidly enough to significantly perturb the electronic structure, or when electronic states become closely spaced, the approximation breaks down, leading to inaccurate predictions of molecular behavior. This failure arises because the assumption of separable nuclear and electronic motion is no longer valid, necessitating more sophisticated methods capable of explicitly accounting for the interplay between these degrees of freedom to achieve a truly accurate depiction of dynamic molecular processes.

Accurately simulating molecular behavior often demands computational frameworks that move beyond approximations like the Born-Oppenheimer method, particularly when electronic and nuclear motions become intertwined – a phenomenon known as non-adiabaticity. These advanced frameworks strive to explicitly account for the coupling between electronic states, allowing for the observation of effects such as energy transfer between states and the possibility of transitions that are forbidden under the Born-Oppenheimer approximation. Such capabilities are crucial for modeling complex chemical reactions, photochemistry, and the behavior of molecules in strong fields, ultimately providing a more complete and realistic depiction of molecular dynamics. These methods frequently involve solving the time-dependent Schrödinger equation on multiple potential energy surfaces simultaneously, a computationally intensive task but essential for understanding processes where electronic and nuclear degrees of freedom are strongly correlated.

The electronic power spectrum of AMIEP coupled to a nanoscale photonic molecule (NPoM) exhibits a modified spectral density <span class="katex-eq" data-katex-display="false">J(\omega)</span> due to light-matter coupling, scaling with the number of molecules <span class="katex-eq" data-katex-display="false">N_{mol}</span> to reveal enhanced excitation energies <span class="katex-eq" data-katex-display="false">\Delta E_{S1}</span> originating from the initial S0 to S1 transition.
The electronic power spectrum of AMIEP coupled to a nanoscale photonic molecule (NPoM) exhibits a modified spectral density J(\omega) due to light-matter coupling, scaling with the number of molecules N_{mol} to reveal enhanced excitation energies \Delta E_{S1} originating from the initial S0 to S1 transition.

RT-NEO-TDDFT: A Framework for Coherent Dynamics

RT-NEO-TDDFT (Real-Time Non-Equilibrium Orbital-dependent Time-Dependent Density Functional Theory) is a quantum chemistry methodology designed to simulate the simultaneous evolution of electronic and nuclear degrees of freedom. Unlike methods relying on the Born-Oppenheimer approximation, RT-NEO-TDDFT treats electrons and nuclei on equal footing, allowing for the investigation of non-adiabatic processes where the electronic and nuclear motions are strongly coupled. This is achieved through the direct propagation of both electronic and nuclear wavepackets in time, offering a robust framework for studying phenomena such as ultrafast photochemistry, charge transfer, and energy dissipation in complex molecular systems. The method’s robustness stems from its ability to accurately describe both ground and excited state dynamics without relying on perturbative treatments of non-adiabatic coupling.

RT-NEO-TDDFT utilizes the Von Neumann equations, a set of equations governing the time evolution of quantum states, to simultaneously propagate the electronic and nuclear wavefunctions. This is achieved through a fully coupled treatment, where the time derivative of the density matrix, \frac{d\rho}{dt} = -\frac{i}{\hbar}[H, \rho], is calculated for both electronic and nuclear degrees of freedom. Unlike approaches that treat nuclei classically, this method allows for a direct calculation of the time-dependent wavefunction for all particles, enabling the simulation of non-adiabatic effects and avoiding the need for separate electronic and nuclear propagation schemes. The resulting equations are discretized and solved self-consistently at each time step to obtain the updated quantum state of the entire system.

The integration of classical cavity modes within the RT-NEO-TDDFT framework utilizes the Velocity Verlet Algorithm for time propagation. This algorithm efficiently solves the classical equations of motion for the electromagnetic field, treating the cavity modes as harmonic oscillators with time-dependent forces induced by the molecular system. The Velocity Verlet Algorithm is particularly suited for this purpose due to its symplectic properties, which preserve energy conservation over extended simulation times – crucial for accurately modeling light-matter interactions. Specifically, the algorithm iteratively updates the positions and velocities of the cavity mode coordinates based on the calculated forces, enabling the simulation of phenomena such as strong coupling regimes and the exchange of energy between the molecule and the electromagnetic field. The time step used in the Velocity Verlet integration is determined by the frequencies of the cavity modes to ensure stability and accuracy.

The Born-Oppenheimer approximation, a cornerstone of computational chemistry, assumes that the motion of nuclei and electrons can be treated separately due to the significant mass difference between them. However, this approximation breaks down when potential energy surfaces become closely spaced or intersect, leading to non-adiabatic effects. RT-NEO-TDDFT circumvents this limitation by directly incorporating nuclear motion into the quantum mechanical propagation, effectively treating electrons and nuclei as a single, coupled quantum system. This is achieved through the use of the Von Neumann equations, which allow for the simultaneous evolution of both electronic and nuclear wavefunctions, thus explicitly accounting for non-adiabatic coupling and enabling the study of phenomena where the Born-Oppenheimer approximation is invalid, such as photochemistry and electron transfer.

Cavity interactions modulate proton transfer dynamics in oHBA, shifting the frequency of maximal emission <span class="katex-eq" data-katex-display="false">\omega^{\text{max}}_{k}</span> and altering H-OD and H-OA distances (red indicates cavity on, blue indicates cavity off) compared to the uncoupled system.
Cavity interactions modulate proton transfer dynamics in oHBA, shifting the frequency of maximal emission \omega^{\text{max}}_{k} and altering H-OD and H-OA distances (red indicates cavity on, blue indicates cavity off) compared to the uncoupled system.

Enhanced Proton Transfer Through Nanophotonic Confinement

Excited-State Proton Transfer (ESIPT) is a fundamental photophysical process occurring in a diverse range of biological and chemical systems, including fluorescent proteins, light-harvesting complexes, and certain organic molecules. This process involves the transfer of a proton from a donor molecule to an acceptor molecule following photoexcitation, typically occurring on a femtosecond timescale. ESIPT is critical for energy dissipation, fluorescence modulation, and the initiation of photochemical reactions. Its efficiency is highly sensitive to the surrounding environment, including solvent polarity, hydrogen bonding networks, and the specific molecular structure of the donor-acceptor pair. The process is frequently observed in molecules possessing an intramolecular hydrogen bond between a hydroxyl or amine group and a carbonyl or nitrogen atom, facilitating rapid proton translocation upon excitation with light.

Excited-state proton transfer (ESIPT) dynamics were investigated through computational modeling utilizing the test molecules o-Hydroxybenzaldehyde and 4-Amino-2-[1-(methylimino)ethyl]phenol. These molecules were selected to represent a range of ESIPT donor-acceptor characteristics and to facilitate comparison of proton transfer rates under varying electromagnetic field conditions. Simulations were performed to analyze the influence of plasmonic nanocavities on the ESIPT process, specifically examining how localized field enhancements affect the speed and efficiency of proton transfer within each molecule. The choice of these two molecules allows for a robust assessment of the general applicability of observed enhancements to diverse ESIPT systems.

Nanoparticle-on-Mirror (NIM) nanocavities function by supporting Localized Surface Plasmon Resonance (LSPR) to achieve significant electromagnetic field enhancement. This configuration consists of a metallic nanoparticle positioned in close proximity to a metallic mirror, creating a cavity where electromagnetic waves are confined and amplified. When light interacts with the structure, the free electrons within the metal collectively oscillate, resulting in a resonant enhancement of the electric field intensity within the cavity. The magnitude of this enhancement is dependent on the size, shape, and material of the nanoparticle, as well as the gap distance between the nanoparticle and the mirror; smaller gaps generally yield larger field enhancements. This intensified electromagnetic field is localized within the cavity volume, increasing the interaction between light and any molecules positioned within that space.

The Gaussian Spectral Density (GSD) model is employed to characterize the distribution of electromagnetic modes within the plasmonic nanocavity, providing a statistical description of the local density of optical states. This approach allows for the simulation of light-matter interactions by treating the electromagnetic field as a collection of harmonic oscillators with frequencies determined by the GSD. The model’s parameters are adjusted to reflect the geometry and material properties of the nanocavity, enabling accurate prediction of the resonant wavelengths and field enhancements. Simulations are performed with a time step of 0.1 atomic units (a.u.), balancing computational cost with the need to resolve the ultrafast dynamics of excited-state proton transfer processes and accurately capture the frequency-dependent response of the plasmonic system.

Simulations of the photoisomerization of HBA reveal that cavity coupling modulates the reaction dynamics, shifting the preferred vibrational modes and altering the H-O bond distances as evidenced by changes in the proton position expectation value between the <span class="katex-eq" data-katex-display="false">	ext{S}_0</span> and <span class="katex-eq" data-katex-display="false">	ext{S}_1</span> states.
Simulations of the photoisomerization of HBA reveal that cavity coupling modulates the reaction dynamics, shifting the preferred vibrational modes and altering the H-O bond distances as evidenced by changes in the proton position expectation value between the ext{S}_0 and ext{S}_1 states.

The Emergence of Polaritonic Control

The interaction between light and matter can yield fascinating hybrid states known as polaritons when a molecular excitation strongly couples with a plasmonic mode. This isn’t simply a shift in energy levels; it’s the creation of entirely new quasiparticles – part photon, part exciton – that inherit properties from both constituents. This strong coupling regime, distinct from weak interactions, dramatically alters the system’s behavior, evidenced by the creation of upper and lower polariton branches in the energy spectrum. These polaritonic states exhibit unique dispersion relations and lifetimes, influencing how energy flows within the molecule-cavity system. The formation of polaritons provides a pathway to control and manipulate molecular excitations, opening opportunities to engineer novel light-matter interactions and explore new regimes of chemical reactivity.

Rabi splitting provides compelling evidence for the strong coupling occurring within the plasmonic nanocavity. This phenomenon, manifested as a distinct separation in the energy spectrum, arises when the interaction between light and matter becomes significantly strong, leading to the creation of mixed light-matter states – polaritons. Observations reveal a substantial Rabi splitting of 30-50 meV, achieved with a surprisingly small ensemble of just nine molecules confined within the cavity. The magnitude of this splitting directly indicates the strength of the light-matter interaction, far exceeding the energy dissipation rates, and confirms that the system is operating in the strong coupling regime where conventional perturbative treatments are no longer valid. This clear spectroscopic signature not only validates the formation of polaritons but also opens possibilities for tailoring molecular excited-state dynamics through cavity quantum electrodynamics.

The creation of polaritons – hybrid light-matter quasiparticles – fundamentally alters how molecules respond to light, with significant consequences for excited-state processes like proton transfer. When a molecule strongly interacts with a plasmonic cavity, the resulting polariton states exhibit a modified energy landscape and altered spatial distribution of electronic density. This impacts the timescale and efficiency of proton transfer, potentially either accelerating the process by providing new reaction pathways or hindering it through destructive interference effects. Investigations reveal that the extent of this modification depends critically on the strength of the light-matter coupling and the specific characteristics of the polariton states, opening avenues for precise control over chemical reaction rates and yields through manipulation of the surrounding nanophotonic environment.

Recent research establishes the capacity to accurately model and anticipate how a challenging nanocavity environment – specifically, one characterized by energy loss and multiple resonant modes – influences the intricate process of excited-state proton transfer. Utilizing a cavity loss rate of 50 femtoseconds, simulations reveal the substantial potential of plasmonic nanocavities not merely as passive observation tools, but as active regulators of fundamental chemical reactions. This level of control opens avenues for manipulating molecular behavior at an unprecedented scale, suggesting future applications in areas such as light-harvesting, photocatalysis, and the design of novel molecular devices where precise control over proton transfer is paramount.

The study illuminates a fundamental truth about complex systems – their inherent impermanence. Just as infrastructure succumbs to erosion, the molecular dynamics within plasmonic nanocavities are not static, but rather a transient state influenced by energy loss and the environment. Kapitsa observed, “It is in the nature of things that everything changes.” This aligns with the research’s focus on non-adiabatic dynamics – the system’s evolution isn’t smooth, but punctuated by shifts responding to the cavity’s characteristics. The tracking of excited-state proton transfer, while a significant achievement, reveals a fleeting moment within a larger cycle of energy exchange and decay, a reminder that even strong coupling is temporary.

The Long View

The demonstrated linkage between excited-state proton transfer and lossy plasmonic cavities does not represent an arrival, but a reorientation. Every abstraction carries the weight of the past, and here, the abstraction is the very notion of controlling quantum dynamics. While tracking these events within a cavity is a technical achievement, the question is not whether observation influences outcome-it always does-but whether the influence is meaningfully durable. The inherent decay within both the molecular system and the plasmonic structure dictates that any imposed ‘control’ is transient, a fleeting alignment before inevitable dissipation.

Future work will undoubtedly refine the theoretical framework, perhaps incorporating more sophisticated treatments of decoherence or environmental effects. However, true progress lies in acknowledging the limitations of intervention. The focus should shift from forcing coherence to designing systems that accept decay gracefully, maximizing the useful lifetime of the quantum state within the constraints of the physical reality. Only slow change preserves resilience.

The real challenge is not to build ever-more-complex cavities, but to understand the fundamental limits of information retention within any dynamic system. The cavity, in the end, is merely a mirror, reflecting the inevitable entropic march forward – a beautiful, fleeting illusion of control against the backdrop of universal decay.


Original article: https://arxiv.org/pdf/2603.12373.pdf

Contact the author: https://www.linkedin.com/in/avetisyan/

See also:

2026-03-16 23:00