Augmented plane wave methods enable an efficient description of atom-centered or localized features of the electronic density, circumventing high energy cutoffs and thus prohibitive computational costs of pure plane wave formulations. To complement existing implementations for ground-state properties and excitation energies, we present the extension of the Gaussian and augmented plane wave method to excited-state nuclear gradients within the Tamm–Dancoff approximation of time-dependent density functional theory and its implementation in the CP2K program package. Benchmarks for a test set of 35 small molecules demonstrate that maximum errors in the nuclear forces for excited states of singlet and triplet spin multiplicity are smaller than 0.1 eV/Å. The method is furthermore applied to the calculation of the zero-phonon line of defective hexagonal boron nitride. This spectral feature is reproduced with an error of 0.6 eV in comparison to GW–Bethe–Salpeter reference computations and 0.4 eV in comparison to experimental measurements. Accuracy assessments and applications thus demonstrate the potential use of the outlined developments for large-scale applications on excited-state properties of extended systems. The data in this record are supplementary data to the published article with DOI 10.1021/acs.jctc.4c00614.