arxiv: 2512.03277 · v1 · submitted 2025-12-02 · ⚛️ physics.chem-ph · cond-mat.mes-hall

On-the-Fly Cavity-Molecular Dynamics of Vibrational Polaritons

Sachith Wickramasinghe , Amirhosein Amini , Arkajit Mandal This is my paper

Pith reviewed 2026-05-17 01:48 UTC · model grok-4.3

classification ⚛️ physics.chem-ph cond-mat.mes-hall

keywords vibrational polaritonscavity molecular dynamicsdensity functional tight-bindinglight-matter interactionMulliken chargesBorn chargesFabry-Perot cavitypolaritonic spectra

0 comments

The pith

Mulliken charges can stand in for Born charges to compute qualitatively accurate linear spectra of vibrational polaritons when molecular nonlinearity stays low.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops an on-the-fly molecular dynamics method that couples density functional tight-binding to a light-matter Hamiltonian beyond the long-wavelength limit inside a Fabry-Perot cavity. It demonstrates that the cheaper Mulliken charges produce linear spectra matching those from expensive Born charges under limited nonlinearity from vibrations. This substitution reduces computational cost while preserving essential features of the hybrid light-matter states. The work tests the scheme on water to generate angle-resolved spectra and releases an open-source implementation for broader use in cavity-modified chemistry.

Core claim

Combining DFTB with a real-space light-matter Hamiltonian allows propagation of vibrational polariton dynamics, and Mulliken charges can replace Born charges to obtain qualitatively accurate linear spectra when the nonlinearity arising from molecular vibrations remains limited; the same replacement produces spurious heating that disqualifies it for studies of energy transport or chemical dynamics.

What carries the argument

The sparse real-space representation of light-matter interactions that supports a parallelized propagation scheme with lightweight inter-CPU communication.

If this is right

Angle-resolved polaritonic spectra become computable for molecules such as water at modest cost.
The method enables routine exploration of cavity effects on vibrational spectra without full Born-charge overhead.
The approximation remains unsuitable for following energy transport or reaction dynamics because it introduces artificial heating.
Open-source code allows direct testing on other molecular systems coupled to confined radiation.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The charge replacement may extend naturally to other semi-empirical electronic-structure methods used in polariton simulations.
Adding explicit higher-order terms to the interaction could remove the heating artifact and broaden applicability to reactive processes.

Load-bearing premise

Nonlinearity in the light-matter interaction term stays small enough that Mulliken charges still capture the essential coupling physics.

What would settle it

A side-by-side calculation of linear spectra for a molecule with strong vibrational nonlinearity, showing clear mismatches in peak positions or intensities between Born-charge and Mulliken-charge versions.

Figures

Figures reproduced from arXiv: 2512.03277 by Amirhosein Amini, Arkajit Mandal, Sachith Wickramasinghe.

**Figure 1.** Figure 1: Schematic illustration of light-matter interactions and the real-reciprocal used here for propagating vibro-polariton dynamics. Schematic illustration of light-matter interaction (a) in the standard multi-mode dipole gauge description and (b) in real space. (c) Schematic illustration of the algorithm used here and implemented in our computational package CavOTF. 48 direction of the cavity field which has b… view at source ↗

**Figure 2.** Figure 2: (a) Born charges distribution of Oxygen and Hydrogen atoms in liquid water and (b) their timedependent fluctuations. obtained as an = N X−1 k=0 r ωk 2N q˜k(t + δt) + ip˜k(t + δt) ωk e ik·Rn . (19) The updated values are then sent to all clients. Each client receives only its corresponding pair (qn, p˜n). Step 5. Repeat. Continue from Step 1 for the next time step. Computational details. Molecular dyna… view at source ↗

**Figure 3.** Figure 3: c shows the simulated molecular spectra featuring a peak corresponding to the H-O-H bending mode in bulk water outside of the cavity (schematically illustrated in Fig. 3a). This particular peak appears at ∼0.19 eV with a narrow profile and is in good agreement with past experimental and theoretical works. 61,62 In Fig. 3d, we illustrate the cavity modified molecular spectra when considering a single rad… view at source ↗

