Initialization with a Fock State Cavity Mode in Real-Time Nuclear--Electronic Orbital Polariton Dynamics

Here is an explanation of the paper using simple language and creative analogies.

The Big Picture: Dancing with Light

Imagine a molecule (like a tiny, vibrating spring) sitting inside a mirrored box (a "cavity"). Usually, scientists study what happens when they shine a bright, steady laser light into this box. The light bounces back and forth, hitting the molecule, making it vibrate, which in turn changes the light. They are "dancing" together. This is called polariton dynamics.

For a long time, scientists have used two ways to simulate this dance on computers:

The "Classical" Way: Treating the light like a smooth, continuous wave (like water in a hose).
The "Quantum" Way: Treating the light as individual packets of energy called photons (like distinct marbles).

The big question this paper asks is: Does it matter if we treat the light as a smooth wave or as individual marbles? Specifically, what happens if we start the dance with exactly one photon (a single marble) instead of a steady stream of light?

The Experiment: The "One-Marble" Start

The researchers used a molecule called HCN (Hydrogen-Cyanide) and put it in a virtual box. They set up two different computer simulations:

Simulation A (The "Mean-Field" Method): This method assumes the molecule and the light are dancing separately but influencing each other. It's like two dancers who can see each other but are never actually holding hands or tangled up.
Simulation B (The "Full-Quantum" Method): This method allows the molecule and the light to get "entangled." Imagine them becoming a single, inseparable unit where the state of one instantly affects the other, no matter how far apart they seem.

They started both simulations with the molecule at rest and the light box containing exactly one photon (a Fock state).

The Results: The Silent vs. The Loud Dance

1. The "Mean-Field" Result (Simulation A): The Silent Room

When they ran the simulation where the light and matter couldn't get entangled, nothing happened.

