Abstract: An efficient time-domain trapezoidal difference based model order reduction(MOR) method is proposed to further improve the efficiency and accuracy of the existing methods for interconnect circuits.A recursive relation is derived from the trapezoidal difference formulation of the time-domain equation of interconnect circuits.The recursive relation formulates a non-homogenous Krylov subspace.A non-homogeneous Arnoldi method is employed to construct the orthonormal basis of this non-homogenous Krylov subspace.The orthonormal basis is then used to project the original interconnect circuits to reduced-order models.The proposed method can guarantee the value matching of the state vector after time-domain trapezoidal difference for the original interconnect circuits and the reduced-order models.The accuracy of the reduced-order models can thus be guaranteed.The MOR procedure is numerically stable and the reduced-order models also preserve the passivity.In comparison with the existing time-domain methods,the computational cost of the proposed method is significantly reduced.Compared with the existing frequency-domain methods,the proposed method can achieve higher accuracy.