In this work, we consider an optimal control problem subject to a nonlinear PDE constraint and apply it to the regularized p-Laplace equation. To this end, a reduced unconstrained optimization problem in terms of the control variable is formulated. Based on the reduced approach, we then derive an a posteriori error representation and mesh adaptivity for multiple quantities of interest.All quantities are combined to one, and then the dual-weighted residual (DWR) method is applied to this combined functional. Furthermore, the estimator allows for balancing the discretization error and the nonlinear iteration error. These developments allow us to formulate an adaptive solution strategy, which is finally substantiated via several numerical examples. numerical example of a quasi-linear PDE constraint, and multiple goal-oriented a posteriori error estimation.In the following, we briefly refer to studies that treat parts of the three topics. Optimal control problems (specifically, a priori estimates and optimality conditions) with quasi-linear (as the p-Laplacian can be classified) elliptic PDE constraints were considered in [16,18,15]. More recently, the extension to optimal control with parabolic PDEs was discussed in [9] and [14].Optimal control problems with (single) goal functionals were investigated in [6,40,5,50,52,43]. The p-Laplacian and a posteriori error estimates were considered in [36,12,20,13], and, more specifically, for goal functional evaluations, we refer to [34,44,25]. To estimate goal functionals, we adopt the dual-weighted residual (DWR) method [7,8] in which an adjoint problem is solved to obtain (local) sensitivity measures that are used for mesh refinement. As is well-known, using a gradientbased approach for the numerical solution of optimal control problems, the same adjoint problem as for the DWR error estimator can be employed. For this reason, it is natural to combine gradient-based optimization with adjoint-based error estimation.We are specifically interested in an extended DWR version in which the discretization and (linear/nonlinear) iteration error are balanced [39,44,37]. As localization technique we employ integration by parts as done in [8] or, for residual based error estimates, in [49]. The extension of [44] to multiple goal functionals was recently undertaken in [25].Three major aims constitute the main contents of this paper: first, the design of a framework for goal-oriented error estimation for optimal control subject to a nonlinear PDE and balancing the discretization and nonlinear iteration error (Section 3). From the optimization point of view, we carefully revisit the important elements for the DWR estimator for optimization problems. The main result in this respect is the a posteriori error representation for the reduced optimal control system for an abstract problem formulation. The second aim is the extension to the simultaneous control of multiple goal functionals (Section 4). As a third goal, based on our theoretical developments, we carefully design an adaptive soluti...