In this response letter, we focus on replying theses three comments so that a proposed auxiliary-differential-equation-based (ADE) Crank-Nicolson finite-difference time-domain method with direct-splitting technique (CNDS-FDTD) can be considered as an original manuscript published in IEEE OJAP. Moreover, related applications in modeling 2D material nanostructure problems are shown. Index Terms-Auxiliary differential equation (ADE) technique, Crank-Nicolson finite-difference time-domain (CN-FDTD), directsplitting (DS), two-dimensional layered material (2DLM). A. Reply to 'No innovations in mathematical derivations' At first, we appreciate that our manuscript can attract your attention. Let us illustrate the differences between this proposal, ref. [2] and [3], so as to prove the novelty and significance of our manuscript. Dr. Sun proposed two unconditionally stable methods in ref. [2], one has the same numerical dispersion relation as the alternating-direction implicit finite-difference time-domain (FDTD) method, and the other has a much more isotropic numerical velocity. However, these two methods did not adopt truncation scheme of the complex-frequency-shifted perfectly matched layer (CFS-PML) as the absorbing boundary condition (ABC) in ref. [2], therefore, it is the key and obvious difference between ref. [2] and ours. As for ref. [3], Dr. Jiang had developed the bilineartransform-based Crank-Nicolson FDTD (CN-FDTD) with the direct-splitting (DS) scheme so that it can further reduce the requirement of the computer resources, which belongs to the performance improvement on the numerical method. Compared with ref. [3], we proposed auxiliary-differential-equation-based CNDS-FDTD method which can directly transfer the equations from frequency domain to time domain so that the transferring relation from frequency domain to Laplace domain to Z domain to time domain in ref. [3] can be avoided. Besides, Dr. Jiang et al focus on the improvement of their numerical method itself. However, we pay close attention to the combination between a novel numerical method and its more interesting applications,