“…Recently, Zhang, Nelson, and Tretiak implemented analytic nuclear gradients on the dipole and simulated the photoexcited dynamics of the stilbene molecule . In this work, the authors modified the NEXMD software package ,,− to include the pQED Hamiltonian (see Section ) at the Jaynes-Cummings level with all proper gradients required for this Hamiltonian (i.e., without DSE and making the rotating wave approximation). Additionally, the gradients on the potential energy surfaces, nonadiabatic couplings, and dipole gradients were achieved analytically at the TD-AM1 level of theory in the collective electronic oscillator (CEO) framework. ,, Most importantly, the nuclear gradient on the bare transition dipole between the ground and excited electronic states was computed as μ 0α = Tr[ μ̂ X̂ 0α ] in the atomic orbital { o , v } basis and can be understood as, ∂ μ 0 α ∂ R j = ∑ o v ∂ μ o v ∂ boldR j X v o 0 α + μ o v ∂ X v o 0 α ∂ R j where X vo 0α is the transition density matrix similar to that found in eq between the ground and α th excited electronic state in the CIS-approximation ,,, (see additional discussion in Sec.…”