In total, 15 chromophore-solvent subsystems were constructed, optimized, and the LSC NAMD simulated for at least 1 fs with a 1 a.u. time step. CAM-B3LYP/6-31G(d) was also used for all simulations and S 1 and S 2 excited states were computed with TDDFT within the TDA. The last 40 time steps of all simulations were used in our benchmarks, averaging over the computational time needed to compute individual gradients and NAC vectors between the three states, as well as the total wall time for a NAMD time step. Each time step within a simulation used the ground state molecular orbitals and and excited state guess coefficients of the previous geometry as the initial guess for the next TeraChem calculation. To ensure that our statistics are not influenced by time steps that do not have these initial guesses available, we eliminated the first few sets of timings of each simulation from our averaging. We compared the performance using TCPBs to distribute work over multiple computers on the same network, running all 15 simulations on both on a single and dual node setup.