The TeraChem software package is built to take advantage of GPU parallelization, and offers significant speed up of quantum electronic structure calculations compared to CPU-based packages. A variety of excited state electronic structure methods have been implemented within TeraChem, including SA-CASSCF, full-CI, CASCI, and TDDFT. These GPUaccelerated implementations have been used with both AIMS and surface hopping NAMD simulations. Recent extension of TeraChem to a protocol buffer clientserver model has allowed for straightforward interface between TeraChem's electronic structure capabilities and any external code for on-the-fly dynamics. Motivated by the success of the recent on-the-fly nonadiabatic LSC simulations, here we introduce the PySCES code (Python code for Semiclassical Dynamics with onthe-fly Electronic Structure) that interfaces the LSC method with the TeraChem electronic structure code. We explain the interface between the LSC Python code and TeraChem, describing how NAMD can be run in parallel over multiple GPU nodes and how we account for the phase difference in wave function that leads to sign swapping of the NAC vector. We compare, for a single NAMD trajectory, the excited state occupations for ethylene using GAMESS and TeraChem, then showcase the performance of the code for two larger test systems: solvated fulvene and the charge transfer in a zinc phthalocyanine-C 60 donor-acceptor complex.