We thank the reviewer for the review and helpful feedback. We respond to your specific concerns and questions below.
**Re: Lipschitz assumption on estimated score:** We need this assumption in the proofs so that we get the property ``when the position $\hat x$ is not too far off from the true position $x$, the estimated score is also not too far off". For instance, this property is used in Lemma A.2.
**Re: proof of Lemma A.2:** You are absolutely right and we apologize for the sloppiness in this proof. Our proof is meant to work with our Assumption 2.4 as stated, and not the stronger point-wise approximation. Specifically, what we meant to say in Line 736 is that 1) at the true position $x$, by Assumption 2.4, the expected difference of the estimated score and true score is bounded; and 2) our estimated $\hat x$ is not too far off from the true $x$, and due to lipschitz continuity assumption on both the estimated score $s$ and true score $\nabla q$, both the estimation score and the true score evaluated at the estimated position $\hat x$ is close to their value evaluated at the true position $x$. However, in this part of the proof, we got confused in the writing, and instead of applying Assumption 2.2, 2.3 and 2.4 together, we applied Assumption 2.2 and the point-wise version of Assumption 2.4. We apologize again for this mistake and have uploaded a corrected version of our submission. We mainly changed Lemma A.2 and Lemma B.2, but also parts of other theorems where these Lemmas are used, the main flow of our argument remain unchanged.
**Re: background discussion:** We added a definition of the log concave function in section C, and will continue to work on improving the clarity of our discussion for the next version of our paper.
**Re: optimality of $d$ dependence for KL and TV:** A crucial difference between [1] and our work is that [1] mainly analyze the DDPM setting with SDE, while we analyze the DDIM setting with probability flow ODE. Also, to the best of our knowledge, although the iteration complexity in [1] has a linear dependence on $d$, their optimality discussion does not rule out sublinear (or better than $\sqrt{d}$) iteration complexity. Specifically, in their calculation, the KL error induced by discretizing time is of order $\tilde{O}(d h)$, where $h$ is the step size. They mentioned that in the KL error bound, the linear $d$ dependence in this discretization error bound is optimal, however, the dependence on the step size $h$ is not necessarily optimal and [1] theorized that this could be improved to $O(d h^\beta)$ for $\beta \geq 2$, especially when there is a smoothness assumption for the score function. In this case, one can set the step size to be $h = \Theta(d^{-1/\beta})$ and obtain $d^{1/\beta}$ iteration complexity.
**Re: iteration complexity in A.10:** We are actually running predictor and corrector step for different amount of time. We are running each predictor step for $1/L$ time (specified in Algorithm 3), but we run each corrector step for $T_{corr} = \Theta(L^{-1/2}d^{-1/18})$ time (as specified in Theorem A.10). Since we take corrector step size $h_{corr} = \Theta(d^{-17/36})$, we need $\Theta(d^{5/12})$ iteration complexity to complete one corrector step, which coincides with the iteration complexity to complete one predictor step.