**Figure 4.** Figure 4: Cavity modified molecular spectra and angle resolved photonic spectra under vibrational strong coupling. Panels (a)–(d) show cavity modified molecular spectra at various light–matter coupling strengths for η0 = 0.0001 with ω0 = 0.19 eV, and panels (f)–(i) show the corresponding spectra for η0 = 0.00034 with ω0 = 0.43 eV. Angle resolved photonic spectra are shown in panel (e) for ω0 = 0.19 eV and in panel (… view at source ↗

**Figure 5.** Figure 5: Angle-resolved photonic spectra under vibrational strong coupling when using Mulliken charges vs Born charges at (a) ω0 = 0.16 eV, ω0 = 0.19 eV, ω0 = 0.22 eV, and ω0 = 0.43 eV. Here we use the light-matter coupling η = ηc · ω 3/2 0 with ηc = 0.00105. with the k = 0 cavity mode tuned to 0.19 eV. Thus, k = 0 cavity mode is in resonance with the H–O–H bending vibration, forming the upper and lower polariton … view at source ↗

read the original abstract

In this work, we combine the density functional tight-binding (DFTB) approach with a light-matter Hamiltonian beyond the long-wavelength approximation to propagate the dynamics of vibrational polaritons formed by coupling molecular vibrations to confined radiation inside a Fabry-P\'{e}rot optical cavity. Here, we develop a parallelized propagation scheme with lightweight inter-CPU communication by exploiting the sparse nature of the light-matter interactions in the real space representation. We find that the computationally expensive Born charges required for our propagation can be replaced with the computationally inexpensive Mulliken charges to obtain qualitatively accurate linear spectra especially when the nonlinearity (arising from molecular vibrations) of the light-matter interaction term is not substantial. However, the same approach may not be suitable to be used for studying cavity modification of energy transport or chemical dynamics as this approximation leads to spurious heating of the light-matter hybrid system. We demonstrate the utility of this on-the-fly approach to compute angle resolved polaritonic spectra of water. We implement our approach as an open-source computational package, CavOTF, which is available on GitHub.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives a practical parallelized DFTB scheme for vibrational polariton dynamics and shows Mulliken charges can work for linear spectra under limited nonlinearity, but the regime for that swap stays loosely defined.

read the letter

The main thing to know is that this work builds an on-the-fly propagation method that pairs density functional tight-binding with a real-space light-matter Hamiltonian beyond the long-wavelength limit, then shows that inexpensive Mulliken charges can stand in for Born charges to get qualitatively decent linear spectra when vibrational nonlinearity stays modest. They also flag that the same swap produces spurious heating and should not be used for transport or reactivity studies. The parallel scheme exploits sparse interactions to keep inter-CPU communication light, which is a sensible engineering choice for scaling. Releasing the CavOTF code on GitHub and running an angle-resolved water example adds concrete value that others can test directly. The distinction between spectra and dynamics use cases is a fair warning. The soft spot sits in the approximation itself. The abstract ties the Mulliken replacement to cases where nonlinearity is not substantial, yet no explicit threshold, deviation metric, or conservation check appears to bound that regime. The heating result hints that energy balance may drift even in milder settings, which could quietly shift spectral features. Without tighter numerical validation or error bars across coupling strengths, it is hard to know how far the shortcut reliably extends. This paper is aimed at computational groups working on cavity-modified molecular dynamics who need something lighter than full ab initio treatments. A reader who builds or adapts simulation tools for polariton chemistry would find the implementation details and charge comparison worth examining. I would send it for peer review so referees can verify the numerical evidence and ask for clearer limits on the approximation.

Referee Report

2 major / 2 minor

Summary. The manuscript develops an on-the-fly cavity-molecular dynamics scheme that couples density-functional tight-binding (DFTB) to a light-matter Hamiltonian beyond the long-wavelength approximation. It exploits sparsity for a parallelized propagation algorithm and reports that computationally cheap Mulliken charges can replace expensive Born charges to produce qualitatively accurate linear polariton spectra provided the vibrational nonlinearity in the light-matter interaction remains limited; the same replacement produces spurious heating that precludes reliable transport or reactive dynamics. The approach is demonstrated on angle-resolved spectra of water and released as the open-source CavOTF package.

Significance. If the central approximation holds under a clearly delineated regime, the method supplies an efficient route to cavity-modified vibrational spectra for systems too large for ab initio treatments, with the open-source implementation and parallelization constituting concrete strengths. The reported distinction between linear-response accuracy and heating artifacts in non-equilibrium observables is potentially useful for the polaritonics community.

major comments (2)

Abstract and §3 (charge-approximation discussion): the central claim that Mulliken charges suffice for qualitatively accurate linear spectra 'especially when the nonlinearity ... is not substantial' is load-bearing yet lacks any quantitative metric, threshold, or deviation measure (e.g., integrated spectral difference or bound on quadratic interaction terms) that would delineate the validity regime; without such a criterion the statement remains difficult to falsify or apply.
§4 (heating and transport results): the observation of spurious heating under the Mulliken approximation is presented as disqualifying for energy-transport or chemical-dynamics studies, but no explicit energy-conservation diagnostic, drift rate, or comparison of conserved quantities versus Born-charge reference trajectories is supplied, leaving open whether the artifact is systematic or controllable.

minor comments (2)

The parallelization scheme and sparsity exploitation are described at a high level; a short pseudocode or timing table would clarify the claimed lightweight inter-CPU communication.
Notation for the light-matter interaction term (real-space representation) should be cross-referenced to the specific equation defining the dipole or charge operator to avoid ambiguity for readers implementing the method.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comments. We address each major comment below and propose revisions to strengthen the presentation of our results.

read point-by-point responses

Referee: Abstract and §3 (charge-approximation discussion): the central claim that Mulliken charges suffice for qualitatively accurate linear spectra 'especially when the nonlinearity ... is not substantial' is load-bearing yet lacks any quantitative metric, threshold, or deviation measure (e.g., integrated spectral difference or bound on quadratic interaction terms) that would delineate the validity regime; without such a criterion the statement remains difficult to falsify or apply.

Authors: We agree that the validity regime would be better delineated with a quantitative metric. In the revised manuscript we will add a quantitative comparison, specifically the integrated absolute difference between the linear spectra obtained with Mulliken and Born charges across the systems studied. We will also report the relative magnitude of the quadratic terms in the light-matter interaction for the cases examined, thereby providing an operational threshold for when the nonlinearity remains limited. revision: yes
Referee: §4 (heating and transport results): the observation of spurious heating under the Mulliken approximation is presented as disqualifying for energy-transport or chemical-dynamics studies, but no explicit energy-conservation diagnostic, drift rate, or comparison of conserved quantities versus Born-charge reference trajectories is supplied, leaving open whether the artifact is systematic or controllable.

Authors: We acknowledge that explicit diagnostics would make the heating artifact more transparent. In the revision we will include time traces of the total energy for representative trajectories using both Mulliken and Born charges, together with the observed energy drift rates. We will also compare additional conserved quantities (such as the cavity-field energy) between the two approximations to illustrate the systematic character of the spurious heating. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation rests on explicit propagation and spectral comparisons

full rationale

The paper develops a new parallelized on-the-fly propagation scheme that couples DFTB to a light-matter Hamiltonian beyond the long-wavelength limit, exploiting sparsity for inter-CPU communication. The central finding—that Mulliken charges can substitute for Born charges to yield qualitatively accurate linear spectra when vibrational nonlinearity in the interaction term remains limited—is presented as an empirical outcome of their simulations on water, not as a definitional equivalence, a fitted parameter renamed as a prediction, or a result forced by self-citation. No load-bearing uniqueness theorems, ansatzes smuggled via prior work, or renamings of known patterns appear in the derivation chain; the results are benchmarked against computed angle-resolved polaritonic spectra and are therefore self-contained against external computational checks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard assumptions of DFTB electronic structure and the validity of the chosen light-matter Hamiltonian form; no new free parameters or invented entities are explicitly introduced in the abstract.

axioms (2)

domain assumption DFTB provides a sufficiently accurate description of molecular vibrations and charges for the systems studied
Invoked implicitly when using DFTB for on-the-fly propagation of polariton dynamics.
domain assumption The light-matter Hamiltonian beyond the long-wavelength approximation accurately captures cavity confinement effects
Central to the propagation scheme described.

pith-pipeline@v0.9.0 · 5497 in / 1374 out tokens · 31317 ms · 2026-05-17T01:48:12.273296+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We find that the computationally expensive Born charges required for our propagation can be replaced with the computationally inexpensive Mulliken charges to obtain qualitatively accurate linear spectra especially when the nonlinearity (arising from molecular vibrations) of the light-matter interaction term is not substantial.
IndisputableMonolith/Foundation/BlackBodyRadiationDeep.lean off_match_positive unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

However, the same approach may not be suitable to be used for studying cavity modification of energy transport or chemical dynamics as this approximation leads to spurious heating of the light-matter hybrid system.

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

7 extracted references · 7 canonical work pages

[1]

F.; Recabal, F.; Herrera, F.; Simpkins, B

(2) Ahn, W.; Triana, J. F.; Recabal, F.; Herrera, F.; Simpkins, B. S. Modification of ground-state chemical reactivity via light–matter coherence in infrared cavities.Science2023,380, 1165–1168. (3) Nagarajan, K.; Thomas, A.; Ebbesen, T. W. Chemistry under vibrational strong coupling. Journal of the American Chemical Society2021, 143, 16877–16889. (4) Yin...

work page 2023
[2]

Modification of thermal chemical rates in a cavity via resonant effects in the collective regime.The Journal of Physical Chemistry Letters2023,14, 8397–8404

(28) Sun, J.; Vendrell, O. Modification of thermal chemical rates in a cavity via resonant effects in the collective regime.The Journal of Physical Chemistry Letters2023,14, 8397–8404. (29) Li, T. E.; Nitzan, A.; Subotnik, J. E. Collective vi- brational strong coupling effects on molecular vi- brational relaxation and energy transfer: Numer- ical insights...

work page 2021
[3]

S.; Keshavamurthy, S

(36) Mondal, S.; Wang, D. S.; Keshavamurthy, S. Dis- sociation dynamics of a diatomic molecule in an optical cavity.The Journal of Chemical Physics 2022,157. (37) Liu, W.; Chen, J.; Dou, W. Polaritons under ex- tensive disordered gas-phase molecular rotations in a Fabry–P´ erot cavity.The Journal of Physical Chemistry C2024,128, 12544–12550. (38) Li, X.; ...

work page 2022
[4]

Theory of vibra- tional polariton chemistry in the collective cou- pling regime.The Journal of Chemical Physics 2022,156

(39) Mandal, A.; Li, X.; Huo, P. Theory of vibra- tional polariton chemistry in the collective cou- pling regime.The Journal of Chemical Physics 2022,156. (40) Campos-Gonzalez-Angulo, J. A.; Ribeiro, R. F.; Yuen-Zhou, J. Resonant catalysis of thermally ac- tivated chemical reactions with vibrational polari- tons.Nature communications2019,10,

work page 2022
[5]

F.; Du, M.; Chen, L.; Yang, Z.; Wang, J.; Yuen-Zhou, J.; Xiong, W

10 (41) Xiang, B.; Ribeiro, R. F.; Du, M.; Chen, L.; Yang, Z.; Wang, J.; Yuen-Zhou, J.; Xiong, W. In- termolecular vibrational energy transfer enabled by microcavity strong light–matter coupling.Sci- ence2020,368, 665–667. (42) Brawley, Z. T.; Pannir-Sivajothi, S.; Yim, J. E.; Poh, Y. R.; Yuen-Zhou, J.; Sheldon, M. Vibra- tional weak and strong coupling m...

work page 2021
[6]

J.; Morokuma, K

(55) Zheng, G.; Lundberg, M.; Jakowski, J.; Vreven, T.; Frisch, M. J.; Morokuma, K. Implementation and benchmark tests of the DFTB method and its application in the ONIOM method.International Journal of Quantum Chemistry2009,109, 1841–1854. (56) Wickramasinghe, S.; Hoehn, A.; Wetthas- inghe, S. T.; Lin, H.; Wang, Q.; Jakowski, J.; Rassolov, V.; Tang, C.; ...

work page 2012
[7]

B.; Fournier, J

(62) Carpenter, W. B.; Fournier, J. A.; Biswas, R.; Voth, G. A.; Tokmakoff, A. Delocalization and stretch-bend mixing of the HOH bend in liquid wa- ter.The Journal of Chemical Physics2017,147. (63) Vinaykin, M.; Benderskii, A. V. Vibrational sum-frequency spectrum of the water bend at the air/water interface.The Journal of Physical Chemistry Letters2012,3...

work page 2023