The Analogy: Imagine you throw a single marble at a wall. If the wall is perfectly rigid and you aren't allowed to "feel" the impact, the marble just sits there, and the wall doesn't move.
The Science: Because the starting state (one photon) has a "zero average position" (the photon isn't really "here" or "there" in a classical sense), the computer thought there was no force to push the molecule. The molecule stayed still, and the light stayed still. No energy was exchanged. No "polariton" (the dance partner) was formed.

2. The "Full-Quantum" Result (Simulation B): The Hidden Dance

When they ran the simulation that allowed entanglement, something magical happened, but it was invisible to the naked eye.

The Analogy: Imagine two dancers spinning so fast and tightly together that if you look at their hands (the average position), they seem to be standing perfectly still. But if you look at their speed or the energy of their spin, they are vibrating wildly.
The Science:
- The "Hands" Stay Still: Just like in the first simulation, the average position of the light and the molecule didn't move. If you only looked at the "average," you'd think nothing was happening.
- The "Energy" Vibrates: However, when the scientists looked at the squares of the movement (energy and intensity), they saw a wild oscillation. The molecule and the light were exchanging energy back and forth rapidly.
- The "Entanglement" Meter: They measured a "quantum entanglement score" (Von Neumann entropy). This score went up and down, proving that the molecule and the light had become a single, quantum-mechanical entity. They were truly dancing together.

The "Even vs. Odd" Mystery

The researchers noticed a weird pattern that only the "Full-Quantum" method could see:

Odd Powers (1st, 3rd, 5th): These represent the "average" position. They stayed flat (no movement).
Even Powers (2nd, 4th, 6th): These represent the "energy" or "intensity." These oscillated wildly.

The Metaphor: Think of a pendulum swinging.

If you look at its position (left or right), it moves back and forth.
But if you look at a Fock state (a specific quantum state), it's like the pendulum is swinging so perfectly symmetrically that its average position is always zero. It's not "left" or "right" on average; it's just "swinging."
However, if you look at the square of the position (how far it is from the center, regardless of direction), you see it moving away and coming back. That's the "even power" oscillation.

Why This Matters

This paper proves that how you start the experiment changes the physics.

Classical vs. Quantum: If you treat light classically (or use the simplified "Mean-Field" method), you might miss the entire phenomenon. You would conclude that "nothing happens" with a single photon.
The Hidden Reality: The "Full-Quantum" method reveals that even with a single photon, the molecule and light are interacting and forming a new hybrid particle (a polariton). They just do it in a way that doesn't look like a simple back-and-forth motion.
The Takeaway: To understand the future of "Polariton Chemistry" (using light to control chemical reactions), we cannot just use classical physics. We need the full quantum toolkit to see the "hidden dance" that happens when we deal with very small amounts of light.

Summary

Imagine a room with a single person (the molecule) and a single ghost (the photon).

Old View: The ghost is invisible and has no weight, so the person doesn't move.
New View: The ghost and the person are holding hands so tightly they become one super-being. They are vibrating with energy, but because they are moving in perfect sync, they don't seem to go anywhere. If you only look at their location, you miss the magic. If you look at their energy, you see the dance.

This paper teaches us that in the quantum world, what you don't see (the average position) isn't always what's happening. The real action is in the hidden connections (entanglement) and the energy fluctuations.

Here is a detailed technical summary of the paper "Initialization with a Fock State Cavity Mode in Real-Time Nuclear–Electronic Orbital Polariton Dynamics" by Millan F. Welman and Sharon Hammes-Schiffer.

1. Problem Statement

The field of molecular polaritons investigates the strong coupling between light (cavity modes) and matter (molecules) to control chemical dynamics. A central challenge is identifying physical phenomena that strictly require a quantum mechanical treatment of the electromagnetic field, as opposed to a semiclassical or classical treatment.

Previous theoretical work often initializes cavity modes in coherent states, which possess non-zero expectation values for position and momentum, mimicking classical oscillators. This allows for a direct comparison with semiclassical methods. However, Fock states (number states), representing a specific integer number of photons, are a fundamentally quantum state with zero expectation values for position and momentum operators. The authors investigate whether initializing a cavity mode in a Fock state (specifically a single-photon state) reveals unique quantum dynamics that are absent in classical or coherent-state simulations, particularly regarding the formation of polaritons and light-matter entanglement.

2. Methodology

The study employs two first-principles dynamical simulation methods developed by the authors, both based on Real-Time Nuclear–Electronic Orbital Time-Dependent Density Functional Theory (RT-NEO-TDDFT). In this framework, specific nuclei (protons) are treated quantum mechanically alongside electrons.

System: A single HCN molecule under Vibrational Strong Coupling (VSC) with an x-polarized cavity mode resonant with the proton bending mode ( $\omega_c \approx 2803 \text{ cm}^{-1}$ ).
Initial Conditions:
- Molecule: Ground state (SCF).
- Cavity Mode: Fock state $|1\rangle$ (single photon), meaning the expectation values of the coordinate ( $\hat{q}$ ) and momentum are zero.
Two Approaches Compared:
1. Mean-Field Quantum (mfq-RT-NEO): Treats the light and matter subsystems with separate density matrices coupled via mean fields. It neglects light-matter entanglement.
2. Full-Quantum (fq-RT-NEO): Propagates a joint light-matter density matrix, explicitly allowing for light-matter entanglement.
Observables: Time evolution of the cavity coordinate ( $q$ ), molecular dipole ( $\mu$ ), their powers ( $q^n, \mu^n$ ), and the von Neumann entropy ( $S(t)$ ) as a measure of entanglement.
Theoretical Analysis: The authors utilize a Quantum Rabi Model (extending the Jaynes-Cummings model) to analytically explain the simulation results, specifically analyzing the block-diagonal structure of the Hamiltonian in dressed-state bases.

3. Key Contributions

Demonstration of Methodological Divergence: The paper highlights a critical failure of mean-field approaches (and by extension, classical field treatments) to describe polariton formation when initialized in a Fock state, whereas full-quantum methods succeed.
Discovery of "Hidden" Polariton Dynamics: The authors identify that while first-order observables ( $q$ and $\mu$ ) remain static in a Fock state initialization, second-order observables ( $q^2$ and $\mu^2$ ) and entanglement entropy exhibit significant oscillations, signaling polariton formation.
Parity Rule for Operator Expectations: The study establishes a general trend where expectation values of odd powers of the operators remain constant, while even powers oscillate, a phenomenon explained by the symmetry properties of the Hamiltonian in the dressed-state basis.
Spectral Analysis: The power spectra of the oscillating even-power operators reveal a Rabi splitting ( $\Omega_R$ ) and a doublet centered at $2\omega_c $, rather than the canonical splitting at$ \omega_c$ seen in coherent state simulations.

4. Key Results

A. Mean-Field Quantum (mfq-RT-NEO) Results

No Dynamics: When initialized with a Fock state, the expectation values of the cavity coordinate $q(t)$ and the molecular dipole $\mu(t)$ remain constant (zero) throughout the simulation.
No Polariton Formation: Because the mean-field equations rely on the coupling of expectation values (which are zero), no energy exchange occurs. The method predicts no polariton formation, effectively failing to capture the physics of the single-photon regime.

B. Full-Quantum (fq-RT-NEO) Results

Static First-Order Observables: Surprisingly, even in the full-quantum treatment, the expectation values $q(t)$ and $\mu(t)$ remain constant. This is because the initial state ( $|1\rangle \otimes |0\rangle_{mol}$ ) is a superposition of only Jaynes-Cummings eigenstates, which do not couple to the counter-rotating terms required to oscillate the first-order operators in this specific basis.
Dynamic Entanglement: The von Neumann entropy $S(t)$ oscillates significantly (reaching ~48% of the theoretical maximum), indicating strong light-matter entanglement and coherent energy exchange.
Oscillating Even Powers: The expectation values of the squared operators ( $q^2(t)$ $q^{2} (t)$ and $\mu^2(t)$ $μ^{2} (t)$ ) exhibit clear oscillations with two timescales:
- A slow oscillation (~59 fs) corresponding to the Rabi frequency ( $\Omega_R \approx 565 \text{ cm}^{-1}$ ).
- A fast oscillation (~6 fs) corresponding to $2\omega_c$.
Odd vs. Even Powers: A systematic analysis shows that expectation values for odd powers ( $n=1, 3, 5...$ ) are constant, while even powers ( $n=2, 4, 6...$ ) oscillate.

C. Spectral Features

The power spectrum of $q^2(t)$ reveals:

A single peak at the Rabi splitting frequency ( $\Omega_R$ ).
A pair of peaks centered at **$2\omega_c $** (twice the cavity frequency) with a splitting of$ \Omega_R$.
No features at the bare cavity frequency $\omega_c$ .

D. Theoretical Explanation

Using a Quantum Rabi Model with an extended basis (treating the cavity as a 3-level system and molecule as a 2-level system), the authors explain:

The Hamiltonian splits into Jaynes-Cummings (JC) and Counter-Rotating (CR) blocks based on the parity of the difference between photon number and molecular excitation ( $F-M$ ).
The initial Fock state $|10\rangle$ lies entirely within the JC block.
First-order operators ( $q, \mu$ ) only have matrix elements between JC and CR blocks; since the state has no CR component, these do not oscillate.
Second-order operators ( $q^2, \mu^2$ ) have non-zero matrix elements within the JC block, allowing for oscillations between the upper and lower polariton states.
The spectral peaks at $2\omega_c $arise from transitions involving higher excited states (e.g.,$ |21\rangle$) accessible via the squared operators.

5. Significance

Beyond Classical Electrodynamics: The results demonstrate that initializing a cavity in a Fock state leads to physical phenomena (entanglement and energy exchange) that are invisible to classical or mean-field treatments. This validates the necessity of full-quantum electrodynamics for certain regimes of polariton chemistry.
Limitations of the Jaynes-Cummings Model: The study shows that the standard two-level, rotating-wave approximation (RWA) picture is insufficient. Capturing the dynamics of even-power operators and the specific spectral features requires going beyond RWA and utilizing larger basis sets (multi-level systems).
New Diagnostic Tools: The paper suggests that in Fock-state initialized systems, entanglement entropy and higher-order correlation functions (like $\langle q^2 \rangle$ ) are superior indicators of polariton formation compared to standard dipole or field coordinate oscillations.
Implications for Experiment: While current experiments often use coherent states (lasers), the findings suggest that future experiments involving single-photon sources or specific quantum state preparations could reveal unique quantum dynamics not predicted by semiclassical models.

In conclusion, this work establishes that the choice of initial quantum state (Fock vs. Coherent) fundamentally alters the observability of polariton dynamics, necessitating full-quantum treatments to capture the rich physics of light-matter entanglement in the single-photon regime.