While the separation of densities can be easily transferred to a separation of local potentials and energies, the treatment of non-local Coulomb contributions furthermore requires to introduce screening charges, constructed such that the corresponding multipoles sum to zero, cancelling spurious interactions between core and interatomic regions. PAW is closely related to ultrasoft pseudo potentials and is state-of-the-art in leading program packages of solid-state material science. The PAW scheme was also used to extend GPW to the all-electron Gaussian and augmented plane wave (GAPW) approach. Starting from an initial implementation for groundstate properties, GAPW was extended to time-dependent linear response perturbation theory and it was shown, that convergence of the plane-wave energy cutoff is accelerated in comparison to GPW, enabling efficiency due to smaller cutoffs while still keeping the GPW scaling for the Kohn-Sham (KS) matrix construction at O(N logN ) with system size N . More precisely, a threshold of 200 Ry was shown to be sufficient to ensure converged energies within ≈ 1µE h and bond lengths fluctuations of less then 20 µÅ. Furthermore, GAPW implementations were generalized for the treatment of core electrons, replacing pseudo potential contributions by corresponding full nuclear potential expressions, and extended to describe X-ray absorption as well as NMR and EPR spectroscopy. Regarding the computation of hyperfine parameters for the latter spectroscopies, Van Speybroeck and Waroquier et al. furthermore suggested a hybrid scheme, which implies an all-electron treatment only for selected nuclei while keeping pseudo potentials for the remaining atoms of the system. Within this contribution, we present the implementation of GAPW excited-state nuclear gradients based on a variational Lagrangian, circumventing the outlined restrictions of the recently published related GPW implementation 38 (section 2.1). Based on the fundamental ideas of Blöchl to separate the hard and soft densities (section 2.2) implying screening charges for non-local energies (section 2.3), the required working equations (section 2.4) involve a discrimination of the GAPW representation for ground-and excited-state densities.