Initialization with a Fock state cavity mode in real-time nuclear-electronic orbital polariton dynamics

Molecular polaritons have drawn great interest in recent years as a possible avenue for providing optical control over chemical dynamics. A central challenge in the field is to identify physical phenomena that require a quantum rather than a classical treatment of electrodynamics. In this work, we use our recently developed mean-field quantum (mfq) and full-quantum (fq) real-time nuclear-electronic orbital (RT-NEO) time-dependent density functional theory methods to simulate polaritonic dynamics for a molecule under vibrational strong coupling when a quantized cavity mode is initialized in a Fock state rather than a coherent state. Our previous work showed that a coherent state initial condition for the cavity mode leads to polariton formation for both the mfq-RT-NEO and fq-RT-NEO methods. Herein, we show that the mfq-RT-NEO method, which does not allow light-matter entanglement, does not predict polariton formation for a Fock state initial condition. Similar to the mfq-RT-NEO method, the fq-RT-NEO method does not predict oscillations of the cavity mode coordinate and molecular dipole operator expectation values for a Fock state initial condition. However, the fq-RT-NEO method does predict oscillations of the expectation values of even powers of these operators as well as light-matter entanglement, implicating polariton formation with a Fock state initial condition. All these observations can be explained with model systems. These results suggest that using a quantized cavity mode initial condition that does not have a direct analogy to an initial condition in classical electrodynamics can lead to physical phenomena that can only be described by a quantum treatment of the cavity mode.