Molecular Vibrations Defy Simple Thermal Description in New Simulations

R. Schwengelbeck and colleagues investigate how disorder impacts vibrational dynamics within systems of molecules strongly coupled to optical cavities. Disorder generates non-Gaussian states of vibrational modes at the single-molecule level, challenging the assumption that these states can be accurately represented as thermal distributions. Matrix product state simulations reveal this effect persists even as the number of molecules increases, highlighting the limitations of semiclassical approximations such as the Ehrenfest and truncated Wigner methods in describing these quantum phenomena.

This establishes the key role of both disorder and inherently quantum behaviour in understanding nuclear dynamics within the emerging field of polaritonic chemistry. Vibrational modes within molecular systems exhibit non-Gaussian behaviour on short timescales, specifically at the single-molecule level. Using exact matrix product state simulations with a Holstein-Tavis-Cummings model, a framework for describing light–matter interactions, the authors found that disorder sustains this effect even as the number of molecules increases.

This means conventional descriptions of nuclear wavepackets, which treat them as thermally distributed, are inadequate; the simulations reveal that the vibrational states are genuinely non-thermal. Furthermore, comparisons with simpler, commonly used approximations such as the Ehrenfest method show that it only accurately predicts average behaviour for very large systems. The truncated Wigner approximation also fails to capture these non-Gaussian effects, highlighting the need for more sophisticated quantum mechanical treatments.

Non-classical vibrational behaviour persists in simplified polaritonic models

The field of polaritonic chemistry is receiving increasing attention as a means of using light to control chemical reactions. It relies on creating hybrid light–matter states called polaritons, which offer the potential to reshape molecular potential energy surfaces and modify reactivity. Polariton formation arises from the strong coupling of molecular excitations with photons within an optical cavity, leading to the creation of new quasiparticles with unique properties. This strong-coupling regime is typically defined by a condition in which the rate of interaction between molecules and the electromagnetic field exceeds the rates of dissipation in the system. Accurately modelling the behaviour of many molecules collectively coupled to light remains a key challenge, because current methods often simplify the quantum nature of these systems.

Detailed simulations reveal that disorder within molecular systems strongly influences vibrational behaviour when coupled with light. Non-classical vibrations, deviating from predictions made by simpler conventional models, were observed and emphasise the importance of fully quantum treatments. A Holstein–Tavis–Cummings model, representing N molecules interacting with light, demonstrates that static disorder sustains non-Gaussian vibrational states within molecular ensembles even as their size increases. This model is a cornerstone of the study, describing N molecules, each with a harmonic vibrational mode, interacting with a single mode of the electromagnetic field. The disorder introduced refers to variations in molecular properties, such as vibrational frequencies or coupling strengths, across the ensemble.

The persistence of these non-classical vibrations indicates that conventional approximations of nuclear wavepackets are insufficient for accurately modelling light–matter interactions. Consequently, this challenges the validity of widely used semiclassical methods such as the Ehrenfest and truncated Wigner approximations in describing cavity-modified nuclear dynamics. The Ehrenfest method, a classical trajectory approach, treats nuclei as classical particles moving on a potential energy surface that incorporates the effects of the electromagnetic field. While computationally efficient, it neglects quantum coherence and entanglement, leading to inaccuracies in strongly coupled systems. The truncated Wigner approximation attempts to incorporate some quantum effects by representing the quantum density matrix as a Gaussian distribution, but it fails to capture non-Gaussian features such as those observed in this study.

Further analysis explored how model parameters, including the strength of light–matter coupling and the degree of static disorder, affect the observed vibrational behaviour, revealing a complex interplay between these factors.

The researchers employed exact matrix product state (MPS) simulations, a powerful numerical technique for studying strongly correlated quantum systems. MPS provides an efficient representation of the many-body wavefunction, enabling accurate calculations of system dynamics. This method is particularly well-suited for one-dimensional systems, such as the linear chain of molecules considered in this study. The simulations focused on the evolution of the vibrational wavepacket following incoherent excitation by a photon. This initial excitation introduces energy into the system, driving the vibrational modes and allowing the researchers to probe their subsequent dynamics.

Despite the model’s simplification of real molecular systems, the persistence of non-classical behaviour—vibrational states deviating from simple thermal descriptions—is significant. Approximations that work well for larger systems fail in this regime, and light-induced changes to molecular vibrations are not easily predicted by conventional, less detailed methods. The implications extend to the design of future experiments, highlighting the need to account for disorder effects when probing polaritonic chemistry. Understanding the role of disorder is crucial for interpreting experimental results and developing strategies to control molecular vibrations using light.

The observed non-Gaussian behaviour suggests that vibrational energy is not evenly distributed among modes, leading to the formation of localized wavepackets that may influence chemical reaction rates. This work provides a foundation for developing more accurate theoretical models and designing experiments to explore the full potential of polaritonic chemistry.

The fact that these non-classical effects remain robust even as the number of molecules N increases is particularly significant. It suggests that these quantum phenomena are not merely artefacts of small system sizes but are inherent features of the underlying physics. This robustness has implications for the scalability of polaritonic chemistry, indicating that meaningful control over chemical reactions may be possible even in larger, more complex systems. The researchers demonstrate this persistence up to a certain value of N, and further investigations are needed to determine the limits of this behaviour.

Researchers found that disorder creates non-classical vibrational states within molecules strongly coupled to light. This is significant because it shows that standard methods for modelling molecular vibrations, such as the Ehrenfest and truncated Wigner approximations, struggle to accurately represent these effects, particularly for smaller systems. The study used simulations of the Holstein–Tavis–Cummings model to show that this non-Gaussian behaviour persists even as the number of molecules increases. The authors suggest further investigation is needed to fully understand the limits of this observed robustness.