Summary
The authors consider parametric stochastic optimization problems of the form $\min_x F(\theta,x)$ where $F(\theta,x)= \mathbb E_\xi [ f(x, \theta, \xi)]$ under the condition that $f$ is strongly convex in $x$ for any fixed $\theta, \xi$. This ensures that for fixed $\theta$, there exists a unique minimizer $x^*(\theta)$. The authors construct a sequence $x_k$ by gradient descent for $F(\theta, \cdot)$ and demonstrate that an associated sequence, denoted by $\partial_\theta x_k$, converges to the parameter derivative $\partial_\theta x^*(\theta)$ of the solution map $x^*(\theta)$. They explore various regimes of stochasticity and various learning rate schedules, obtaining rates in different settings. Strong continuity assumptions are placed on $f, \nabla f, D^2f$.
Strengths
The article is generally well-written and explains its results intuitively. The statements of results are precise, yet clear.
Weaknesses
* I find the title almost misleading. The authors do not consider derivatives (modifications) of the algorithm, but they consider the derivatives of iterates with respect to an additional parameter. I would propose something along the lines of "Derivatives in parametric optimization and their behavior along SGD iteration".
* I do not follow the numerical illustration in Section 4 at all. The independent parameter $\theta$ becomes a random quantity in this section, and it is unclear what stochastic gradient estimates the authors use: Only the deterministic objective function $F$ is specified (line 291), and $\xi$ is never mentioned in this section. I am not sure why the authors draw $\theta$ from a random distribution. I believe that the experiments either do not match the setting considered above or the presentation needs to be clarified considerably.
Questions
* Is there a conflict between the statement in line 131 that "[...] the initialization of the algorithm $x_0(\theta)$ depend[s] on some parameter $\theta$" and the assumption that $\partial_\theta x_0(\theta) =0$ in Theorem 2.2?
* In Lemma 2.1, the restriction that $\eta_k \leq \mu/L^2$ is generally much more severe than the usual bound $\eta_k \leq 1/L$. Perhaps by considering the quadratic case, could the authors speculate whether it is necessary to ensure the convergence of derivatives or whether it could be relaxed?
* In Remark 2.3, if we use a lower estimate for $\mu$, this also changes the estimate for $\kappa$: The constants $c, u$ are not independent. In fact, $\eta_0 = 1/(4\mu \kappa^2) = \mu/(4L^2)$ satisfies the necessary condition above. If we increased $c$ without adjusting $u$, we may easily enter a regime where $\eta_0> \mu/L^2$. A more careful consideration appears to be needed.