The character of a conical intersection, as well as the nonadiabatic dynamics of chromophores through such intersections, will be affected by a solvent environment. However, such effects are rarely considered within NAMD simulations due to the computational expense of including explicit solvent in on-the-fly dynamics. With the increased computational efficiency of Ter-aChem for large-scale electronic structure calculations, solvated chromophore NAMD, although still computationally formidable, becomes feasible within PySCES. We next show the performance of PySCES employing TeraChem as we increase the amount of explicit solvent included in the electronic structure calculations.